Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,449,201 A→T 5.9% intergenic (+44/+27) 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,2010AT 5.9% 236.7 / 4.6 85intergenic (+44/+27)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 = 6.86e-01

CCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCGT  >  NC_000913/2449053‑2449350
                                                                                                                                                    |                                                                                                                                                     
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                                                             <  1:450653/150‑1 (MQ=255)
                                                                                         gccgGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTa                                                             >  2:16148/1‑150 (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                                                   >  1:428443/1‑147 (MQ=255)
                                                                                                      tGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgca                                                   <  2:428443/147‑1 (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                       >  2:682145/1‑150 (MQ=255)
                                                                                                                               gTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGt                       <  1:1606988/150‑1 (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)
                                                                                                                                                    |                                                                                                                                                     
CCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCGT  >  NC_000913/2449053‑2449350

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: