Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,775,503 C→T 5.7% intergenic (+482/‑42) yfjW → / → ypjI CP4‑57 prophage; putative inner membrane protein/pseudogene, CP4‑57 putative prophage remnant;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,775,5030CT 5.7% 445.0 / 10.1 140intergenic (+482/‑42)yfjW/ypjICP4‑57 prophage; putative inner membrane protein/pseudogene, CP4‑57 putative prophage remnant;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (62/69);  new base T (2/6);  total (65/75)
Fisher's exact test for biased strand distribution p-value = 2.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.12e-01

Only 100 of 142 total aligned reads displayed.

GTAAATGCCCCTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGCCGGTATCACC  >  NC_000913/2775356‑2775652
                                                                                                                                                   |                                                                                                                                                     
gTAAATGCCCCTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAg                                                                                                                                                     <  1:1609777/150‑1 (MQ=255)
gTAAATGCCCCTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAg                                                                                                                                                     >  2:1345069/1‑150 (MQ=255)
gTAAATGCCCCTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTACGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCTATTAACGGCAGAATCTCAg                                                                                                                                                     >  2:949438/1‑150 (MQ=255)
   aaTGCCCCTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGAt                                                                                                                                                  >  1:1522062/1‑150 (MQ=255)
         ccTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCCTGAGATTCTGcc                                                                                                                                            <  2:69730/150‑1 (MQ=255)
          cTACATAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCTCATTTGGCACAGAACGGACTTTCAGCAATATTAAAGCATTCAGATATATATCTTTTCTACAATCATACCCGATTAACGGCAGAATCTAAGGATTCTGCCg                                                                                                                                           >  1:1847664/1‑150 (MQ=255)
          cTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAGCTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCg                                                                                                                                           >  1:172929/1‑150 (MQ=255)
           tACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGc                                                                                                                                          >  2:1083372/1‑150 (MQ=255)
            aCAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTTGCACAGAACGGCATATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCt                                                                                                                                         <  1:1508116/150‑1 (MQ=255)
            aCAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCCTGAGATTCTGCCgtt                                                                                                                                         <  2:431497/150‑3 (MQ=255)
              aaaaTGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCttt                                                                                                                                       >  2:1090534/1‑150 (MQ=255)
                 aTGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTTAGGATTCTGCCGCtttttt                                                                                                                                    <  1:353336/150‑1 (MQ=255)
                 aTGAACCATTATTTTTAGAGGGAGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACAGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCtttttt                                                                                                                                    <  2:152363/150‑1 (MQ=255)
                  tGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTtta                                                                                                                                   <  1:537536/150‑1 (MQ=255)
                     aCCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTttatta                                                                                                                                <  1:809863/150‑1 (MQ=255)
                     aCCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTttatta                                                                                                                                <  2:1199623/150‑1 (MQ=255)
                      ccATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCGGGATTCTGCCGCTTTTTtattat                                                                                                                               >  2:740792/1‑150 (MQ=255)
                      ccATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTtattat                                                                                                                               <  1:929109/150‑1 (MQ=255)
                       cATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTattatt                                                                                                                              >  2:654297/1‑150 (MQ=255)
                       cATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTattatt                                                                                                                              >  1:1196881/1‑150 (MQ=255)
                        aatattTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTAttt                                                                                                                             <  1:1120771/148‑1 (MQ=255)
                         ttattTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTAtttt                                                                                                                            <  2:921676/150‑1 (MQ=255)
                          tattTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTAttttt                                                                                                                           >  2:1675991/1‑150 (MQ=255)
                           attTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCTCAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACACGATTAACGGCAGAATCTCAGGATTCTGCCGGTATTTAATAATTTTTc                                                                                                                          >  1:286964/1‑150 (MQ=255)
                           attTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTAttttta                                                                                                                          >  2:1590479/1‑149 (MQ=255)
                            tttttAGAGGCCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCa                                                                                                                         <  1:1057586/150‑1 (MQ=255)
                                tAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAACAGCATTCAGAGATATATCATTTCTACAAACATCCCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAgggg                                                                                                                     >  1:121773/1‑150 (MQ=255)
                                     ggCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGagtag                                                                                                                <  2:407681/150‑1 (MQ=255)
                                      gCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGagtagt                                                                                                               >  2:1346404/1‑150 (MQ=255)
                                         tCTTGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATAAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAAt                                                                                                            <  1:350187/150‑1 (MQ=255)
                                         tCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAAt                                                                                                            >  2:1729208/1‑150 (MQ=255)
                                          cTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATg                                                                                                           <  2:1318889/150‑1 (MQ=255)
                                           tGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGa                                                                                                          >  1:1828358/1‑150 (MQ=255)
                                           tGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTATGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGa                                                                                                          >  2:1418334/1‑150 (MQ=255)
                                           tGGCACTACAACTTCCGCTTTTGGCACAGAACGGACTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTCAGGGGAGTAGTAATGAg                                                                                                         <  2:411583/150‑1 (MQ=255)
                                             gCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGc                                                                                                        <  1:1110776/150‑1 (MQ=255)
                                                cTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGCTTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGCGCAAc                                                                                                     >  2:266283/1‑150 (MQ=255)
                                                  ccAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGAAActc                                                                                                   >  2:942297/1‑150 (MQ=255)
                                                   cAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAActct                                                                                                  <  2:1522062/150‑1 (MQ=255)
                                                        tCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAgg                                                                                             >  1:1780632/1‑150 (MQ=255)
                                                         ccGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGt                                                                                            <  1:2922/150‑1 (MQ=255)
                                                           ccgtttGGCACAGAACGGCCTATCAGCAATATAAAAGCAAGCAGATATATATCATTTCTCCAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCGTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTg                                                                                          <  2:106508/147‑1 (MQ=255)
                                                            cTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTgg                                                                                         <  2:1127303/150‑1 (MQ=255)
                                                            cTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATAAATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTgg                                                                                         <  2:1927677/150‑1 (MQ=255)
                                                             ttttGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGt                                                                                        <  1:1913537/150‑1 (MQ=255)
                                                             ttttGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGt                                                                                        <  2:922314/150‑1 (MQ=255)
                                                               ttGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAg                                                                                      <  1:438737/150‑1 (MQ=255)
                                                                   cacaGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCAt                                                                                  >  2:1403087/1‑150 (MQ=255)
                                                                   cacaGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCAGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCAt                                                                                  >  2:837037/1‑150 (MQ=255)
                                                                     caGAACGGCCTATCAGCAATATAAAAGCATTCATATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCAttt                                                                                <  2:825759/150‑1 (MQ=255)
                                                                      aGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCAtttt                                                                               <  2:1563017/150‑1 (MQ=255)
                                                                      aGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCAtttt                                                                               <  1:1933887/150‑1 (MQ=255)
                                                                        aaCGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAg                                                                             <  1:1030449/150‑1 (MQ=255)
                                                                         aCGGCCTATCATCAATATAAAATCATTTAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCACATGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGc                                                                            <  2:329703/150‑1 (MQ=255)
                                                                           ggCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCaa                                                                          >  1:1039243/1‑150 (MQ=255)
                                                                           ggCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTAATTTTCAGGGTAGTAGTAATGAGCAACTCTCAAGGTTGGTAGTCATTTTAGCac                                                                          >  1:265469/1‑149 (MQ=255)
                                                                           ggCCTATCAGCAATATAAAAGCATTAAGATATATATCATTTCTACAAACATACCCGATTAACGTCAGAATCTCAGGATTCTGCCGCGTTTTTATTATTTTTCAGGGGAGTAGTAATTAGCAACTCTGAAGGTTGGTAGTCATTTTAGCaa                                                                          >  2:1199085/1‑150 (MQ=255)
                                                                            gCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCaaa                                                                         >  2:1945273/1‑150 (MQ=255)
                                                                            gCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCaaa                                                                         >  1:453845/1‑150 (MQ=255)
                                                                             ccTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAAc                                                                        <  1:1660524/150‑1 (MQ=255)
                                                                                aTCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCt                                                                     <  1:1571414/150‑1 (MQ=255)
                                                                                  cAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCttt                                                                   >  1:1113499/1‑150 (MQ=255)
                                                                                  cAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCttt                                                                   >  1:1121516/1‑150 (MQ=255)
                                                                                  cAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATAATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGCTAGTCATTATAGCAAACGCttt                                                                   >  2:243782/1‑150 (MQ=255)
                                                                                    gCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCt                                                                 <  2:1171762/150‑1 (MQ=255)
                                                                                        tataAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGta                                                             >  2:1151868/1‑149 (MQ=255)
                                                                                        tataAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTc                                                             <  1:1931437/150‑1 (MQ=255)
                                                                                         ataAAATCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCt                                                            <  2:1846703/150‑1 (MQ=255)
                                                                                            aaaGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGcc                                                         >  1:1162531/1‑150 (MQ=255)
                                                                                                          tatatCTTTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTTGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCgg                                           <  1:933076/150‑1 (MQ=255)
                                                                                                          tatatCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCgg                                           <  2:1034043/150‑1 (MQ=255)
                                                                                                             atCATTTCTACAAACATTCCCGATTAACGGCAGATTCTCAGGATTCTGCCGCTTTTTTTTTATTTTTCAGTGGAGTAGTAATGAGCAACTCTGAAGTTTGGTAGACAGTTTAGCAAACGCTTTCTGGTCTGCCACAGTGGGCATCGGCtt                                        >  2:57929/1‑149 (MQ=255)
                                                                                                               cATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGAt                                      >  2:51200/1‑150 (MQ=255)
                                                                                                               cATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCGTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGGCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGAt                                      <  1:132240/150‑1 (MQ=255)
                                                                                                                aTTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGAtt                                     <  1:654297/150‑1 (MQ=255)
                                                                                                                   tCTACAAACATACCCGATTATCGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCg                                  <  1:1345069/150‑1 (MQ=255)
                                                                                                                   tCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCg                                  <  1:1478530/150‑1 (MQ=255)
                                                                                                                    cTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTAAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGt                                 >  2:1697287/1‑150 (MQ=255)
                                                                                                                      aCAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAg                               >  2:471987/1‑150 (MQ=255)
                                                                                                                      aCAAACATACACGATTAACGGAAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTTCCGCAGTGGGCATCGGCTGATTGCGTat                               >  1:86444/1‑149 (MQ=255)
                                                                                                                       cAAACATACCCGATTAACGGCCGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGc                              <  1:1832352/150‑1 (MQ=255)
                                                                                                                       cAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGGGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGc                              <  1:2003736/150‑1 (MQ=255)
                                                                                                                       cAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGc                              >  1:1716381/1‑150 (MQ=255)
                                                                                                                             tACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGTGCATCGGCTGATTGCGTAGCAGGGcc                        >  1:1449383/1‑150 (MQ=255)
                                                                                                                             tACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGcc                        >  1:1444766/1‑150 (MQ=255)
                                                                                                                               cccGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCAc                      <  1:1665223/150‑1 (MQ=255)
                                                                                                                               cccGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATAGGCTGATTGCGTAGCAGGGCCAc                      >  2:51693/1‑150 (MQ=255)
                                                                                                                                 cGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTg                    >  2:2968/1‑150 (MQ=255)
                                                                                                                                  gATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTgg                   >  1:1258554/1‑150 (MQ=255)
                                                                                                                                        cGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTg             <  2:1593/150‑1 (MQ=255)
                                                                                                                                         ggCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCGGCAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGc            >  1:1465404/1‑150 (MQ=255)
                                                                                                                                         ggCAGAATCCTGAGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGc            >  1:628318/1‑150 (MQ=255)
                                                                                                                                         ggCAGAATCCTGAGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGc            >  1:1831802/1‑150 (MQ=255)
                                                                                                                                          gCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGcc           <  1:1513299/150‑1 (MQ=255)
                                                                                                                                            agtaTCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGATGGTTGGTAGTCTTTTTAGCAAACTCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATGGCGTAGCAGGGCCACTGGTATCTGCCgg         <  1:957658/147‑1 (MQ=255)
                                                                                                                                            aGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGCCgg         <  2:824524/150‑1 (MQ=255)
                                                                                                                                            aGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGCCgg         <  2:12606/150‑1 (MQ=255)
                                                                                                                                              aaTCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGCCGGTa       <  2:1751175/150‑1 (MQ=255)
                                                                                                                                                  tcAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGATTTCTGGTCTGCCGCAGTGGGCATCGGATGATTGCGTAGCAGGGCCACTGGTATCTGCCGGTATCAc   >  2:1261392/1‑150 (MQ=255)
                                                                                                                                                   cAGGATTCTGCCTCTTTTTAATTATATATCAGGGGAGTAGTACTGAGCAACGCTGAAGGGGGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATAGCGTAGCAGGGCCACTGGTATCTGCCGGTATCAcc  <  1:248777/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTAAATGCCCCTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGCCGGTATCACC  >  NC_000913/2775356‑2775652

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