Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,151,956 A→C 7.3% L10R (CTT→CGT) ‡ exbB ← membrane spanning protein in TonB‑ExbB‑ExbD complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,151,9560AC 7.3% 348.7 / 11.8 116L10R (CTT→CGT) ‡exbBmembrane spanning protein in TonB‑ExbB‑ExbD complex
Reads supporting (aligned to +/- strand):  ref base A (62/45);  new base C (2/7);  total (64/52)
Fisher's exact test for biased strand distribution p-value = 7.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 119 total aligned reads displayed.

GCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTTCTCATTACTATTTGCATACTGCCGTACCTTTGCTTTCTTTTCCTTG  >  NC_000913/3151807‑3152103
                                                                                                                                                      |                                                                                                                                                   
gctgctCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATAA‑CGGACTTAACGACAATATCGGCGTGCTGATACATACACCAGACGGaaa                                                                                                                                                     >  1:129245/1‑150 (MQ=255)
   gctCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAATTAAGCCCATTCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAATGt                                                                                                                                                  >  2:1513188/1‑150 (MQ=255)
   gctCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGt                                                                                                                                                  >  2:223888/1‑150 (MQ=255)
         gcTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCt                                                                                                                                            <  2:615283/150‑1 (MQ=255)
          cTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGACCTTTCCGTCTg                                                                                                                                           <  2:1482615/150‑1 (MQ=255)
          cTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTg                                                                                                                                           <  1:78951/150‑1 (MQ=255)
          cTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTg                                                                                                                                           <  1:119446/150‑1 (MQ=255)
          cTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTg                                                                                                                                           >  2:875291/1‑150 (MQ=255)
          cTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTg                                                                                                                                           >  2:1536019/1‑150 (MQ=255)
           ttAAGGCGTCGCTTCTGAGTGAAGAACTCTACGCTCTTACTGATGAAGCTTGCCCAGGTGACTTCGGCGGCCAAAATAAGCCGATTCATCA‑CGGTCTTTACGACAATTTCGGCGCGCTGATAGATACCCCAGACGGACCTTTCCGTCTgg                                                                                                                                          <  1:311055/150‑2 (MQ=255)
           ttAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGACCTTTCCGTCTgg                                                                                                                                          <  1:12712/150‑2 (MQ=255)
           ttAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGACCTTTCCGTCTgg                                                                                                                                          <  1:1476601/150‑2 (MQ=255)
           ttAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGc                                                                                                                                          <  2:344448/150‑1 (MQ=255)
            tAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGACCATTCCGTCTggg                                                                                                                                         <  2:753565/150‑3 (MQ=255)
              aGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGACCTTTCCGTCTggggt                                                                                                                                       <  1:696365/150‑5 (MQ=255)
              aGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCAtt                                                                                                                                       >  1:1436843/1‑150 (MQ=255)
              aGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCAtt                                                                                                                                       >  1:342119/1‑150 (MQ=255)
               ggCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTa                                                                                                                                      >  2:1544960/1‑150 (MQ=255)
                  gACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTATGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAat                                                                                                                                   <  1:736235/150‑1 (MQ=255)
                      cTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAattatt                                                                                                                               >  1:124745/1‑150 (MQ=255)
                        tCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTAc                                                                                                                             <  1:839466/150‑1 (MQ=255)
                          tGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTAccc                                                                                                                           <  2:549911/150‑1 (MQ=255)
                               gaagccctctACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAATTAAGCCCATTCATCA‑CGCACTTACCGACAATATCGGCGTGCTGATACATACCCCAGACGGAACGGTCCGTCTGCATTAATTTATTACCCCCtgg                                                                                                                      >  2:1243201/7‑149 (MQ=255)
                               gaagaaCTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACtgt                                                                                                                      >  2:390812/1‑150 (MQ=255)
                                  gaaCTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTa                                                                                                                   <  2:1429123/150‑1 (MQ=255)
                                  gaaCTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTa                                                                                                                   <  1:195969/150‑1 (MQ=255)
                                    aCTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTAtc                                                                                                                 >  1:1063108/1‑150 (MQ=255)
                                      tctACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGACCTTTCCGTCTgg                                                                                                                                          >  2:302465/1‑122 (MQ=255)
                                      tctACGCTATTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCGGTCTGCATTAAATTATTACCCACTGTGTAtctc                                                                                                               <  1:352833/150‑1 (MQ=255)
                                       ctACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTcc                                                                                                              >  2:556113/1‑150 (MQ=255)
                                       ctACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTcc                                                                                                              >  2:1062951/1‑150 (MQ=255)
                                            ctctTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGAc                                                                                                         <  2:4373/150‑1 (MQ=255)
                                             tctTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACg                                                                                                        <  2:124745/150‑1 (MQ=255)
                                               ttACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCa                                                                                                      <  1:225999/150‑1 (MQ=255)
                                                tACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTTTTACCCACTGTGTATCTCCAGGACGCaa                                                                                                     >  1:1548381/1‑150 (MQ=255)
                                                 aCTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAg                                                                                                    <  2:880181/150‑1 (MQ=255)
                                                 aCTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAg                                                                                                    <  1:628280/150‑1 (MQ=255)
                                                  cTGAAGTAGATTGCCCAGGTGACTACGGAGGCCATAATAAGCCGAATCATCA‑CGCACTTAACGACTATATCGGGGTGCTGATACATTCGCCAGACGGAAAGGTCCGTCTGCATTAATTTATTACCCACTGTGTATCTCCTGGACGCAAGt                                                                                                   <  1:898141/150‑1 (MQ=255)
                                                  cTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCATCGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAg                                                                                                    >  2:1450296/1‑150 (MQ=255)
                                                  cTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGt                                                                                                   <  2:43731/150‑1 (MQ=255)
                                                   tGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTc                                                                                                  <  1:1256922/150‑1 (MQ=255)
                                                   tGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTc                                                                                                  >  2:1491107/1‑150 (MQ=255)
                                                    gaagaaGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTca                                                                                                 >  2:609184/1‑150 (MQ=255)
                                                       gaagaTTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACaa                                                                                              >  2:812561/1‑150 (MQ=255)
                                                                 cAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCa                                                                                    >  1:1039263/1‑150 (MQ=255)
                                                                   ggTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCata                                                                                  >  1:1392996/1‑150 (MQ=255)
                                                                    gTGACTACGGAGGCCAAAATAAGCCCAATCATCC‑CGAACATAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCataa                                                                                 >  1:1089383/1‑150 (MQ=255)
                                                                         tACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCataataata                                                                            <  1:494105/150‑1 (MQ=255)
                                                                          aCGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAatat                                                                           >  1:756137/1‑150 (MQ=255)
                                                                           cGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATc                                                                          <  2:189849/150‑1 (MQ=255)
                                                                             gAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCaa                                                                        <  2:861314/150‑1 (MQ=255)
                                                                                 ccAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAacg                                                                    >  1:1419758/1‑150 (MQ=255)
                                                                                 ccAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAacg                                                                    >  1:1250432/1‑150 (MQ=255)
                                                                                 ccAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCCACAGCATTAAATTATTACCCACTCTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCAGAATAATATCACAacg                                                                    >  2:338569/1‑150 (MQ=255)
                                                                                  cAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAacga                                                                   >  2:876312/1‑150 (MQ=255)
                                                                                   aaaaTAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAacgac                                                                  <  1:364380/150‑1 (MQ=255)
                                                                                       tAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTTAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCg                                                              >  1:1288379/1‑150 (MQ=255)
                                                                                       tAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCg                                                              <  1:875291/150‑1 (MQ=255)
                                                                                         aGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGaa                                                            >  1:656645/1‑150 (MQ=255)
                                                                                         aGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGaa                                                            >  1:574635/1‑150 (MQ=255)
                                                                                          gCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAAt                                                           >  2:391195/1‑150 (MQ=255)
                                                                                           cccAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAAtt                                                          >  2:1296768/1‑150 (MQ=255)
                                                                                               atcatcA‑CGCACTTAACGACAATATCGGCGTGCTGATACAGACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATa                                                      >  1:1602595/1‑150 (MQ=255)
                                                                                                tcttca‑CGCACTTAACGACAATATCGGCGTGCTGATAGATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACGCACTGTGTATCTCCAGGACGCAAGTCTCAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAg                                                     <  2:1089383/147‑1 (MQ=255)
                                                                                                tcatca‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAg                                                     >  1:347342/1‑150 (MQ=255)
                                                                                                tcatca‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAg                                                     >  2:138039/1‑150 (MQ=255)
                                                                                                   tca‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCg                                                  <  2:107405/150‑1 (MQ=255)
                                                                                                     a‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGtt                                                >  2:388874/1‑150 (MQ=255)
                                                                                                        gCACTTAACGACAATATCGGCGTGCTGATATATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTtct                                              >  2:1256504/1‑150 (MQ=255)
                                                                                                         cACTTAACGACAATATCGGCGTGCTGCTACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTtctc                                             >  2:503324/1‑150 (MQ=255)
                                                                                                         cACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTtctc                                             >  1:525157/1‑150 (MQ=255)
                                                                                                             tAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTTCTCATTa                                         >  1:1059654/1‑150 (MQ=255)
                                                                                                             tAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTTCTCATTa                                         <  1:876312/150‑1 (MQ=255)
                                                                                                             tAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTTCTCATTa                                         >  1:724003/1‑150 (MQ=255)
                                                                                                                cGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTTCTCATTACTa                                      >  2:791886/1‑150 (MQ=255)
                                                                                                                    aaTATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTTCTCATTACTATTTg                                  >  2:1013108/1‑150 (MQ=255)
                                                                                                                    aaTATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTTCTCATTACTATTTg                                  <  2:787717/150‑1 (MQ=255)
                                                                                                                      tatCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTTCTCATTACTATTTGCa                                <  1:1536019/150‑1 (MQ=255)
                                                                                                                       atCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTTCTCATTACTATTTGCAt                               <  1:1662599/150‑1 (MQ=255)
                                                                                                                        tCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTTCTCATTACTATTTGCATa                              <  2:1436843/150‑1 (MQ=255)
                                                                                                                         cGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATATTAATATCAAAACGACGTCGAATTGATAGTCGTTCTCATTACTATTTGCATAc                             <  1:1544960/150‑1 (MQ=255)
                                                                                                                               gCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTTCTCATTACTATTTGCATACTGCCGt                       >  2:1415390/1‑150 (MQ=255)
                                                                                                                                  gATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTTCTCATTACTATTTGCATACTGCCGTAcc                    <  2:1039263/150‑1 (MQ=255)
                                                                                                                                     acatacCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTTCTCATTACTATTTGCATACTGCCGTACCttt                 <  1:1450296/150‑1 (MQ=255)
                                                                                                                                     acatacCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCAAAATAATATCAAAACGACGTCGAATTGATAGTCGTTCTCATTACTATTTGCATACTGCCGTACCttt                 >  2:1227120/1‑150 (MQ=255)
                                                                                                                                      catacCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGCCGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTTCTCATTACTATTTGTATACTGCCGTACCTTTg                >  1:902307/1‑150 (MQ=255)
                                                                                                                                       atacCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTTCTCATTACTATTTGCATACTGCCGTACCTTTGc               >  1:383004/1‑150 (MQ=255)
                                                                                                                                       atacCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTTCTCATTACTATTTGCATACTGCCGTACCTTTGc               >  2:753786/1‑150 (MQ=255)
                                                                                                                                        tacCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTTCTCATTACTATTTGCATACTGCCGTACCTTTGct              <  1:556113/150‑1 (MQ=255)
                                                                                                                                         acCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACTACGTCGAATTGATAGTCGTTCTCATTACTATTTGCATACTGCCGTACCTTTGctt             <  2:96532/150‑1 (MQ=255)
                                                                                                                                           cccAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTTCTCATTACTATTTGCATACTGCCGTACCTTTGctttc           <  1:390812/150‑1 (MQ=255)
                                                                                                                                            ccAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTTCTCATTACTATTTGCATACTGCCGTACCTTTGctttct          <  1:223888/150‑1 (MQ=255)
                                                                                                                                            ccAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTTCTCATTACTATTTGCATACTGCCGTACCTTTGctttct          <  1:1653711/150‑1 (MQ=255)
                                                                                                                                            ccAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTTCTCATTACTATATGCATACTGCCGTACCTTTGctttct          <  1:1513188/150‑1 (MQ=255)
                                                                                                                                             cAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTTCTCATTACTATTTGCATACTGCCGTACCTTTGctttctt         >  1:1379031/1‑150 (MQ=255)
                                                                                                                                              ataCGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAACATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTTCTCATTACTATTTGCATACTGCCGTACCTTTGctttcttt        <  2:129245/148‑1 (MQ=255)
                                                                                                                                                 cGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTTCTCATTACTATTTGCATACTGCCGTACCTTTGCTTTCTTTTcc     >  1:331306/1‑150 (MQ=255)
                                                                                                                                                  gcaaaGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATTTCAAAACGACGTCGAATTGATAGTCGTTCTCATTACTATTTGCATACTGCCGTACCTTTGCTTTCTTTTCCt    <  2:240987/148‑1 (MQ=255)
                                                                                                                                                   gAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTTCTCATTACTATTTGCATACTGCCGTACCTTTGCTTTCTTTTCCtt   >  1:573862/1‑150 (MQ=255)
                                                                                                                                                   gAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTTCTCATTACTATTTGCATACTGCCGTACCTTTGCTTTCTTTTCCtt   >  1:580682/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCA‑CGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTTCTCATTACTATTTGCATACTGCCGTACCTTTGCTTTCTTTTCCTTG  >  NC_000913/3151807‑3152103

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