Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,449,200 G→A 5.9% intergenic (+43/+28) 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,2000GA 5.9% 256.8 / 3.9 85intergenic (+43/+28)smrB/yfcOputative DNA endonuclease/DUF2544 family putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (41/39);  new base A (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 = 5.88e-01

ATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCG  >  NC_000913/2449051‑2449349
                                                                                                                                                     |                                                                                                                                                     
aTCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCg                                                                                                                                                       >  2:450653/1‑150 (MQ=255)
  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:1161476/1‑150 (MQ=255)
                                                                                      cTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATActgc                                                                 <  1:211182/150‑1 (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                                                  >  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                      <  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               <  1:815747/150‑1 (MQ=255)
                                                                                                                                        tAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAggg               <  2:298738/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)
                                                                                                                                                     |                                                                                                                                                     
ATCCGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCG  >  NC_000913/2449051‑2449349

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: