Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,616,073 T→C E331E (GAA→GAG)  yjjI ← DUF3029 domain‑containing protein YjjI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,616,0730TC95.9% 142.8 / NA 49E331E (GAA→GAG) yjjIDUF3029 domain‑containing protein YjjI
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (22/25);  minor base A (2/0);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

ATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTG‑CCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCTACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAA‑CGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTGCGCTCGGCAATGGCTTTCAGGTTAAGGCGTACC  >  NC_000913/4615830‑4616319
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
atgatgCGGGGCGCCAGTTTGCAGATGGGTGTCCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGGGTGGTG‑CCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTcacaca                                                                                                                                                                                                                                                  <  2:419847/251‑1 (MQ=255)
   atgCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTG‑CCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGc                                                                                                                                                                                                                                                >  1:127977/1‑250 (MQ=255)
               attttGCAGATGGGTGACCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTG‑CCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCtt                                                                                                                                                                                                                                   <  1:499346/249‑1 (MQ=255)
                            gTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTG‑CCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTCTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGACAGTTAACCGCTTCCGCCAGCCCAt                                                                                                                                                                                                                       >  1:299156/1‑250 (MQ=255)
                             tGATCGGATCTGGCTCATCGCCATACGGCACCCGCGCGCCCGGCGTGGGG‑CCGGTCTCGGAACTGATCACCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACCAACGCCGCCAGTGGCGCGCTGATGAGATAACCCACTTCGTTTGCGTCGGACTCTTTACCGTAAGGCGCGGCAATCCCTTCTTTCTCACACACCACGTTAATCGCTTCCGCCacaccaccc                                                                                                                                                                                                                     >  2:66501/1‑243 (MQ=255)
                                  ggATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTG‑CCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGcc                                                                                                                                                                                                                <  1:504940/251‑1 (MQ=255)
                                         gCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTG‑CCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGc                                                                                                                                                                                                                                                >  1:357381/1‑212 (MQ=255)
                                           tCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTG‑CCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTgg                                                                                                                                                                                                       >  1:75004/1‑251 (MQ=255)
                                            cATCGCCAGACGGCAAACGCGCGCCCGGCGGGGTG‑CCGATATCGGAACTGATCCCCGACGGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGc                                                                                                                                                                                                      <  2:131344/251‑1 (MQ=255)
                                              ccgcctgtcGGCAAACGCGCGCGCGGC‑TGGTGGCCGAAATCGGACTTGTCCCCCGTCTGGGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCTCGGGGTTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCGTTCTCACACAGCAAGTTAAACGCTTCCGCCAGCCCATACATGCCAAACATTGGc                                                                                                                                                                                                      <  2:49859/241‑1 (MQ=255)
                                               cGCCATACGGCAAACGCGCGCCCGGCGTGGTG‑CCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTAATTTGCGGCGGCCTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCGCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGACAAACCTTGGccc                                                                                                                                                                                                    >  2:466398/1‑248 (MQ=255)
                                                          aaaCGCGCGCCCGGGGGGGTG‑CCGATACGGGAACGGTGCCCCGCGTGGGCGTGTAACAGGGGGCGTTGTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAA‑CGTTCAgg                                                                                                                                                                                       <  2:382404/251‑1 (MQ=255)
                                                                  gcCCGGCGTGGTG‑CCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGGAGCGCGCGGCCATCCCTTCTTTCTAACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACCTTGGCCCACAA‑CGTTCAggggtgatc                                                                                                                                                                                >  2:8638/1‑244 (MQ=255)
                                                                   cccGGCGTGGTG‑CCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAA‑CGTTCAGGGTTAATc                                                                                                                                                                                >  1:499706/1‑249 (MQ=255)
                                                                              g‑CCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAA‑CGTTCAGGGTTAAGCAGCCCTTCTTTCa                                                                                                                                                                   <  1:388946/251‑1 (MQ=255)
                                                                                      tCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAACCTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGAGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGATTCCGCCAGCCCATACATGCCAACCATTGGCACCAAC‑AGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGa                                                                                                                                                             >  2:216113/1‑250 (MQ=255)
                                                                                        ggAACTGCTCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCCGCCATATTTCACGGGGGTATTGGCGACACACTCCGCCAGTTGCGCGCTGATGCGATAACCCAATTCATTTGCGGCGGCTTCTTTCACGTAGCGCGCGGCAATCCCTTCTTTCTCACACCGCAAGTTTACCGCTTCCGCCAGCCCATACCTGCCAAACATTCGAACAAAA‑CGTTCAGGGTTAATCCGCCCTTCTTTCAcccgggcg                                                                                                                                                           >  2:487102/1‑244 (MQ=255)
                                                                                                      actgggcgggaaaCATGGCGGGTTTTTGCCAGCCATATGTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACGTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAA‑CGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAaagaag                                                                                                                                            <  2:75004/240‑1 (MQ=255)
                                                                                                            cGTGTAACATGGCTCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAA‑CGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGAt                                                                                                                                      <  1:293252/251‑1 (MQ=255)
                                                                                                            cGTGAAACCTGGCGCGTTTTTGCCCGCAATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGACTCCGCCAGCCCATCCATGCCAAAAATTGG‑ACAAAAACGTTCCGGGTTAATCAGCCCTTCTTTCAACAGGAAGCTATTCCCACAGAAGTGGGAt                                                                                                                                      >  2:432551/1‑251 (MQ=255)
                                                                                                               gtaaaatgggtggttggtGCCAGCCATATTCCAAGGGGTTCTGGGCGACAACCGCCGTCAGGTGCGCTCTGATGCGATAACCCGCTGTATTCGCGGCGACTTCTTTACCGTAGCGCGCGGCAATCCCTTCGCTCTCACACAGCAAGTTAACCGCTTCCGCCATCCCATACATGCCAAACATTGGCACAAAA‑CGATCAGGGTTAATCAGCCCTTCTCTCATCAGGAAGCTATTCTCAAAGAAGTGTGAttg                                                                                                                                    <  2:298935/233‑1 (MQ=255)
                                                                                                                 aaCATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGGATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACCTGCCAAACCTTGGCACAAAA‑CGTTCAGGGGTAATCCGCCCTTCGTTCACCAGGATGCTATTCTCAACGAAGGGTGAttgttt                                                                                                                                 >  2:346813/1‑250 (MQ=255)
                                                                                                                   aaTGGCGTGTGTTGGCCAGCCATAGTTCACGGGGGTAATGGCGACTAACTCCGTCAGTTGCGCGCTGATGCGATAACCCACGTCAGTTGCGGCGGCTTCTGTACCGCAGCGCGCGGCAATCCCTTCTTTCTCACCCAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCAAAAAA‑CGTTCAGGGTTAATCAGCCCTTCTTGCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGAt                                                                                                                               <  2:129278/250‑1 (MQ=255)
                                                                                                                            tgttcgccagccaGATTTCCCGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAACCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAA‑CGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAAt                                                                                                                       <  2:31070/245‑1 (MQ=255)
                                                                                                                              tgtGCCAGCCAATTTTCACGGGGGTTTTGGCGACAAACGCCGCCAGTTGCGCTCTTATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAA‑CGTGCAGGGTTAATCAGCCCTTCTTGCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTc                                                                                                                     <  2:241176/248‑1 (MQ=255)
                                                                                                                                 gccagccaTATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAA‑CGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCacac                                                                                                                 <  1:500677/251‑1 (MQ=255)
                                                                                                                                  ccagccaTATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATACGATAACCCACTTCATACGCGGCGGCTAATTTACCGTAGCGCGCGGCAATCCCTTCTTTATCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACCTTGGCACAAAA‑CGTTCAGGGTTAATCAGACCTTCTTTCACCAGGAAGCTATTCTCACAGAAGTGTGATTGTTGATAGAGGAATTCACAcc                                                                                                                >  1:451613/1‑251 (MQ=255)
                                                                                                                                          gtttCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCACTGATGCGCTAACCCACTTCATTTGCGGCGGCTTCTTTACCCTAGCGCGCGGCAATCCCTTCTTTCTCACACATCAAGTTAACCGCTTCCGCCAACCCATAAATGCCCAACATTGGCACAATA‑CGTTCAGGGATAATCAGCCCTTCTTTCACCATGAAGCTATTCTCACGGAAGTATGAATGTTGGTAGAGGAATTCACACCGCGCATCg                                                                                                        >  2:375780/2‑251 (MQ=255)
                                                                                                                                               aCGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAA‑CGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCgatgat                                                                                                   >  1:192309/1‑251 (MQ=255)
                                                                                                                                                cGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAA‑CGTTCAGGGTTAATCAGCCCTTCGTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGatgatg                                                                                                  <  1:487102/251‑1 (MQ=255)
                                                                                                                                                  ggggTAGTGGTGACAAACTCCGCCAGTTGCGCGCGGATTCGATAACCCACTTCATTTGCGGCGGCTGCTTTACCGTAGCGCGCGTCAATCCCTACTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACAGGCCAAACATTGGCACAAAA‑CGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTTATAGAGGAATTCACCCCGCGCATCGATGATGGc                                                                                                <  2:470704/251‑1 (MQ=255)
                                                                                                                                                      tATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGAGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCATCCGCCAGCCCATACATGCCACACATTGGCACAAAA‑CGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGATATTCTCAAAGACGTGTGATTGTTTATAGAGGAATTCACACCTCGCATCGATGATGGCGATc                                                                                            >  2:414426/1‑251 (MQ=255)
                                                                                                                                                                  acaccGCCAGTTGCGCGCTGCTGCGATAACCCACTTCATTTGCGGCGTCTTCTTTACCGTAGCGCGCGGCAAGCCCTTCTGTCTCACACAGCAAGGTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAA‑CGTTCAGGGGTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCag                                                                                <  2:62132/248‑1 (MQ=255)
                                                                                                                                                                   cTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAA‑CGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCagg                                                                               >  2:170417/1‑250 (MQ=255)
                                                                                                                                                                      cGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCACGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAA‑CGCTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCACAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGCTGATGGCGATCTGCTGCTGGCagtagt                                                                            >  1:72414/1‑251 (MQ=255)
                                                                                                                                                                       gcccgTTGCGCGCGTATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAA‑CGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTg                                                                           <  2:127977/247‑1 (MQ=255)
                                                                                                                                                                                  ggtGATGCGATAACCCACTTCATGTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAA‑CGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAAGTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTgc                                                                <  2:122228/249‑1 (MQ=255)
                                                                                                                                                                                             aaCCCACTTCATTTGCGGCGTCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAA‑CGTGCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTAGTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAaagaag                                                     <  1:472943/251‑1 (MQ=255)
                                                                                                                                                                                                   cTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACCTGCCAAACATTGGCACAAAC‑CGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAAAGTGCGAGTAAAGAAGTCAtac                                               >  1:63652/1‑249 (MQ=255)
                                                                                                                                                                                                     tCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAA‑CGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATcc                                               >  1:13469/1‑249 (MQ=255)
                                                                                                                                                                                                            cggcggCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAA‑CGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAg                                             <  2:35176/244‑1 (MQ=255)
                                                                                                                                                                                                            cggcggCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAA‑CGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAg                                             >  1:35176/1‑244 (MQ=255)
                                                                                                                                                                                                             ggcggcTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAA‑CGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGc                                     <  1:171448/251‑1 (MQ=255)
                                                                                                                                                                                                                ggcTTCTTTACCGTAGCGCGCGGCAAGCCCTTCTTTCTCTCACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAT‑CGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCGGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTgc                                  <  1:293828/251‑1 (MQ=255)
                                                                                                                                                                                                                         aCCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAA‑CGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTGCGCTCGGCa                          >  1:219946/1‑250 (MQ=255)
                                                                                                                                                                                                                             tAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCGTCCGCCAGCCCATACATGCCAAACATTGGCACAAAA‑CGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACAACGCGCATCGATGATGGCGATCTGCTGCTGGCAGGCGCGCGGGAGAGTGCGCGTAAAGAAGTCGTTCAGCGCTTCGCTGGGCTCGGCACTGGc                     >  2:344493/1‑251 (MQ=255)
                                                                                                                                                                                                                                cgcgcgCCAATCCCTTCTCTCCCCCACAGCAAGTTAACCGCGTCCGCCAGCCCATACATGCCAAACATTGGCACAAAA‑CGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCTCGCATCGATGATGGCGATCTGCTGCTGGCAGGAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTGCGCTCGGCAATGGCttg                  <  2:415085/251‑2 (MQ=255)
                                                                                                                                                                                                                                               tcttTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAA‑CGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGATGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTGCGCTCGGCAATGGCTTTCAGGTTAAGGCGt     <  1:466398/249‑1 (MQ=255)
                                                                                                                                                                                                                                                cttTATCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAA‑CGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTGCGCTCGGCAATGGCTTTCAGGTTAAGGCGTAcc  >  2:184241/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
ATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTG‑CCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCTACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAA‑CGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTGCGCTCGGCAATGGCTTTCAGGTTAAGGCGTACC  >  NC_000913/4615830‑4616319

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑