Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,775,505 G→A 5.1% intergenic (+484/‑40) 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,5050GA 5.1% 409.8 / 6.9 140intergenic (+484/‑40)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 G (62/69);  new base A (2/5);  total (64/76)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.61e-01

Only 100 of 146 total aligned reads displayed.

GTAAATGCCCCTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGCCGGTATCACCGA  >  NC_000913/2775356‑2775654
                                                                                                                                                     |                                                                                                                                                     
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)
      gCCCCTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCt                                                                                                                                                 >  1:1171762/1‑150 (MQ=255)
         ccTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCCTGAGATTCTGcc                                                                                                                                              <  2:67812/150‑1 (MQ=255)
         ccTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCCTGAGATTCTGcc                                                                                                                                              <  2:774775/150‑1 (MQ=255)
         ccTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCCTGAGATTCTGcc                                                                                                                                              <  2:69730/150‑1 (MQ=255)
          cTACATAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCTCATTTGGCACAGAACGGACTTTCAGCAATATTAAAGCATTCAGATATATATCTTTTCTACAATCATACCCGATTAACGGCAGAATCTAAGGATTCTGCCg                                                                                                                                             >  1:1847664/1‑150 (MQ=255)
          cTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAGCTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCg                                                                                                                                             >  1:172929/1‑150 (MQ=255)
           tACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCCTGAGATTCTGCCgt                                                                                                                                            <  1:1647005/150‑2 (MQ=255)
            aCAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTTGCACAGAACGGCATATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCt                                                                                                                                           <  1:1508116/150‑1 (MQ=255)
             cAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGCACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACACGATTAACGGCAGAATCTCAGGATTCTGCCGCtt                                                                                                                                          >  1:1800626/1‑150 (MQ=255)
              aaaaTGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCttt                                                                                                                                         >  2:1090534/1‑150 (MQ=255)
               aaaTGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCtttt                                                                                                                                        >  2:106891/1‑150 (MQ=255)
                 aTGAACCATTATTTTTAGAGGGAGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACAGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCtttttt                                                                                                                                      <  2:152363/150‑1 (MQ=255)
                  tGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTtta                                                                                                                                     <  1:537536/150‑1 (MQ=255)
                     aCCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTttatta                                                                                                                                  <  2:1199623/150‑1 (MQ=255)
                     aCCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTttatta                                                                                                                                  <  1:799044/150‑1 (MQ=255)
                       cATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTattatt                                                                                                                                >  1:1196881/1‑150 (MQ=255)
                       cATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTattatt                                                                                                                                >  2:407956/1‑150 (MQ=255)
                        attattTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTAttt                                                                                                                               <  2:1580390/150‑1 (MQ=255)
                        aatattTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTAttt                                                                                                                               <  1:1120771/148‑1 (MQ=255)
                         ttattTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTAtttt                                                                                                                              <  2:921676/150‑1 (MQ=255)
                          tattTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTAttttt                                                                                                                             >  2:1675991/1‑150 (MQ=255)
                           attTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTAttttta                                                                                                                            >  2:1590479/1‑149 (MQ=255)
                            tttttAGAGGCCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCa                                                                                                                           <  1:1057586/150‑1 (MQ=255)
                              tttAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATAACCGATTAACGTCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAgg                                                                                                                         >  1:1237294/1‑150 (MQ=255)
                                  gagGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGag                                                                                                                     <  1:1699328/150‑1 (MQ=255)
                                  gagGGAGTCTGGCACTCCAACTTCCGCTTTTGGAACAGAACGTCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTCACGGCAGAATCTCAGGATTCTGCCGCGTTTTTATTATTTTTCAGGGGag                                                                                                                     <  1:1498616/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)
                                          cTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGTATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATg                                                                                                             <  1:1831932/150‑1 (MQ=255)
                                          cTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATg                                                                                                             <  2:1318889/150‑1 (MQ=255)
                                           tGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTATGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGa                                                                                                            >  2:1418334/1‑150 (MQ=255)
                                           tGGCACTACAACTTCCGCTTTTGGCACAGAACGGACTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTCAGGGGAGTAGTAATGAg                                                                                                           <  2:411583/150‑1 (MQ=255)
                                             gCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGc                                                                                                          <  1:1110776/150‑1 (MQ=255)
                                                  ccAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGAAActc                                                                                                     >  2:942297/1‑150 (MQ=255)
                                                   cAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAActct                                                                                                    <  2:1522062/150‑1 (MQ=255)
                                                      cTTACGCTTTTTGAACAGAACGGCCTATCAGCAATATAATATCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCATAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGTGAGTAGTAATGAGCAACTCTGaa                                                                                                 <  1:132426/150‑1 (MQ=255)
                                                        tCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAgg                                                                                               >  1:1780632/1‑150 (MQ=255)
                                                         ccGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGt                                                                                              >  2:28669/1‑150 (MQ=255)
                                                         ccGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGt                                                                                              <  1:2922/150‑1 (MQ=255)
                                                           ccgtttGGCACAGAACGGCCTATCAGCAATATAAAAGCAAGCAGATATATATCATTTCTCCAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCGTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTg                                                                                            <  2:106508/147‑1 (MQ=255)
                                                            cTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTgg                                                                                           >  2:1662735/1‑150 (MQ=255)
                                                            cTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATAAATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTgg                                                                                           <  2:1927677/150‑1 (MQ=255)
                                                             ttttGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGt                                                                                          <  2:922314/150‑1 (MQ=255)
                                                               ttGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAg                                                                                        <  1:438737/150‑1 (MQ=255)
                                                                 ggCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTc                                                                                      >  1:1527798/1‑150 (MQ=255)
                                                                   cacaGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCAt                                                                                    >  2:1403087/1‑150 (MQ=255)
                                                                   cacaGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCAGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCAt                                                                                    >  2:837037/1‑150 (MQ=255)
                                                                     caGAACGGCCTATCAGCAATATAAAAGCATTCATATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCAttt                                                                                  <  2:825759/150‑1 (MQ=255)
                                                                      aGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCAtttt                                                                                 <  1:1933887/150‑1 (MQ=255)
                                                                        aaCGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAg                                                                               <  1:1030449/150‑1 (MQ=255)
                                                                           ggCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCaa                                                                            >  1:1039243/1‑150 (MQ=255)
                                                                           ggCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCaa                                                                            >  1:249431/1‑150 (MQ=255)
                                                                           ggCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTAATTTTCAGGGTAGTAGTAATGAGCAACTCTCAAGGTTGGTAGTCATTTTAGCac                                                                            >  1:265469/1‑149 (MQ=255)
                                                                           ggCCTATCAGCAATATAAAAGCATTAAGATATATATCATTTCTACAAACATACCCGATTAACGTCAGAATCTCAGGATTCTGCCGCGTTTTTATTATTTTTCAGGGGAGTAGTAATTAGCAACTCTGAAGGTTGGTAGTCATTTTAGCaa                                                                            >  2:1199085/1‑150 (MQ=255)
                                                                                  cAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCttt                                                                     >  1:1113499/1‑150 (MQ=255)
                                                                                  cAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCttt                                                                     >  1:1121516/1‑150 (MQ=255)
                                                                                  cAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATAATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGCTAGTCATTATAGCAAACGCttt                                                                     >  2:243782/1‑150 (MQ=255)
                                                                                   aGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCttta                                                                    >  1:533230/1‑149 (MQ=255)
                                                                                    gCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCt                                                                   <  2:1171762/150‑1 (MQ=255)
                                                                                    gCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCt                                                                   <  2:451854/150‑1 (MQ=255)
                                                                                        tataAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGta                                                               >  2:1151868/1‑149 (MQ=255)
                                                                                        tataAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTc                                                               <  1:1931437/150‑1 (MQ=255)
                                                                                        tataAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTc                                                               <  2:1716017/150‑1 (MQ=255)
                                                                                         ataAAATCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCt                                                              <  2:1846703/150‑1 (MQ=255)
                                                                                              aGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGGCTGCCGc                                                         <  1:949438/150‑1 (MQ=255)
                                                                                                cATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAg                                                       <  1:357889/150‑1 (MQ=255)
                                                                                                         atatatCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTCAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCg                                              <  2:2014020/150‑1 (MQ=255)
                                                                                                          tatatCTTTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTTGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCgg                                             <  1:933076/150‑1 (MQ=255)
                                                                                                          tatatCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCgg                                             <  2:1034043/150‑1 (MQ=255)
                                                                                                          tatatCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCgg                                             >  1:1451910/1‑150 (MQ=255)
                                                                                                             atCATTTCTACAAACATTCCCGATTAACGGCAGATTCTCAGGATTCTGCCGCTTTTTTTTTATTTTTCAGTGGAGTAGTAATGAGCAACTCTGAAGTTTGGTAGACAGTTTAGCAAACGCTTTCTGGTCTGCCACAGTGGGCATCGGCtt                                          >  2:57929/1‑149 (MQ=255)
                                                                                                                aTTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGAtt                                       <  1:654297/150‑1 (MQ=255)
                                                                                                                   tCTACAAACATACCCGATTATCGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCg                                    <  1:1345069/150‑1 (MQ=255)
                                                                                                                     tACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTa                                  >  2:112029/1‑150 (MQ=255)
                                                                                                                      aCAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAg                                 >  2:471987/1‑150 (MQ=255)
                                                                                                                      aCAAACATACACGATTAACGGAAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTTCCGCAGTGGGCATCGGCTGATTGCGTat                                 >  1:86444/1‑149 (MQ=255)
                                                                                                                       cAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGGGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGc                                <  1:2003736/150‑1 (MQ=255)
                                                                                                                             tACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGcc                          >  1:1444766/1‑150 (MQ=255)
                                                                                                                              aCCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCa                         >  1:1873066/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)
                                                                                                                                         ggCAGAATCTCTCGATCCTGCCGCTTTTTTATTATTTGTCAGGCCAGAAGTAATCAGCAACTCTGAAGGTTGGTAGTCATTTTAGCCAACGCTTTCTGGTCAGCCGCCGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGc              <  1:341443/150‑1 (MQ=255)
                                                                                                                                         ggCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGc              >  2:189799/1‑150 (MQ=255)
                                                                                                                                         ggCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCGGCAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGc              >  1:1465404/1‑150 (MQ=255)
                                                                                                                                         ggCAGAATCCTGAGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGc              >  1:1831802/1‑150 (MQ=255)
                                                                                                                                            aGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGCCgg           <  2:523850/150‑1 (MQ=255)
                                                                                                                                            aGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGCCgg           <  2:12606/150‑1 (MQ=255)
                                                                                                                                            aGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGCCgg           <  1:885527/150‑1 (MQ=255)
                                                                                                                                                  tcAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGATTTCTGGTCTGCCGCAGTGGGCATCGGATGATTGCGTAGCAGGGCCACTGGTATCTGCCGGTATCAc     >  2:1261392/1‑150 (MQ=255)
                                                                                                                                                    aGGATTCTGCCGCTTTTTAATTATTTATCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGCCGGTATCACCg   <  2:388030/150‑1 (MQ=255)
                                                                                                                                                    aGGATGCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACGCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGCCGGTATCACCg   <  2:787935/150‑1 (MQ=255)
                                                                                                                                                     ggATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTTGTCATTTTAGCAAAGGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGCCGGTATCACCGa  >  2:343311/1‑150 (MQ=255)
                                                                                                                                                     ggATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGATTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGCCGGTATCACCGa  >  1:244540/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTAAATGCCCCTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGCCGGTATCACCGA  >  NC_000913/2775356‑2775654

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: