Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,561,206 A→T 9.0% intergenic (+308/‑162) yffL → / → yffM CPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,561,2060AT 9.0% 406.5 / 26.0 122intergenic (+308/‑162)yffL/yffMCPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (55/56);  new base T (5/6);  total (60/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-01

Only 100 of 127 total aligned reads displayed.

TTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTAAA  >  NC_000913/2561059‑2561353
                                                                                                                                                   |                                                                                                                                                   
ttAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAg                                                                                                                                                    >  3:948523/1‑149 (MQ=255)
 tAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAg                                                                                                                                                    >  2:207035/1‑148 (MQ=255)
 tAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAg                                                                                                                                                    >  1:826887/1‑148 (MQ=255)
   aGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCt                                                                                                                                                  >  1:111293/1‑148 (MQ=255)
   aGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCt                                                                                                                                                  <  2:687190/148‑1 (MQ=255)
   aGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCt                                                                                                                                                  >  2:517830/1‑148 (MQ=255)
    gAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTg                                                                                                                                                 >  2:119766/1‑148 (MQ=255)
       ttCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCTTGCTGGaa                                                                                                                                              <  1:141692/148‑1 (MQ=255)
       ttCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGaa                                                                                                                                              <  1:360748/148‑1 (MQ=255)
          tAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACg                                                                                                                                           <  1:181972/148‑1 (MQ=255)
          tAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATCTCCAGCAAGCTGGAAACGc                                                                                                                                          <  4:281130/149‑1 (MQ=255)
           aGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCTTGCTGGaaatga                                                                                                                                          <  2:562376/148‑4 (MQ=255)
               gACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTg                                                                                                                                                 >  3:110279/1‑137 (MQ=255)
                  gggAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtca                                                                                                                                   <  1:662225/148‑1 (MQ=255)
                    gAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcact                                                                                                                                 >  2:763261/1‑148 (MQ=255)
                     aGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGcc                                                                                                                                         <  2:906473/139‑1 (MQ=255)
                     aGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGcc                                                                                                                                         >  1:906473/1‑139 (MQ=255)
                       ttttttAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACg                                                                                                                             <  4:567242/149‑1 (MQ=255)
                        tttttAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACg                                                                                                                             >  1:539726/1‑148 (MQ=255)
                             aGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACaa                                                                                                                        <  1:517830/148‑1 (MQ=255)
                             aGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACaa                                                                                                                        <  2:826887/148‑1 (MQ=255)
                              gCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAAc                                                                                                                       <  1:185036/148‑1 (MQ=255)
                                cTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTg                                                                                                                    >  4:533341/1‑149 (MQ=255)
                                 tCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACg                                                                                                                                           <  1:317799/125‑1 (MQ=255)
                                 tCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACg                                                                                                                                           >  2:317799/1‑125 (MQ=255)
                                  cccGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGa                                                                                                                   >  2:530348/1‑148 (MQ=255)
                                    cggtcggtCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAAc                                                                                                                 >  1:97091/1‑148 (MQ=255)
                                      gtcggtcATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAAcagc                                                                                                              <  4:913600/149‑1 (MQ=255)
                                       tcggtcATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAAcagc                                                                                                              <  1:698580/148‑1 (MQ=255)
                                         ggtcATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGa                                                                                                           >  4:270660/1‑149 (MQ=255)
                                              tAGAGAGGCGCCTAGAGTTAAACAGAAGAATTAGTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACcaa                                                                                                       >  2:567282/1‑148 (MQ=255)
                                               agagagTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACcaac                                                                                                      >  2:961303/1‑148 (MQ=255)
                                                gagagTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAAca                                                                                                                >  2:842494/1‑137 (MQ=255)
                                                gagagTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACcaacaa                                                                                                    >  4:471474/1‑149 (MQ=255)
                                                     tCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGaaa                                                                                                                                             >  4:260071/1‑103 (MQ=255)
                                                     tCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGaaa                                                                                                                                             <  3:260071/103‑1 (MQ=255)
                                                      cGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACcat                                                                                               <  1:530348/148‑1 (MQ=255)
                                                       gCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACcatca                                                                                             <  4:102363/149‑1 (MQ=255)
                                                       gCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCCGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACcatca                                                                                             <  4:597420/149‑1 (MQ=255)
                                                        ccTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACcatca                                                                                             <  1:703856/148‑1 (MQ=255)
                                                         cTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACcaac                                                                                                      <  4:90276/138‑1 (MQ=255)
                                                         cTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACcaac                                                                                                      >  3:90276/1‑138 (MQ=255)
                                                         cTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACcatcat                                                                                            <  1:893824/148‑1 (MQ=255)
                                                          tAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATg                                                                                           >  1:749842/1‑148 (MQ=255)
                                                                 aaaCAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCg                                                                                    <  1:119766/148‑1 (MQ=255)
                                                                  aaCAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGtt                                                                                  >  3:945187/1‑149 (MQ=255)
                                                                  aaCAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGt                                                                                   >  2:289911/1‑148 (MQ=255)
                                                                   aCAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCTTGCTGGaaat                                                                                                                                            >  3:627740/1‑89 (MQ=255)
                                                                       aagaaTTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATgg                                                                             >  3:361563/1‑149 (MQ=255)
                                                                       aagaaTTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATgg                                                                             <  4:812038/149‑1 (MQ=255)
                                                                          aaTTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATgggg                                                                           <  1:207035/148‑1 (MQ=255)
                                                                           attattTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGa                                                                          >  1:48311/1‑148 (MQ=255)
                                                                                ttCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACcat                                                                                               >  4:229105/1‑122 (MQ=255)
                                                                                ttCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACcat                                                                                               <  3:229105/122‑1 (MQ=255)
                                                                                 tCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCgg                                                                    <  1:846509/148‑1 (MQ=255)
                                                                                    ccgccTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGtata                                                                >  3:565254/1‑149 (MQ=255)
                                                                                      gccTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTatat                                                               <  2:111293/148‑1 (MQ=255)
                                                                                       ccTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATAtc                                                              >  1:405019/1‑148 (MQ=255)
                                                                                         tGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCa                                                           >  4:391687/1‑149 (MQ=255)
                                                                                           aTAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCaaa                                                         <  4:1947/149‑1 (MQ=255)
                                                                                             aGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAgg                                                       <  3:499189/149‑1 (MQ=255)
                                                                                             aGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAg                                                        >  2:147606/1‑148 (MQ=255)
                                                                                               gTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGt                                                     >  3:809296/1‑149 (MQ=255)
                                                                                                 cGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTg                                                    <  2:162734/148‑1 (MQ=255)
                                                                                                 cGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTg                                                    >  1:322062/1‑148 (MQ=255)
                                                                                                      ggTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCt                                              <  3:428719/149‑1 (MQ=255)
                                                                                                         tGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAgcgc                                                                      >  2:468639/1‑122 (MQ=255)
                                                                                                         tGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAgcgc                                                                      <  1:468639/122‑1 (MQ=255)
                                                                                                            gAGCATTAACCATGCTGATTAACGCAATCCCTTCCGGCAAGCTGGAAACGCCCATTCACTCACGCACACCCTGAACAGCAGACCACCAAGCCATCATGGCATTCGTTCATGGCGAGCGCAGTCGATCTCAAAGGGTGATTTCTTATtgg                                        >  4:592162/1‑149 (MQ=255)
                                                                                                             aGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtgg                                        >  1:104772/1‑148 (MQ=255)
                                                                                                             aGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtgg                                        <  1:669274/148‑1 (MQ=255)
                                                                                                             aGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtgg                                        >  2:355987/1‑148 (MQ=255)
                                                                                                               cATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtggtgg                                     <  3:222299/149‑1 (MQ=255)
                                                                                                                 ttAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGt                                                     <  3:796882/131‑1 (MQ=255)
                                                                                                                 ttAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGt                                                     >  4:796882/1‑131 (MQ=255)
                                                                                                                 ttAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTggg                                    >  1:28539/1‑148 (MQ=255)
                                                                                                                 ttAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGc                                   >  4:267715/1‑149 (MQ=255)
                                                                                                                    aCCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTa                                            >  4:676111/1‑137 (MQ=255)
                                                                                                                    aCCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTa                                            <  3:676111/137‑1 (MQ=255)
                                                                                                                     ccATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCtt                                <  2:549387/148‑1 (MQ=255)
                                                                                                                         gCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGaa                           <  4:131746/149‑1 (MQ=255)
                                                                                                                           tGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAGCAAAACATCATGGCATTCGTTCATTGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAgg                         <  4:687200/149‑1 (MQ=255)
                                                                                                                             aTTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAAcagcag                                                                                                            >  2:407963/1‑64 (MQ=255)
                                                                                                                             aTTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAAcagcag                                                                                                            <  1:407963/64‑1 (MQ=255)
                                                                                                                             aTTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGc                        <  1:961303/148‑1 (MQ=255)
                                                                                                                             aTTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCt                       >  4:634250/1‑149 (MQ=255)
                                                                                                                               tAGCGCAATCATTTCCAGCAAGCTGGAACCGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTc                     <  4:808284/149‑1 (MQ=255)
                                                                                                                                 gcgcAATCATTTCCAGCAAGCTGGAAACGCCCCTTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGa                   <  4:777269/149‑1 (MQ=255)
                                                                                                                                       tCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATg                                                                                           >  4:847918/1‑71 (MQ=255)
                                                                                                                                       tCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATg                                                                                           <  3:847918/71‑1 (MQ=255)
                                                                                                                                        cATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtg                                         >  2:201073/1‑120 (MQ=255)
                                                                                                                                        cATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtg                                         <  1:201073/120‑1 (MQ=255)
                                                                                                                                          tttCCAGCTTGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTAtt                                          >  4:414576/1‑117 (MQ=255)
                                                                                                                                          tttCCAGCTTGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTAtt                                          <  3:414576/117‑1 (MQ=255)
                                                                                                                                          tttCCAGCTTGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGcc          >  3:850832/1‑149 (MQ=255)
                                                                                                                                           ttCCAGCTTGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGAAATTGcc          >  1:858279/1‑148 (MQ=255)
                                                                                                                                             ccAGCTTGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTc                               >  4:714635/1‑125 (MQ=255)
                                                                                                                                             ccAGCTTGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTc                               <  3:714635/125‑1 (MQ=255)
                                                                                                                                                   aGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCGTGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTaaa  >  2:627198/1‑148 (MQ=255)
                                                                                                                                                   aGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACAAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTaaa  <  1:837925/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTAAA  >  NC_000913/2561059‑2561353

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: