Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,550,578 C→G 6.1% intergenic (‑33/+63) yfeX ← / ← yfeY porphyrinogen oxidase, cytoplasmic/RpoE‑regulated lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,550,5780CG 6.1% 510.3 / 16.2 158intergenic (‑33/+63)yfeX/yfeYporphyrinogen oxidase, cytoplasmic/RpoE‑regulated lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (76/72);  new base G (5/5);  total (81/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

Only 100 of 160 total aligned reads displayed.

GTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTCC  >  NC_000913/2550429‑2550725
                                                                                                                                                     |                                                                                                                                                   
 tCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATcc                                                                                                                                                    >  1:1068749/1‑150 (MQ=255)
    gctaaTGTTTTACTGGCCGCACGCAGGGCGGCAACTTCCCCTTTCACGTTGGCTTCGATCCTAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTg                                                                                                                                                 <  2:698699/147‑1 (MQ=255)
     cAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTgg                                                                                                                                                <  2:1280898/150‑1 (MQ=255)
     cAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTgg                                                                                                                                                <  1:1668609/150‑1 (MQ=255)
      aaaTGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGc                                                                                                                                               <  1:924253/150‑1 (MQ=255)
      aaaTGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGc                                                                                                                                               <  2:1083354/150‑1 (MQ=255)
      aaaTGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGc                                                                                                                                               <  2:243710/150‑1 (MQ=255)
         tGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCAtt                                                                                                                                            <  2:1760530/150‑1 (MQ=255)
          gTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCCCGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGACACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCAttg                                                                                                                                           >  2:274327/1‑150 (MQ=255)
           ttttACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCAttgt                                                                                                                                          <  1:536590/150‑1 (MQ=255)
           ttttACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCAttgt                                                                                                                                          >  2:1161193/1‑150 (MQ=255)
            tttACTGGCCGCTCGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAACAATGCCACGGATTGCGTGGCattatt                                                                                                                                         <  1:1940553/150‑1 (MQ=255)
            tttACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCTCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGGATTGCGTGGCattatt                                                                                                                                         <  2:642949/150‑1 (MQ=255)
             ttACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCAgtgttg                                                                                                                                        >  2:1529286/1‑150 (MQ=255)
              tACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGa                                                                                                                                       <  2:83648/150‑1 (MQ=255)
               aCTGGCCGCACGCAGGGCGTCAACTTCACCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGAAATGTTATGGCAAAATGCCACTCTGAACCTGAGACCTTGTTCCTCCTGAAAATAATAATGCCACTCAATCCGTGGCATTGTTGAg                                                                                                                                      <  1:1973122/150‑1 (MQ=255)
                cTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAgc                                                                                                                                     <  2:16157/150‑1 (MQ=255)
                 tGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAgcg                                                                                                                                    >  1:1577260/1‑150 (MQ=255)
                  ggCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAACTGCCTCTCTGACCCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACCGATTGCGTGGCattattattttc                                                                                                                                   <  1:292351/150‑7 (MQ=255)
                       cacgcaGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTAttt                                                                                                                              <  2:486943/150‑1 (MQ=255)
                         cgcaGGGCGTCAACTTCCCCTTTCACGTTGGCATCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTa                                                                                                                            <  2:2014943/150‑1 (MQ=255)
                         cgcaGGGCGTCAACTTCCCCTTTCACGGTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTa                                                                                                                            <  2:1998241/150‑1 (MQ=255)
                               gCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGtttt                                                                                                                      >  1:316898/1‑150 (MQ=255)
                               gCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGtt                                                                                                                      >  1:1949549/1‑150 (MQ=255)
                                      cTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTcgc                                                                                                               >  2:502993/1‑150 (MQ=255)
                                      cTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTcgc                                                                                                               >  1:1646397/1‑150 (MQ=255)
                                       ttCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTcgcg                                                                                                              <  2:483110/150‑1 (MQ=255)
                                        tCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCgcgg                                                                                                             >  1:238829/1‑150 (MQ=255)
                                         ccccTTTCACGTTGGCTTCGATCAAAATAGCCGAGCGGCAATGTTCTGGCAAAATGCCAATCTGAACCTGAGACATTGTTCCTCCTGAAAATAATATTGCCACGCAATCCGTGGCATTGTTGAGAGCTATTTTACGTGTTTTTTCtcggc                                                                                                            >  2:69653/1‑150 (MQ=255)
                                                cACGGTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGtt                                                                                                     <  2:1676423/150‑1 (MQ=255)
                                                  cGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTg                                                                                                   >  1:712347/1‑150 (MQ=255)
                                                   gTTGGCTTCGATCCACATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTCCCTCCTCTAAATAATAATTCCACGCAATCCGGGGCTTTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGt                                                                                                  <  1:1790481/150‑1 (MQ=255)
                                                   gTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATTCCACTCTGAACCTGAGTCATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGt                                                                                                  <  1:1530532/150‑1 (MQ=255)
                                                     tGGCTGCGATCCAAATCGCCGCGCGGCAATGTTCTGGAAAAAGGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTc                                                                                                <  2:188441/150‑1 (MQ=255)
                                                      ggCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATATTAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCa                                                                                               >  1:1500943/1‑150 (MQ=255)
                                                       gCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAg                                                                                              <  2:1286721/150‑1 (MQ=255)
                                                         ttCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGAc                                                                                            <  2:1917654/150‑1 (MQ=255)
                                                          tCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACg                                                                                           >  2:808692/1‑150 (MQ=255)
                                                            gATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGACCCTGAGACATTGTTCATCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCa                                                                                         >  2:32685/1‑150 (MQ=255)
                                                            gATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCa                                                                                         >  2:275888/1‑150 (MQ=255)
                                                            gATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCa                                                                                         >  1:1278146/1‑150 (MQ=255)
                                                            gATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCa                                                                                         >  1:211920/1‑150 (MQ=255)
                                                              tCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCaaa                                                                                       >  2:898659/1‑150 (MQ=255)
                                                               ccAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAAt                                                                                      >  1:1919299/1‑150 (MQ=255)
                                                               ccAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAAt                                                                                      >  1:454601/1‑150 (MQ=255)
                                                                cAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAAtt                                                                                     >  1:1128885/1‑150 (MQ=255)
                                                                cAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAAtt                                                                                     >  2:811448/1‑150 (MQ=255)
                                                                  aaTCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTaa                                                                                   >  2:741846/1‑150 (MQ=255)
                                                                   aTCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGCACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAAc                                                                                  >  1:592422/1‑150 (MQ=255)
                                                                   aTCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGCAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAAc                                                                                  <  1:701645/150‑1 (MQ=255)
                                                                    tCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAAcg                                                                                 <  2:877640/150‑1 (MQ=255)
                                                                      gccgcGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAAcgcc                                                                               <  2:1133139/150‑1 (MQ=255)
                                                                       ccgcGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAAcgccg                                                                              <  1:1909409/150‑1 (MQ=255)
                                                                          cgcgGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCCTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCa                                                                           <  1:422141/150‑1 (MQ=255)
                                                                           gcgGCAATGTTCTGGCAAAATGCCACTCTGACCCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCTTTTTTGAGCGCTATTTTACTTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAg                                                                          >  2:1324011/1‑150 (MQ=255)
                                                                              gCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATa                                                                       >  1:1421075/1‑150 (MQ=255)
                                                                               cAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATat                                                                      >  1:62763/1‑149 (MQ=255)
                                                                                  tGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAAttt                                                                   >  2:405199/1‑150 (MQ=255)
                                                                                      gtGGCAAAATGCCACTCTGAACCTTAGACATTGTTACTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACt                                                               <  2:1936089/149‑1 (MQ=255)
                                                                                      cTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACt                                                               <  1:1149650/150‑1 (MQ=255)
                                                                                       tGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTg                                                              <  1:1633849/150‑1 (MQ=255)
                                                                                        ggCAAAATGCCACTCTTAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGCGCGCTATTTTACGTGTTTTTTCTCGGCGGCGGTTTGTTCAGACGCAAATTAACGCAGCAAGATAATTTTACTGa                                                             >  1:147244/1‑150 (MQ=255)
                                                                                          cAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGTGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACt                                                           >  2:1104622/1‑150 (MQ=255)
                                                                                          cAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACt                                                           >  2:1403707/1‑150 (MQ=255)
                                                                                          cAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACt                                                           >  1:1724538/1‑150 (MQ=255)
                                                                                          cAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACt                                                           >  1:1759207/1‑150 (MQ=255)
                                                                                              aTGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTcc                                                       <  2:1915176/150‑1 (MQ=255)
                                                                                                   aCTCTGATCCTGAGACATTGTTCCTCCTTAAAATTATAATGCCACGCTATCCGTGGCATTGTTGAGCGCTTTTTTACTTGTTTTTTCGCGGCGGCGGGTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTAATGACTTTCCAGttt                                                  <  2:1565770/150‑1 (MQ=255)
                                                                                                     tctGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGttttt                                                <  1:1621891/150‑1 (MQ=255)
                                                                                                     tctGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGttttt                                                <  2:1444071/150‑1 (MQ=255)
                                                                                                      ctGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTc                                               <  1:809800/150‑1 (MQ=255)
                                                                                                        gAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAg                                             <  1:1056337/150‑1 (MQ=255)
                                                                                                        gAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACATGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAg                                             <  2:1494904/150‑1 (MQ=255)
                                                                                                         aaCCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAgg                                            <  2:1982528/150‑1 (MQ=255)
                                                                                                         aaCCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAgg                                            <  1:719303/150‑1 (MQ=255)
                                                                                                             tGAGACATTGTTCCTCCTGAAAATTATAATGCCACGCAATCCGTGGCATTGTTGAGCGGTTTTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTAAGTCGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTAt                                        >  2:1350295/1‑150 (MQ=255)
                                                                                                             tGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTAt                                        <  2:1184774/150‑1 (MQ=255)
                                                                                                              gataCATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATc                                       >  1:1377293/4‑150 (MQ=255)
                                                                                                                 aCATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTc                                    <  1:1889244/150‑1 (MQ=255)
                                                                                                                     tGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGaa                                <  2:1562886/150‑1 (MQ=255)
                                                                                                                      gTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAg                               >  1:915104/1‑150 (MQ=255)
                                                                                                                       ttCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAgg                              <  2:125383/150‑1 (MQ=255)
                                                                                                                         cctcctGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCa                            <  2:1199519/150‑1 (MQ=255)
                                                                                                                         cctccaGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCa                            <  2:1292341/150‑1 (MQ=255)
                                                                                                                           tcctGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTCACGCCGCCAGATTATTTTACTGACTTTCCAGTTTTTCATGTTATCTTCCGAAGGCAtt                          >  1:338663/1‑150 (MQ=255)
                                                                                                                           tcctGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCAtt                          >  1:152415/1‑150 (MQ=255)
                                                                                                                            cctGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTTCGAAGGCATTa                         >  2:1982947/1‑150 (MQ=255)
                                                                                                                            cctGAAAATAATAATGCCACGCAATCCGTAGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTa                         >  2:680203/1‑150 (MQ=255)
                                                                                                                             ctGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTaa                        >  2:1776376/1‑150 (MQ=255)
                                                                                                                              tGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAAc                       <  1:132131/150‑1 (MQ=255)
                                                                                                                              tGAAAATAATAATGCCACGCAATCCGTGGAATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAAc                       <  2:1254934/150‑1 (MQ=255)
                                                                                                                               gAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAAcc                      >  2:1545539/1‑150 (MQ=255)
                                                                                                                               gAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAAcc                      >  2:1690486/1‑150 (MQ=255)
                                                                                                                                 aaaTAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGAGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGTTATCGTCAGAAGGCATTAACCCt                    >  2:377032/1‑150 (MQ=255)
                                                                                                                                    taataatGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCa                 <  1:1322406/150‑1 (MQ=255)
                                                                                                                                        aatGCCACGGATTGCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGAc             >  2:423108/1‑150 (MQ=255)
                                                                                                                                        aatGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGAc             <  2:316898/150‑1 (MQ=255)
                                                                                                                                        aatGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGAc             >  1:808367/1‑150 (MQ=255)
                                                                                                                                             cACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGc        >  2:1872179/1‑150 (MQ=255)
                                                                                                                                                   aTCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTcc  >  2:1783341/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTCC  >  NC_000913/2550429‑2550725

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: