Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,616,022 T→C V348V (GTA→GTG)  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,0220TC98.1% 154.4 / ‑4.4 52V348V (GTA→GTG) yjjIDUF3029 domain‑containing protein YjjI
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (21/30);  total (21/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.37e-01

TGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTG‑CCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCTACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAA‑CGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAG  >  NC_000913/4615778‑4616268
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       
tgagggTTTCGGCGAGCGTCAGATTGGCGCTGATGCCGAAATAATAATAAGCATGAGGCGGGGCGACAGTTTGCAGATGGGGGTCCGGAGCTGGCTCATCGCCTTACGGCAAACGCGCGCCCGGGGTGGTG‑CCGATAGCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTc                                                                                                                                                                                                                                                    <  2:437161/246‑1 (MQ=255)
tGATGGTTTCGTCGAGCGTCAGAAGGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTG‑CCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTc                                                                                                                                                                                                                                                    <  2:120568/250‑1 (MQ=255)
        tcgccgAGCGTCAGAATGCCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGAGCGGAGCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTG‑CCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCTACTTCATTTgcgg                                                                                                                                                                                                                                            <  2:128321/246‑1 (MQ=255)
          gtcgAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTG‑CCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGcggcgg                                                                                                                                                                                                                                         >  1:407916/1‑251 (MQ=255)
             gAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTG‑CCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAAACCACTTCATTTGcggcggatt                                                                                                                                                                                                                                      >  1:62132/1‑248 (MQ=255)
               gCGTCAGAATGTCGCTGATGCCGGAATAATAAAAAGCATGAGGCGGGGCGACAGTTTGCAGATGGGTGATCGGAGCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGGGTGGTG‑CCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACTGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGcttct                                                                                                                                                                                                                                    <  1:474209/251‑1 (MQ=255)
               gCGTCAGAATGGCGGTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTGGCAGATGGGTTATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTG‑CCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGcttct                                                                                                                                                                                                                                    <  2:446736/251‑1 (MQ=255)
                         gtccgttaTGCCGGAATAAGAATAAGCATGATGCGGGGCGACAGGTTGCAGATGGGTGAGCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCTTGGGG‑CCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAgc                                                                                                                                                                                                                          <  2:382439/244‑1 (MQ=255)
                         gTCGCTGATGCCGGAATAATAAGAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTG‑CCGATAGCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAg                                                                                                                                                                                                                           <  1:232039/250‑1 (MQ=255)
                                             aataaGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATTTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTG‑CCGATATCGGAACTGATCCCCGACTGTGCGTTTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTATCGCGCGGCAATCCCTTCttt                                                                                                                                                                                                       <  1:441163/250‑1 (MQ=255)
                                             aataaGCATGAGGCGGGGCGACAGTTTGCAGATGGGTGATCGGAGCGGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTG‑CCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCtttc                                                                                                                                                                                                      <  2:234820/251‑2 (MQ=255)
                                                    atgatgCGGGGCGCCAGTTTGCAGATGGGTGTCCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGGGTGGTG‑CCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTcacaca                                                                                                                                                                                               <  2:419847/251‑1 (MQ=255)
                                                    atgatgAGGGGCGACAGTTCGCAGATTGGTGCACGGATCTGGCTCTTCGCCATACGGCAAACGCGCGCCCGGCGGGGTG‑CCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCcttctctctcccaca                                                                                                                                                                                               <  2:62386/251‑11 (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)
                                                                                                                                                                            cgcgTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCcttct                                                                                                                                                                                                         >  1:119964/1‑122 (MQ=255)
                                                                                                                                                                            cgcgTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCcttct                                                                                                                                                                                                         <  2:119964/122‑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)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       
TGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTG‑CCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCTACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAA‑CGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAG  >  NC_000913/4615778‑4616268

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