Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,449,202 A→T 5.9% intergenic (+45/+26) smrB → / ← yfcO putative DNA endonuclease/DUF2544 family putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,449,2020AT 5.9% 244.1 / 4.7 85intergenic (+45/+26)smrB/yfcOputative DNA endonuclease/DUF2544 family putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (41/39);  new base T (3/2);  total (44/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.91e-01

CCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCGTG  >  NC_000913/2449053‑2449351
                                                                                                                                                     |                                                                                                                                                     
ccGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGaa                                                                                                                                                       <  1:19935/150‑1 (MQ=255)
  gCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTTCGGGTAATGGCAGGGGCGAATg                                                                                                                                                     >  2:1102411/1‑150 (MQ=255)
      gtgATGGCCTTTCATCAGTCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGcccc                                                                                                                                                 <  1:380169/150‑1 (MQ=255)
       tgATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCcct                                                                                                                                                >  2:651684/1‑150 (MQ=255)
           ggCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGCGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCATTCGCCcctgcca                                                                                                                                            <  2:1706288/150‑2 (MQ=255)
            gCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCATTCGCCcctgccat                                                                                                                                           <  2:1439526/150‑3 (MQ=255)
               tttCATCAGGCACCAAATGAGTATGGCGGTGATGCGGCGTTGGTGGTGTTGATTGAAGTCGAAGAGTGGCCGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGgtg                                                                                                                                        <  2:878820/150‑1 (MQ=255)
               tttCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGgtg                                                                                                                                        >  2:1373609/1‑150 (MQ=255)
               tttCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGGTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGgtg                                                                                                                                        <  1:1391844/150‑1 (MQ=255)
                ttCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGAGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGgtgt                                                                                                                                       >  1:1206605/1‑150 (MQ=255)
                   cacaGGCACCAAACGAGTATGGCGGTGATGCGGCGCGGGTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGGTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATc                                                                                                                                    <  1:1210783/148‑1 (MQ=255)
                   atcaGGCACCAAAAGAGTATGGCGGTGATGCGGCGGTGTTGGTGTTGATTGAAGTCGAAGTGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATc                                                                                                                                    <  2:660005/150‑1 (MQ=255)
                    tcaGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCt                                                                                                                                   >  1:323266/1‑150 (MQ=255)
                     caGGCACCAAAAGAGGATGGCGGTGATGCGGCGTTGTTGGTGTTGAGTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCAATTAGGAGTTGAGCGCGTTGCGCCAGGTACGGGGAATGGCAGGGGCGAATGACCCTGCCTTgtgtttggt                                                                                                                                  >  2:1044729/1‑145 (MQ=255)
                          aCCAAAAGAGTATCGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGc                                                                                                                             <  1:1157092/150‑1 (MQ=255)
                            cAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTa                                                                                                                           >  2:726110/1‑150 (MQ=255)
                             aaaaGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTAc                                                                                                                          <  2:766628/150‑1 (MQ=255)
                                 gagTAGGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGACGAGGGGCTGCCGCCGGAGTGGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACt                                                                                                                      <  1:1696189/150‑1 (MQ=255)
                                   gTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTTCGGGTAATGGCAGGGGCGAATGCCACTGCCTGGTGTATCTAAGGGCTACAACTCa                                                                                                                    >  1:339040/1‑150 (MQ=255)
                                   gTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTTGTGTATCTAAGGGCTACAACTCa                                                                                                                    >  2:302842/1‑150 (MQ=255)
                                   gTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCa                                                                                                                    >  2:1474133/1‑150 (MQ=255)
                                   gTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAACAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGAGAGTTGCGACAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGGGTATCTTAGGGCTACAACTCa                                                                                                                    >  1:257488/1‑150 (MQ=255)
                                     aTGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCTGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATTGCCCTGCCTGGTGTTTATAAGGGGTTCaaaagctt                                                                                                                  >  1:1530922/1‑144 (MQ=255)
                                     aTGGCGGTGATGCGGCGTGGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCAtt                                                                                                                  <  1:1437306/150‑1 (MQ=255)
                                      tGGCGGTGATGAGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCAtta                                                                                                                 >  2:1695773/1‑149 (MQ=255)
                                       ggCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCa                                                                                                                >  1:796237/1‑150 (MQ=255)
                                        gCGGTTAGGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGATGCCCTGAGGAGTTGAGCGCGGTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAg                                                                                                               <  1:1137056/150‑1 (MQ=255)
                                                       ggttgGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCAGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCg                                                                                                <  1:1082637/149‑1 (MQ=255)
                                                         ttgGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGc                                                                                              >  1:1364096/1‑150 (MQ=255)
                                                          tgGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCa                                                                                             <  1:594951/150‑1 (MQ=255)
                                                               ttgattgaAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAg                                                                                        <  2:568009/150‑1 (MQ=255)
                                                                   ttgaAGTCGCAGAGTGGCTGCCGCCGGAGTGGCCCTGAGGAGGTTAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATc                                                                                    <  2:578527/150‑1 (MQ=255)
                                                                   tggaAGTCGAAGTGGGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGAAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATc                                                                                    <  2:1136530/148‑1 (MQ=255)
                                                                     gaAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATAGa                                                                                  >  1:347253/1‑150 (MQ=255)
                                                                       aggcGAAGAGGGGCTGCAGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGAtt                                                                                <  2:1240172/147‑1 (MQ=255)
                                                                        gTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTg                                                                               >  1:923087/1‑150 (MQ=255)
                                                                          cGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGc                                                                             >  2:410485/1‑150 (MQ=255)
                                                                           gAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAAAGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCCAGTATAATCGATTGGcg                                                                            >  2:1558119/1‑149 (MQ=255)
                                                                           gAAGAGGGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCt                                                                            <  1:1172580/150‑1 (MQ=255)
                                                                             agagTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCCGGGGCGAATGCCCCTGCCTGGTGTATCGAAGGGCTTCAACTCATTCAGTTTGTACTGGAACCGGCAAGTAGCATCGATTGGCTGc                                                                          >  1:555156/1‑150 (MQ=255)
                                                                              gtgTGGCTGCCTGCGGAGGTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCa                                                                         <  1:323536/148‑1 (MQ=255)
                                                                                    cTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATActgc                                                                   <  1:211182/150‑1 (MQ=255)
                                                                                    cTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATActgc                                                                   >  1:1161476/1‑150 (MQ=255)
                                                                                     tcccgccgGAGTTTCCATGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATTGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGTGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATActgct                                                                  <  2:1096552/148‑1 (MQ=255)
                                                                                       ccgccgGTGGTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTgg                                                                <  1:222831/150‑1 (MQ=255)
                                                                                         gccgGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTa                                                              >  2:16148/1‑150 (MQ=255)
                                                                                         gccgGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTa                                                              <  1:450653/150‑1 (MQ=255)
                                                                                             gAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCaa                                                          <  2:323266/150‑1 (MQ=255)
                                                                                              aGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAac                                                         >  2:282498/1‑150 (MQ=255)
                                                                                                  gCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgc                                                     >  1:1189891/1‑150 (MQ=255)
                                                                                                    ccTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgcac                                                   >  2:809512/1‑150 (MQ=255)
                                                                                                     cggaggagGGGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGGATCTAACGGCGACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgcacg                                                  <  1:839750/148‑1 (MQ=255)
                                                                                                     cTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgcacg                                                  >  2:1195592/1‑150 (MQ=255)
                                                                                                      tGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgca                                                    <  2:428443/147‑1 (MQ=255)
                                                                                                      tGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgca                                                    >  1:428443/1‑147 (MQ=255)
                                                                                                      tGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACAGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgcacgc                                                 >  1:1234450/1‑149 (MQ=255)
                                                                                                              tGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATCCCCCTGCCTGGTGTATCTAAGGGCTACATCTCATTCAGTATGTACTGGATCCGGCAAGTACAAACGATTGGCTGCATACTGCAGGTATCAACACGCACTACCATATCg                                         >  1:270423/1‑150 (MQ=255)
                                                                                                               gAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGCCCATATCgt                                        >  1:969397/1‑150 (MQ=255)
                                                                                                                 gcgcGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTCCAACTCAGTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCgtgt                                      >  2:1045968/1‑150 (MQ=255)
                                                                                                                   gcgTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCAAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGCTTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTaa                                    >  2:391433/1‑150 (MQ=255)
                                                                                                                    cgTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAAc                                   >  2:35563/1‑150 (MQ=255)
                                                                                                                     gTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAAGTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACt                                  >  2:826347/1‑150 (MQ=255)
                                                                                                                     gTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACt                                  <  1:15359/150‑1 (MQ=255)
                                                                                                                      ttGCGCCAGGTACGGGTAATGGCAGGGGCGAAGCCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATACAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTg                                 <  1:882262/150‑1 (MQ=255)
                                                                                                                        gcgcCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTTGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGtt                               <  1:302842/150‑1 (MQ=255)
                                                                                                                            cAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTc                           >  2:1595944/1‑150 (MQ=255)
                                                                                                                               gTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGt                        <  1:1606988/150‑1 (MQ=255)
                                                                                                                               gTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGt                        >  2:682145/1‑150 (MQ=255)
                                                                                                                                tACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGtt                       <  1:1383849/150‑1 (MQ=255)
                                                                                                                                  cGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTAATGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATAAACACGCACGACCATATCGTGTAACTGTTTGTCCGTtat                     >  2:30803/1‑150 (MQ=255)
                                                                                                                                  cGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCGAAGGGCGACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGCTTGGGTGCATACTGCTGGTTTAAACACGCCCGACCATATCGTGGTACTGTTTGTCCGTgat                     >  2:721788/1‑150 (MQ=255)
                                                                                                                                   gggTAATTGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTtata                    <  2:62021/150‑1 (MQ=255)
                                                                                                                                    ggTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGCCAAGTAGAAACGATTGGCTGCATACTGCTGGTATCAACACGCACTACCATATCGTGTAACTGTTTGTCCGTTATAg                   >  2:58626/1‑150 (MQ=255)
                                                                                                                                    ggTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAAACGGCAAGTAGAATCGAGTGGCTACATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTTTAg                   >  2:953456/1‑150 (MQ=255)
                                                                                                                                     gTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAgg                  <  2:1364096/150‑1 (MQ=255)
                                                                                                                                     caccaggCAGGGGCATTCGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAgg                  >  2:960104/6‑150 (MQ=255)
                                                                                                                                      tAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTTAGGGCTACAACTCATTCAGTATGTACTGGACCCGGCAAGTAGAATCGATTGGCTGCTTACTGCGGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAggg                 <  2:260468/150‑1 (MQ=255)
                                                                                                                                      tAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAggg                 <  2:298738/150‑1 (MQ=255)
                                                                                                                                      tAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAggg                 <  1:815747/150‑1 (MQ=255)
                                                                                                                                       aaTGGCAGGGGCGAATTCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATAATGCTGGTATCAACAAGCACGACCATATCGTGTTACTTTTTGTCCGTTATAGGgt                >  2:1146675/1‑150 (MQ=255)
                                                                                                                                        aTGGCAGGGGTGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGgtg               <  1:472023/150‑1 (MQ=255)
                                                                                                                                          ggCAGGGGCATTCGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGgtgta             >  1:228057/1‑149 (MQ=255)
                                                                                                                                          ggCAGGGGCATTCGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTc             >  1:590357/1‑150 (MQ=255)
                                                                                                                                           gCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCg            <  2:1134323/150‑1 (MQ=255)
                                                                                                                                                  cGAATGCCCCTTCCTGGTGTATCTAAGGGCTACAACTCATTCTGTTTGTACTGGAACCTGCATGTTGATTGGATTGGCTGCATAGTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGc     <  2:1417120/150‑1 (MQ=255)
                                                                                                                                                   gAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCg    <  1:649237/150‑1 (MQ=255)
                                                                                                                                                    aaTGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCGt   >  1:239050/1‑150 (MQ=255)
                                                                                                                                                     aTGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCGTg  <  1:1573004/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCGTG  >  NC_000913/2449053‑2449351

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: