Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,775,502 T→C 5.0% intergenic (+481/‑43) 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,5020TC 5.0% 425.1 / 8.8 140intergenic (+481/‑43)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 T (63/69);  new base C (2/5);  total (65/75)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

Only 100 of 143 total aligned reads displayed.

ACTGTAAATGCCCCTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGCCGGTATCAC  >  NC_000913/2775353‑2775651
                                                                                                                                                     |                                                                                                                                                     
aCTGTAAATGCCCCTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAAtct                                                                                                                                                       <  1:1841564/150‑1 (MQ=255)
aCTGTAAATGCCCCTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAAtct                                                                                                                                                       >  2:63678/1‑150 (MQ=255)
   gTAAATGCCCCTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTACGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCTATTAACGGCAGAATCTCAg                                                                                                                                                    >  2:949438/1‑150 (MQ=255)
         gCCCCTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCt                                                                                                                                              >  1:1171762/1‑150 (MQ=255)
            ccTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCCTGAGATTCTGcc                                                                                                                                           <  2:69730/150‑1 (MQ=255)
            ccTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCCTGAGATTCTGcc                                                                                                                                           <  2:774775/150‑1 (MQ=255)
             cTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAGCTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCg                                                                                                                                          >  1:172929/1‑150 (MQ=255)
              tACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGc                                                                                                                                         >  2:1083372/1‑150 (MQ=255)
              tACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCCTGAGATTCTGCCgt                                                                                                                                         <  1:1647005/150‑2 (MQ=255)
               aCAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTTGCACAGAACGGCATATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCt                                                                                                                                        <  1:1508116/150‑1 (MQ=255)
               aCAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCCTGAGATTCTGCCgtt                                                                                                                                        <  2:431497/150‑3 (MQ=255)
                cAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGCACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACACGATTAACGGCAGAATCTCAGGATTCTGCCGCtt                                                                                                                                       >  1:1800626/1‑150 (MQ=255)
                 aaaaTGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCTTTGAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCttt                                                                                                                                      <  2:1235340/150‑1 (MQ=255)
                  aaaTGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCtttt                                                                                                                                     >  2:106891/1‑150 (MQ=255)
                  aaaTGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAAAGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCtttt                                                                                                                                     >  1:1967604/1‑150 (MQ=255)
                     tGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTtta                                                                                                                                  <  1:537536/150‑1 (MQ=255)
                        aCCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTttatta                                                                                                                               <  1:809863/150‑1 (MQ=255)
                        aCCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTttatta                                                                                                                               <  1:799044/150‑1 (MQ=255)
                         ccATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCGGGATTCTGCCGCTTTTTtattat                                                                                                                              >  2:740792/1‑150 (MQ=255)
                         ccATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTtattat                                                                                                                              <  1:929109/150‑1 (MQ=255)
                          cATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTattatt                                                                                                                             >  2:407956/1‑150 (MQ=255)
                          cATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTattatt                                                                                                                             >  1:1196881/1‑150 (MQ=255)
                          cATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTattatt                                                                                                                             >  2:654297/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)
                              attTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCTCAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACACGATTAACGGCAGAATCTCAGGATTCTGCCGGTATTTAATAATTTTTc                                                                                                                         >  1:286964/1‑150 (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)
                                             cTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGTATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATg                                                                                                          <  1:1831932/150‑1 (MQ=255)
                                             cTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATg                                                                                                          <  2:1318889/150‑1 (MQ=255)
                                              tGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGa                                                                                                         >  1:1828358/1‑150 (MQ=255)
                                               ggCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACTTACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAg                                                                                                        <  2:1375763/150‑1 (MQ=255)
                                                gCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGc                                                                                                       <  1:1110776/150‑1 (MQ=255)
                                                   cTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAAc                                                                                                    >  2:1180285/1‑150 (MQ=255)
                                                     ccAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGAAActc                                                                                                  >  2:942297/1‑150 (MQ=255)
                                                      cAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAActct                                                                                                 <  2:1522062/150‑1 (MQ=255)
                                                       aaCTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTg                                                                                                <  2:1507406/150‑1 (MQ=255)
                                                            ccGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGt                                                                                           >  2:28669/1‑150 (MQ=255)
                                                               cTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTgg                                                                                        <  2:1127303/150‑1 (MQ=255)
                                                               cTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTgg                                                                                        >  2:1662735/1‑150 (MQ=255)
                                                                ttttGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGt                                                                                       <  1:1913537/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)
                                                                         aGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCAtttt                                                                              <  2:1563017/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)
                                                                               gCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCaaa                                                                        >  2:1945273/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)
                                                                                      aGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCttta                                                                 >  1:533230/1‑149 (MQ=255)
                                                                                       gCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCt                                                                <  2:1171762/150‑1 (MQ=255)
                                                                                       gCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCt                                                                <  2:451854/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)
                                                                                                             tatatCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCgg                                          >  1:1451910/1‑150 (MQ=255)
                                                                                                             tatatCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCgg                                          <  2:1034043/150‑1 (MQ=255)
                                                                                                                  cATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCGTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGGCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGAt                                     <  1:132240/150‑1 (MQ=255)
                                                                                                                      tCTACAAACATACCCGATTATCGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCg                                 <  1:1345069/150‑1 (MQ=255)
                                                                                                                      tCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCg                                 <  1:1478530/150‑1 (MQ=255)
                                                                                                                        tACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTa                               >  2:112029/1‑150 (MQ=255)
                                                                                                                         aCAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAg                              >  2:471987/1‑150 (MQ=255)
                                                                                                                          cAAACATACCCGATTAACGGCCGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGc                             <  1:1832352/150‑1 (MQ=255)
                                                                                                                          cAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGc                             >  1:1716381/1‑150 (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)
                                                                                                                                           cGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTg            <  2:1593/150‑1 (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:628318/1‑150 (MQ=255)
                                                                                                                                             gCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGcc          <  1:1513299/150‑1 (MQ=255)
                                                                                                                                              cAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGCCg         <  2:1307358/150‑1 (MQ=255)
                                                                                                                                               agtaTCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGATGGTTGGTAGTCTTTTTAGCAAACTCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATGGCGTAGCAGGGCCACTGGTATCTGCCgg        <  1:957658/147‑1 (MQ=255)
                                                                                                                                               aGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGCCgg        <  2:12606/150‑1 (MQ=255)
                                                                                                                                               aGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGCCgg        <  2:523850/150‑1 (MQ=255)
                                                                                                                                               aGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGCCgg        <  1:885527/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACTGTAAATGCCCCTACAAAATGAACCATTATTTTTAGAGGGCGTCTGGCACTCCAACTTCCGCTTTTGGCACAGAACGGCCTATCAGCAATATAAAAGCATTCAGATATATATCATTTCTACAAACATACCCGATTAACGGCAGAATCTCAGGATTCTGCCGCTTTTTTATTATTTTTCAGGGGAGTAGTAATGAGCAACTCTGAAGGTTGGTAGTCATTTTAGCAAACGCTTTCTGGTCTGCCGCAGTGGGCATCGGCTGATTGCGTAGCAGGGCCACTGGTATCTGCCGGTATCAC  >  NC_000913/2775353‑2775651

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