Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,414,240 A→T 6.6% intergenic (+49/‑35) yjfC → / → aidB ATP‑Grasp family ATPase/DNA alkylation damage repair protein; flavin‑containing DNA binding protein, weak isovaleryl CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,414,2400AT 6.6% 479.2 / 15.2 164intergenic (+49/‑35)yjfC/aidBATP‑Grasp family ATPase/DNA alkylation damage repair protein; flavin‑containing DNA binding protein, weak isovaleryl CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (83/69);  new base T (6/5);  total (90/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 167 total aligned reads displayed.

GGATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGG  >  NC_000913/4414091‑4414389
                                                                                                                                                      |                                                                                                                                                     
ggATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATa                                                                                                                                                       >  2:1454512/1‑150 (MQ=255)
 gATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATAc                                                                                                                                                      >  2:1730514/1‑150 (MQ=255)
     gTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCAAACTAAt                                                                                                                                                  >  2:887811/1‑150 (MQ=255)
      tCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATg                                                                                                                                                 >  1:401236/1‑150 (MQ=255)
      tCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATg                                                                                                                                                 >  1:1283654/1‑150 (MQ=255)
       cGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTGTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATgg                                                                                                                                                >  2:1295415/1‑150 (MQ=255)
         aggatgtATCGTGTGGTATGGGGAGCCGTGAAGATAACCAACTGATGACCAAAGACACCTCACGGTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTg                                                                                                                                              <  1:1575267/148‑1 (MQ=255)
               aaGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGAGACTGca                                                                                                                                        >  1:1588306/1‑149 (MQ=255)
                  cGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCAtt                                                                                                                                     <  1:874194/150‑1 (MQ=255)
                  cGTGCGGAATGGGGAGCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTTGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTCACCATTAGTATGGTCACTGCCAtg                                                                                                                                     <  2:1664973/150‑2 (MQ=255)
                    tGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGa                                                                                                                                   >  2:1485384/1‑150 (MQ=255)
                     gCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACAATACTAATGGTGACTGCCATTGAt                                                                                                                                  >  1:1148339/1‑150 (MQ=255)
                            ggggATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTT‑GCAATCTCTTTCTGTCATGAATCAATGGCAGTGACCATACTAATGGTGACTGCCATTGATggaggga                                                                                                                           >  2:319693/1‑150 (MQ=255)
                            ggggATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATggaggga                                                                                                                           >  2:377948/1‑150 (MQ=255)
                             gggATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGCCCATACTAATGGTGACTGCCATTGATggagggag                                                                                                                          >  2:32063/1‑150 (MQ=255)
                                   cGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAg                                                                                                                    >  1:1188610/1‑150 (MQ=255)
                                    gTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGTATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGt                                                                                                                   <  1:1065086/150‑1 (MQ=255)
                                      gAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGc                                                                                                                 >  2:618796/1‑150 (MQ=255)
                                       aaGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCa                                                                                                                >  2:106134/1‑150 (MQ=255)
                                          ctAACCCACTGCTCCCCAAAGCCACCTCACCTTTCGGTCCGCAGTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACAACTGGGGACTGCCATTGATGGAGGGAGACACAGTGCACTg                                                                                                             <  2:1694906/149‑1 (MQ=255)
                                           tAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATggagggag                                                                                                                          <  2:460779/136‑1 (MQ=255)
                                           tAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATggagggag                                                                                                                          >  1:460779/1‑136 (MQ=255)
                                           tAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTgg                                                                                                            <  2:691285/150‑1 (MQ=255)
                                            aaCACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGc                                                                                                           >  2:1425791/1‑150 (MQ=255)
                                                    gATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTcac                                                                                                   >  2:464603/1‑150 (MQ=255)
                                                     aTCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTcaca                                                                                                  <  2:1014743/150‑1 (MQ=255)
                                                       cACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACAcc                                                                                                >  2:1300158/1‑150 (MQ=255)
                                                        aCCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCCCTGGCAAACTCGCACCg                                                                                               >  1:1125080/1‑150 (MQ=255)
                                                        aCCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCg                                                                                               >  2:1537028/1‑150 (MQ=255)
                                                         ccAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGt                                                                                              >  2:636821/1‑150 (MQ=255)
                                                          cAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGtt                                                                                             >  2:814970/1‑150 (MQ=255)
                                                            aaGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGtttt                                                                                           <  2:1591598/150‑1 (MQ=255)
                                                              gACACCTCACGTTTCGGTCCGCAGTACATTGATTGATAAGAA‑GTGATACGCAATATCCTTCTGTCATGAATCCATTGCAGTGACCATACTAATGGTGAATTGCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTa                                                                                         <  1:84388/150‑1 (MQ=255)
                                                                 acCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATc                                                                                      <  1:1635479/150‑1 (MQ=255)
                                                                  ccTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCa                                                                                     >  1:646960/1‑150 (MQ=255)
                                                                  ccTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCACTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCa                                                                                     >  1:1026069/1‑150 (MQ=255)
                                                                     cACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAAcc                                                                                  >  1:197189/1‑150 (MQ=255)
                                                                     cACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAAcc                                                                                  <  2:565728/150‑1 (MQ=255)
                                                                         tttCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCtata                                                                              <  1:284988/150‑1 (MQ=255)
                                                                         tttCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCtata                                                                              <  2:126223/150‑1 (MQ=255)
                                                                          ttCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATAc                                                                             >  1:1212454/1‑150 (MQ=255)
                                                                          ttCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATAc                                                                             >  1:563663/1‑150 (MQ=255)
                                                                           tCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATAcc                                                                            <  1:1389528/150‑1 (MQ=255)
                                                                             gTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCAt                                                                          >  2:95725/1‑150 (MQ=255)
                                                                             gTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCAt                                                                          >  1:554868/1‑150 (MQ=255)
                                                                                ccGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTaa                                                                       >  1:1236823/1‑150 (MQ=255)
                                                                                  gCATTACATTGCTGGATTAGAATTTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATAATAATTGTGACTGACATTGATGGAGGTATACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAAt                                                                     <  1:260972/150‑1 (MQ=255)
                                                                                  gCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAAt                                                                     >  2:1707084/1‑150 (MQ=255)
                                                                                   cATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATa                                                                    <  1:713483/150‑1 (MQ=255)
                                                                                   cATTACATTGCTGGATAAGAATGTTTTA‑GAAATCTCTATCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATa                                                                    <  1:1604567/150‑1 (MQ=255)
                                                                                     ttACATTGCTTGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAAc                                                                  <  1:861517/150‑1 (MQ=255)
                                                                                     ttACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAAc                                                                  <  1:763844/150‑1 (MQ=255)
                                                                                      tACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACa                                                                 >  2:532757/1‑150 (MQ=255)
                                                                                       aCATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAg                                                                >  1:194464/1‑150 (MQ=255)
                                                                                       aCATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAg                                                                <  2:1566377/150‑1 (MQ=255)
                                                                                       aCATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAg                                                                <  2:1199886/150‑1 (MQ=255)
                                                                                        cATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGc                                                               >  2:1635040/1‑150 (MQ=255)
                                                                                             cTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACtt                                                          >  1:453263/1‑150 (MQ=255)
                                                                                             cTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACtt                                                          <  1:558524/150‑1 (MQ=255)
                                                                                               ggATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGAACTGGCAAACTAACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTtat                                                        >  2:1084567/1‑150 (MQ=255)
                                                                                                  tAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATAcc                                                     <  1:1485384/150‑1 (MQ=255)
                                                                                                    aGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTg                                                   <  1:1528705/150‑1 (MQ=255)
                                                                                                      aaTGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTc                                                 <  2:1001928/150‑1 (MQ=255)
                                                                                                      aaTGTTTTA‑GCAAACTCTTTCTGTCATGAATACATGGCAGTGACCATACTAATTGTGACTGACATTTTTGGATGGATACACATTTCACTTGCAAACTAACACCGTTTTTTATCAACCTATACCATTAAATAACAGCAACTTATACCTGTc                                                 <  1:1122670/150‑1 (MQ=255)
                                                                                                         gTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGa                                              >  1:1272845/1‑150 (MQ=255)
                                                                                                          ttttA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGTGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATAAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGAt                                             >  1:90628/1‑150 (MQ=255)
                                                                                                          tattA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGAt                                             <  1:42573/148‑1 (MQ=255)
                                                                                                            ttA‑GCATTCTCTTTATGTCATGTATCCATGGCAGTGACCATACTAATGGTGACTGCCTTTGATGGAGGGAGACACAGTGCACTTGCAAACTCACACCGTTTTTAATCAACCTATACCATTATATAACAGCAACTTATACCTGTCTGATgg                                           <  2:1090528/150‑1 (MQ=255)
                                                                                                            ttA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATgg                                           <  1:259698/150‑1 (MQ=255)
                                                                                                                  aaTCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGgcgcgc                                      <  2:400408/150‑1 (MQ=255)
                                                                                                                   aTCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGct                                     >  1:1816057/1‑150 (MQ=255)
                                                                                                                    tctcTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGctc                                    <  1:377948/150‑1 (MQ=255)
                                                                                                                       ctTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGc                                 <  2:428368/150‑1 (MQ=255)
                                                                                                                        tttCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCg                                <  1:923481/150‑1 (MQ=255)
                                                                                                                            tGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACCCCGTTTTTAATCACCCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGc                            >  2:43904/1‑150 (MQ=255)
                                                                                                                             gTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCg                           >  1:1765129/1‑150 (MQ=255)
                                                                                                                              tCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCgg                          <  1:618796/150‑1 (MQ=255)
                                                                                                                               cATGAATCCATGGCAGTGACCATACTTATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATTCCTTTAAATAACAGCAACTTATACCTGTCTGTTGGCGCGCTCTGCGAAGCTGt                         <  1:1650307/150‑1 (MQ=255)
                                                                                                                               cATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGt                         >  1:1725939/1‑150 (MQ=255)
                                                                                                                               cATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGAAGGAGGGAGACACAGTGCCCTGGCAAACTCACACCGTTTTTCATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGt                         <  2:1103258/150‑1 (MQ=255)
                                                                                                                               cATGAATCCATGGCAGTGACCATACTAATGGTGACGGCCATAGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTGTTTAATCAACCAAGACCATTAAATAACAGCAACTTATCCCTGTCTGCTGGCCCGCTCTGAGAAGCGGt                         >  1:920351/1‑150 (MQ=255)
                                                                                                                                ctccaTCAATGGCAGTCACCATTAGTATGGTCACTGCCATGGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTa                        >  1:1748928/5‑150 (MQ=37)
                                                                                                                                 tGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTaa                       >  2:22638/1‑150 (MQ=255)
                                                                                                                                   aaTCCATGGCAGTGGCCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAAcg                     <  2:1220395/150‑1 (MQ=255)
                                                                                                                                    aTCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACCGCAACTTATACCTGTATGATGGCGCGCTCTGCGAAGCGGTAAcgc                    >  1:349739/1‑150 (MQ=255)
                                                                                                                                     tcaaTGGCAGTCACCATTAGTATGGTCACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAAcgcg                   >  2:1211218/4‑150 (MQ=39)
                                                                                                                                     tCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAAcgcg                   <  2:1798888/150‑1 (MQ=255)
                                                                                                                                       aaTGGCAGTCACCATTAGTATGGTCACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTg                 >  2:118435/2‑150 (MQ=255)
                                                                                                                                          ggCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAg              >  2:899184/1‑150 (MQ=255)
                                                                                                                                           gCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAgg             >  2:205761/1‑150 (MQ=255)
                                                                                                                                            cAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGt            >  2:569865/1‑150 (MQ=255)
                                                                                                                                              gTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTgc          >  2:1471087/1‑150 (MQ=255)
                                                                                                                                                 aCCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACCGGCAAACTCCCACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTgctgg       >  1:952315/1‑150 (MQ=255)
                                                                                                                                                   cATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTgctggct     >  1:512105/1‑150 (MQ=255)
                                                                                                                                                   cATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTgctggct     <  2:966492/150‑1 (MQ=255)
                                                                                                                                                    aTACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTgctggctg    >  2:280045/1‑150 (MQ=255)
                                                                                                                                                     tACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGctggctgg   <  2:1283654/150‑1 (MQ=255)
                                                                                                                                                     tACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGctggctgg   >  1:1009265/1‑150 (MQ=255)
                                                                                                                                                     tACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGctggctgg   <  1:1537028/150‑1 (MQ=255)
                                                                                                                                                      aCTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTggg  <  1:83128/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GGATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGG  >  NC_000913/4414091‑4414389

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