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.

GCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTTCTCATTACTATTTGCATACTGCCGTACCTTTGCTTTCTTTTCCTTG  >  NC_000913/3151807‑3152103
                                                                                                                                                      |                                                                                                                                                   
gctgctCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATAACGGACTTAACGACAATATCGGCGTGCTGATACATACACCAGACGGaaa                                                                                                                                                     >  1:129245/1‑150 (MQ=255)
   gctCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAATTAAGCCCATTCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAATGt                                                                                                                                                  >  2:1513188/1‑150 (MQ=255)
   gctCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGt                                                                                                                                                  >  2:223888/1‑150 (MQ=255)
         gcTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCt                                                                                                                                            <  2:615283/150‑1 (MQ=255)
          cTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGACCTTTCCGTCTg                                                                                                                                           <  2:1482615/150‑1 (MQ=255)
          cTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTg                                                                                                                                           <  1:78951/150‑1 (MQ=255)
          cTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTg                                                                                                                                           <  1:119446/150‑1 (MQ=255)
          cTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTg                                                                                                                                           >  2:875291/1‑150 (MQ=255)
          cTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTg                                                                                                                                           >  2:1536019/1‑150 (MQ=255)
           ttAAGGCGTCGCTTCTGAGTGAAGAACTCTACGCTCTTACTGATGAAGCTTGCCCAGGTGACTTCGGCGGCCAAAATAAGCCGATTCATCACGGTCTTTACGACAATTTCGGCGCGCTGATAGATACCCCAGACGGACCTTTCCGTCTgg                                                                                                                                          <  1:311055/150‑2 (MQ=255)
           ttAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGACCTTTCCGTCTgg                                                                                                                                          <  1:12712/150‑2 (MQ=255)
           ttAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGACCTTTCCGTCTgg                                                                                                                                          <  1:1476601/150‑2 (MQ=255)
           ttAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGc                                                                                                                                          <  2:344448/150‑1 (MQ=255)
            tAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGACCATTCCGTCTggg                                                                                                                                         <  2:753565/150‑3 (MQ=255)
              aGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGACCTTTCCGTCTggggt                                                                                                                                       <  1:696365/150‑5 (MQ=255)
              aGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCAtt                                                                                                                                       >  1:1436843/1‑150 (MQ=255)
              aGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCAtt                                                                                                                                       >  1:342119/1‑150 (MQ=255)
               ggCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTa                                                                                                                                      >  2:1544960/1‑150 (MQ=255)
                  gACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTATGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAat                                                                                                                                   <  1:736235/150‑1 (MQ=255)
                      cTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAattatt                                                                                                                               >  1:124745/1‑150 (MQ=255)
                        tCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTAc                                                                                                                             <  1:839466/150‑1 (MQ=255)
                          tGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTAccc                                                                                                                           <  2:549911/150‑1 (MQ=255)
                               gaagccctctACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAATTAAGCCCATTCATCACGCACTTACCGACAATATCGGCGTGCTGATACATACCCCAGACGGAACGGTCCGTCTGCATTAATTTATTACCCCCtgg                                                                                                                      >  2:1243201/7‑149 (MQ=255)
                               gaagaaCTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACtgt                                                                                                                      >  2:390812/1‑150 (MQ=255)
                                  gaaCTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTa                                                                                                                   <  2:1429123/150‑1 (MQ=255)
                                  gaaCTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTa                                                                                                                   <  1:195969/150‑1 (MQ=255)
                                    aCTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTAtc                                                                                                                 >  1:1063108/1‑150 (MQ=255)
                                      tctACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGACCTTTCCGTCTgg                                                                                                                                          >  2:302465/1‑122 (MQ=255)
                                      tctACGCTATTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCGGTCTGCATTAAATTATTACCCACTGTGTAtctc                                                                                                               <  1:352833/150‑1 (MQ=255)
                                       ctACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTcc                                                                                                              >  2:556113/1‑150 (MQ=255)
                                       ctACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTcc                                                                                                              >  2:1062951/1‑150 (MQ=255)
                                            ctctTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGAc                                                                                                         <  2:4373/150‑1 (MQ=255)
                                             tctTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACg                                                                                                        <  2:124745/150‑1 (MQ=255)
                                               ttACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCa                                                                                                      <  1:225999/150‑1 (MQ=255)
                                                tACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTTTTACCCACTGTGTATCTCCAGGACGCaa                                                                                                     >  1:1548381/1‑150 (MQ=255)
                                                 aCTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAg                                                                                                    <  2:880181/150‑1 (MQ=255)
                                                 aCTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAg                                                                                                    <  1:628280/150‑1 (MQ=255)
                                                  cTGAAGTAGATTGCCCAGGTGACTACGGAGGCCATAATAAGCCGAATCATCACGCACTTAACGACTATATCGGGGTGCTGATACATTCGCCAGACGGAAAGGTCCGTCTGCATTAATTTATTACCCACTGTGTATCTCCTGGACGCAAGt                                                                                                   <  1:898141/150‑1 (MQ=255)
                                                  cTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCATCGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAg                                                                                                    >  2:1450296/1‑150 (MQ=255)
                                                  cTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGt                                                                                                   <  2:43731/150‑1 (MQ=255)
                                                   tGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTc                                                                                                  <  1:1256922/150‑1 (MQ=255)
                                                   tGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTc                                                                                                  >  2:1491107/1‑150 (MQ=255)
                                                    gaagaaGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTca                                                                                                 >  2:609184/1‑150 (MQ=255)
                                                       gaagaTTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACaa                                                                                              >  2:812561/1‑150 (MQ=255)
                                                                 cAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCa                                                                                    >  1:1039263/1‑150 (MQ=255)
                                                                   ggTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCata                                                                                  >  1:1392996/1‑150 (MQ=255)
                                                                    gTGACTACGGAGGCCAAAATAAGCCCAATCATCCCGAACATAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCataa                                                                                 >  1:1089383/1‑150 (MQ=255)
                                                                         tACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCataataata                                                                            <  1:494105/150‑1 (MQ=255)
                                                                          aCGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAatat                                                                           >  1:756137/1‑150 (MQ=255)
                                                                           cGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATc                                                                          <  2:189849/150‑1 (MQ=255)
                                                                             gAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCaa                                                                        <  2:861314/150‑1 (MQ=255)
                                                                                 ccAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAacg                                                                    >  1:1419758/1‑150 (MQ=255)
                                                                                 ccAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAacg                                                                    >  1:1250432/1‑150 (MQ=255)
                                                                                 ccAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCCACAGCATTAAATTATTACCCACTCTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCAGAATAATATCACAacg                                                                    >  2:338569/1‑150 (MQ=255)
                                                                                  cAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAacga                                                                   >  2:876312/1‑150 (MQ=255)
                                                                                   aaaaTAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAacgac                                                                  <  1:364380/150‑1 (MQ=255)
                                                                                       tAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTTAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCg                                                              >  1:1288379/1‑150 (MQ=255)
                                                                                       tAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCg                                                              <  1:875291/150‑1 (MQ=255)
                                                                                         aGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGaa                                                            >  1:656645/1‑150 (MQ=255)
                                                                                         aGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGaa                                                            >  1:574635/1‑150 (MQ=255)
                                                                                          gCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAAt                                                           >  2:391195/1‑150 (MQ=255)
                                                                                           cccAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAAtt                                                          >  2:1296768/1‑150 (MQ=255)
                                                                                               atcatcACGCACTTAACGACAATATCGGCGTGCTGATACAGACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATa                                                      >  1:1602595/1‑150 (MQ=255)
                                                                                                tcttcaCGCACTTAACGACAATATCGGCGTGCTGATAGATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACGCACTGTGTATCTCCAGGACGCAAGTCTCAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAg                                                     <  2:1089383/147‑1 (MQ=255)
                                                                                                tcatcaCGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAg                                                     >  1:347342/1‑150 (MQ=255)
                                                                                                tcatcaCGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAg                                                     >  2:138039/1‑150 (MQ=255)
                                                                                                   tcaCGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCg                                                  <  2:107405/150‑1 (MQ=255)
                                                                                                     aCGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGtt                                                >  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)
                                                                                                                                                      |                                                                                                                                                   
GCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTCTGCATTAAATTATTACCCACTGTGTATCTCCAGGACGCAAGTCACAAAATCTGCGCATAATAATATCAAAACGACGTCGAATTGATAGTCGTTCTCATTACTATTTGCATACTGCCGTACCTTTGCTTTCTTTTCCTTG  >  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: