Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,426,178 A→T 7.2% intergenic (+96/‑276) ydaY → / → tmpR pseudogene, Rac prophage;Phage or Prophage Related/Rac prophage; pseudogene, tail protein family;Phage or Prophage Related; putative alpha helix protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,426,1780AT 7.2% 181.9 / 3.0 66intergenic (+96/‑276)ydaY/tmpRpseudogene, Rac prophage;Phage or Prophage Related/Rac prophage; pseudogene, tail protein family;Phage or Prophage Related; putative alpha helix protein
Reads supporting (aligned to +/- strand):  ref base A (28/33);  new base T (3/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

CAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGAT  >  NC_000913/1426029‑1426324
                                                                                                                                                     |                                                                                                                                                  
cAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGa                                                                                                                                                    <  2:850805/150‑1 (MQ=255)
       gAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGCCCCGGTGCTTAAGGAGTTGTTGTACGCCTACTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGAt                                                                                                                                             >  2:331631/1‑150 (MQ=255)
          gCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATaaa                                                                                                                                          >  1:526702/1‑150 (MQ=255)
            gagaAAAACTTAGGTTGAGAGAGATACCTGAACCAGAAATATACCCGGTGCTTAACGAGTTGTTGTACGCCTTTTCGTTAATTACCCTTGCCCGACGTTATGCTGTAATGGCGGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATg                                                                                                                                        <  2:554653/150‑1 (MQ=255)
            gagaAAAACTCTGGGTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTGGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGTCAGAGATaaagg                                                                                                                                        <  2:1241002/150‑3 (MQ=255)
            gagaAAAACTCAGGTTGCCACCGATACCGGATCCAGAAATAGACCCGGTGCTTAAGGAGGTGTTGTACGCCTATTCGGTAAGATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATg                                                                                                                                        <  1:851125/150‑1 (MQ=255)
            gagaAAAACTCAGGTGGCCCCCGATACCGCAACCAGAAATAGACCCGGTGCTCCAGGAGTTGTTGTCCGCCTATTCGGTAATATCCCGTGCCCGACGATATGCTGGAATGGCGGGGGTGCCTTTGCCTTTATCTCTGTCAGAGATaaagg                                                                                                                                        <  1:618475/150‑3 (MQ=255)
              gaAAAACTCAGGTTGCCACCGATACCGGAACCAGAATTAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATCGCTGGGGTGCCTTTGCCTTTATCTCTGTCAGAGATAAATGaa                                                                                                                                      >  2:1007895/1‑150 (MQ=255)
               aaacaCTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGGGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCGTTGCCTTTATCTCTGACAGAGATAAATGAat                                                                                                                                     <  1:900248/150‑1 (MQ=255)
               aaaaaCTCAGGTTTCCACCGATACCGGAACCAGAAATACACCCGGTGCTTAAGGATTTGTTGTACGCCTATTCGGTAAAATCCCGTTCCCGACGTTAAGCTGGAATGGCTTGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAat                                                                                                                                     <  1:759562/150‑1 (MQ=255)
               aaaaaCTCAGGTTGCCAACGATACCGTAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATTCTTGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAat                                                                                                                                     <  1:887583/150‑1 (MQ=255)
                aaaaCTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAata                                                                                                                                    <  2:1075438/150‑1 (MQ=255)
                    cTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTa                                                                                                                                >  2:1262340/1‑150 (MQ=255)
                       cggTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGcc                                                                                                                             <  1:677024/149‑1 (MQ=255)
                           tGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACtt                                                                                                                         <  2:739243/150‑2 (MQ=255)
                                  gATACCGGAACCAGAAATATCCCCGGGGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTTTGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGt                                                                                                                  <  1:678439/150‑1 (MQ=255)
                                   aTACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTa                                                                                                                 <  1:950615/150‑1 (MQ=255)
                                      ccGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattg                                                                                                              >  1:912922/1‑150 (MQ=255)
                                        ggAACCAGCAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATACCGTTCCCGACGTTATGCTGGAATGTCTGGGGTTCCTTTGCCTTTAGCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattgat                                                                                                            <  1:197213/150‑1 (MQ=255)
                                        ggAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattgat                                                                                                            >  1:483643/1‑150 (MQ=255)
                                         gAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattgatt                                                                                                           <  1:120337/150‑1 (MQ=255)
                                          aaCCAGAAATAGACCCGGTGCTTAATGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattgattg                                                                                                          <  1:798452/150‑1 (MQ=255)
                                             cAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAgc                                                                                                       >  1:305950/1‑150 (MQ=255)
                                               gAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAgcgc                                                                                                     <  1:1247929/150‑1 (MQ=255)
                                                  aTAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATGACGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGAt                                                                                                  <  1:458111/150‑1 (MQ=255)
                                                  aTAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGAt                                                                                                  >  1:588644/1‑150 (MQ=255)
                                                       cccGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAAtt                                                                                             <  2:526702/150‑1 (MQ=255)
                                                           gTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGaa                                                                                         >  1:1238353/1‑150 (MQ=255)
                                                             gCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGc                                                                                       <  2:610531/150‑1 (MQ=255)
                                                                    gatttgttgTTCGCATATTCGGTAATATCCCGTGCCCTACGTTATGCTGGAATTGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATc                                                                                <  1:173359/147‑1 (MQ=255)
                                                                    gAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATc                                                                                <  2:596371/150‑1 (MQ=255)
                                                                      gttgttGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCtt                                                                              <  1:393405/150‑1 (MQ=255)
                                                                       ttgttgTACGCCTATTCGGTAATATCCCGTGGCCGACGTTATGCTGTAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACTTAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCttt                                                                             <  1:795077/150‑1 (MQ=255)
                                                                       ttgttgTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGGAGCAGTGATCttt                                                                             <  2:1146586/150‑1 (MQ=255)
                                                                            gtACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCaat                                                                        >  1:703416/1‑148 (MQ=255)
                                                                            gtACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACt                                                                        >  1:409280/1‑150 (MQ=255)
                                                                             tACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTg                                                                       <  1:1092540/150‑1 (MQ=255)
                                                                                 ccTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGAtt                                                                   >  1:502690/1‑149 (MQ=255)
                                                                                 ccTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATg                                                                   >  1:1292252/1‑150 (MQ=255)
                                                                                    aTTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGAcc                                                                <  1:748487/150‑1 (MQ=255)
                                                                                       cGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGt                                                             >  2:1286898/1‑150 (MQ=255)
                                                                                        ggTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTa                                                            >  1:1003889/1‑150 (MQ=255)
                                                                                        ggTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTa                                                            >  1:377099/1‑150 (MQ=255)
                                                                                               tCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAg                                                     <  2:1298059/150‑1 (MQ=255)
                                                                                                     gCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTAAGCCACTCAACCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAg                                               >  2:279275/1‑150 (MQ=255)
                                                                                                      cccGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGAAgt                                              >  1:136452/1‑150 (MQ=255)
                                                                                                        cGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAgtgt                                            >  1:237063/1‑150 (MQ=255)
                                                                                                              tATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAg                                      <  2:315301/150‑1 (MQ=255)
                                                                                                                 gCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAgttg                                   <  1:259283/150‑1 (MQ=255)
                                                                                                                 gCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACGGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAgttg                                   <  1:303563/150‑1 (MQ=255)
                                                                                                                    ggAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTa                                >  1:334609/1‑150 (MQ=255)
                                                                                                                     gAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTaa                               <  2:409280/150‑1 (MQ=255)
                                                                                                                      aaTGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAAt                              >  1:323165/1‑150 (MQ=255)
                                                                                                                           cTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGt                         <  1:434567/150‑1 (MQ=255)
                                                                                                                             ggggTGCCTTTGCCTTAATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTTAGCGCGATGAATTTGAAGCAGTTATCTTTGCAATGGATGACCAGTATTTTCATGAGCAGTGTGTGTAGTTGTTAATTACGTac                       <  1:615490/150‑1 (MQ=255)
                                                                                                                               ggTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTacaa                     >  1:836755/1‑149 (MQ=255)
                                                                                                                                    ctttgcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGt                >  1:312578/1‑150 (MQ=255)
                                                                                                                                    ctttgcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGt                >  2:846850/1‑150 (MQ=255)
                                                                                                                                     tttgcctttTTCTCTGACAGAGATAAATGAATATTTATCCACTCATCCTGTAATTATTGAGCGCGATGAATTTGAAGCCGTGAGCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTTTAGTTGTTAATTACGTACACTCTGtt               <  1:533848/150‑1 (MQ=255)
                                                                                                                                     tattcatttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGGTTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGtt               >  2:707844/7‑150 (MQ=255)
                                                                                                                                       tgcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTAc             >  2:983715/1‑150 (MQ=255)
                                                                                                                                         cctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACag           >  2:986982/1‑150 (MQ=255)
                                                                                                                                          atttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACaga          >  1:738385/2‑150 (MQ=255)
                                                                                                                                            ttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATACGTTATTGATTGAGCGCGATGAATTTGAAGCAATGATCATTCCACTGCAAGCCCAGTATTTTCCGGAGACGTGTGTGTAGTTGTGAATTAAGTACACTCGTTTACAtaga        >  2:784219/1‑150 (MQ=255)
                                                                                                                                             tATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGAt       >  2:890939/1‑150 (MQ=255)
                                                                                                                                               tctcTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGAtgt     <  2:912922/150‑1 (MQ=255)
                                                                                                                                                  ctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGAt  >  1:816078/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGAT  >  NC_000913/1426029‑1426324

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: