Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,426,178 A→T 7.6% intergenic (+198/‑276) ydaY → / → ynaA Rac prophage; putative uncharacterized protein YdaY/Rac prophage; putative prophage tail length tape measure domain‑containing protein YnaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,426,1780AT 7.6% 183.3 / 6.3 67intergenic (+198/‑276)ydaY/ynaARac prophage; putative uncharacterized protein YdaY/Rac prophage; putative prophage tail length tape measure domain‑containing protein YnaA
Reads supporting (aligned to +/- strand):  ref base A (30/31);  new base T (2/3);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.70e-01

AAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGA  >  NC_000913/1426037‑1426316
                                                                                                                                              |                                                                                                                                          
aaGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGAc                                                                                                                                           >  1:97806/1‑143 (MQ=255)
   cGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACaga                                                                                                                                        >  1:57280/1‑143 (MQ=255)
    gagaAAATCTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACagag                                                                                                                                       <  1:167602/143‑1 (MQ=255)
    gagaAAAGCTCAGGTTCCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACagag                                                                                                                                       <  1:436800/143‑1 (MQ=255)
     agaAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAgaga                                                                                                                                      <  1:139164/143‑1 (MQ=255)
      gaAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTGAAGGAGTTGTTGTACGCCTATACGGTAATATCCCGTGCCCCACGGTATGCTGGGAAAGCTGAGGTGGCCTTGCCTTAAGCACGGACAGAGAt                                                                                                                                     >  2:569750/1‑143 (MQ=255)
          atcgcgggTCGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGTCAGAGATaaag                                                                                                                                 <  2:316496/137‑2 (MQ=255)
          aaCTCAGCGTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTGAAGGAGTTGTTGTACGCCTATTCGGGAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGTCAGAGATaaag                                                                                                                                 <  1:495486/143‑2 (MQ=255)
           aCTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATg                                                                                                                                >  1:545806/1‑143 (MQ=255)
           aCTCAGGGTGCGAGCGATACCGGAACCAGAAATGGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGTCAGAGATaaagg                                                                                                                                <  2:21693/143‑3 (MQ=255)
            cTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTACGGAGTTGTTGTACGCCTATTCGGTAATATCACGTGCCCGACGATACGCTGCGTTTGCTTGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGa                                                                                                                               >  1:506773/1‑143 (MQ=255)
               cggTGGCCACCGATACCTGAACCTGAAATAGCCCCGGTGCTTCAGTGCTTGCTCGGCGCCTCTTCGACAATAGTCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAata                                                                                                                            <  2:120458/142‑1 (MQ=255)
               aGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAGTGGCTGGCGTGCGTGTTCCTTTATCTCTGACAGAGCTAGATGAata                                                                                                                            >  1:153294/1‑143 (MQ=255)
                    gCCACGGATACCGGAACCAGAAATAGCCCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAg                                                                                                                       <  1:529767/143‑1 (MQ=255)
                     ccACCGAGACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGc                                                                                                                      <  2:563864/143‑1 (MQ=255)
                      cACCGATACCGGAACCAGAAATCGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGcc                                                                                                                     <  1:86904/143‑1 (MQ=255)
                       aCCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCa                                                                                                                    >  1:513878/1‑143 (MQ=255)
                       aCCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCa                                                                                                                    >  1:217939/1‑143 (MQ=255)
                           gtACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGAAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCGTTTATCTATGACAGAGCTAAATGAAAATTTAGCACCTCa                                                                                                                >  2:542889/2‑143 (MQ=255)
                                ggAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGCTAAATGAATATTTAGCCACTCATCCgg                                                                                                           >  1:461802/1‑143 (MQ=255)
                                ggAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCAATCATCCgg                                                                                                           >  1:284288/1‑143 (MQ=255)
                                    acagacagaGTCCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTatt                                                                                                       <  1:181113/142‑1 (MQ=255)
                                      aGAAATAGGCCCGGTGCTTAAGGAGGTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattga                                                                                                     <  2:505038/143‑1 (MQ=255)
                                          aTAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTAttgattga                                                                                                 <  1:282673/143‑1 (MQ=255)
                                           tAGAGCCGGTGGTTAAGGACGTGTGGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAg                                                                                                <  1:116764/143‑1 (MQ=255)
                                            aGACCCGGTGCGTAGGGAGTTGTGGGACGGCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATCGCTGGGGTGCCCGTGGAGAGATCTCTGACAGAGATAAATGAATATTGGGCCACTCATCCGGTATTGATTGAgc                                                                                               <  2:253886/143‑1 (MQ=255)
                                             gACCCGGTGCTTAAGGAGGTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAgcg                                                                                              <  1:100129/143‑1 (MQ=255)
                                                     ttttAAGGAGTTGTTGTACGCCTATTCGGTAATATGCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAAt                                                                                      <  1:437458/141‑1 (MQ=255)
                                                      cTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAAtt                                                                                     <  1:519916/143‑1 (MQ=255)
                                                              gttgttGTGCGCCTATTCGGTAATATCCCGTGCCCGACTTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATGAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAg                                                                             <  2:384486/143‑1 (MQ=255)
                                                                 gttgtACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGa                                                                          >  2:78994/1‑143 (MQ=255)
                                                                 ggtgtACGCCTAGTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGTCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGa                                                                          <  2:50800/141‑1 (MQ=255)
                                                                  ttgtACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGAt                                                                         >  2:537292/1‑143 (MQ=255)
                                                                  ttgtACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGAGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGAt                                                                         >  1:433560/1‑143 (MQ=255)
                                                                  ttgtACGCCTATTCGGTAATATCCAGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCGTTGCCTTTGTCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGGATTGAATGAGCGCGATGAATTTGAAGCAGTGAt                                                                         >  2:121740/1‑143 (MQ=255)
                                                                      aggCCTATTCGGTAATATCCCGTGCCCGACGTTATGCGGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAGATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCttt                                                                     <  1:495493/141‑1 (MQ=255)
                                                                       cGCCTATTGGGTAATAGCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACATAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTg                                                                    <  2:144844/143‑1 (MQ=255)
                                                                        gCCTATTCGGTAATATCCCGTGCACAAAGATATACTGGACTGGCTGGGGTGCCTTTGCCTTTATCTATGACAGAGATAAACGAATAGTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAAGTTGACACATTGATCTTAGc                                                                   >  2:45106/1‑143 (MQ=255)
                                                                         ccTATTCGGTAATATCCCTTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTGTGATAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCa                                                                  <  2:316123/143‑1 (MQ=255)
                                                                             ttCGGTAACATCCCGTGCCCGACGTTATGCTGGACCGGCTGGGGTGCCTTTGCCTTTATCTCTGCACAAGTGGAAGTAGTATTTTGTCTCTCATCCGGTATTGATTGAGCGCAATGAATTTGAAGCAGTGATCTTTGCACTgg                                                              <  2:554045/143‑1 (MQ=255)
                                                                               cGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGAt                                                            <  1:8434/143‑1 (MQ=255)
                                                                                   aagatCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGAcc                                                        <  1:246797/140‑1 (MQ=255)
                                                                                    atatCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAAGATTTAGCCACTCATCCGGTATTGCTTGAGCGCGATGAAGTTGACGCAGGTATCGTTGCACAGGATGAccc                                                       >  2:561408/1‑142 (MQ=255)
                                                                                      atCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGt                                                     <  1:72166/143‑1 (MQ=255)
                                                                                      atCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGt                                                     >  1:177804/1‑143 (MQ=255)
                                                                                         ccGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTAtt                                                  >  1:455351/1‑143 (MQ=255)
                                                                                         ccGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTAtt                                                  >  1:270913/1‑143 (MQ=255)
                                                                                            tGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTc                                               >  1:329882/1‑143 (MQ=255)
                                                                                               ccGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGCTGAATTTGAAGCAGTGATCTTTGCACTGGGTGACCCCTATTTTCAgg                                            >  1:133704/1‑143 (MQ=255)
                                                                                                  aCGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGc                                         >  1:435840/1‑143 (MQ=255)
                                                                                                         gCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCCgtgtgt                                  >  1:79287/1‑143 (MQ=255)
                                                                                                          cTGGAATGGCTGGGGTGCCTTTGCCTTTAGCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAgtgtgtg                                 <  1:233931/143‑1 (MQ=255)
                                                                                                             gAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAg                              >  2:233479/1‑143 (MQ=255)
                                                                                                                 ggCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAgttgt                          <  1:303299/143‑1 (MQ=255)
                                                                                                                   cgggggTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTCGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATCTTCAGGAGCAGTGTGTGTAgtcgata                        >  2:234593/3‑138 (MQ=255)
                                                                                                                   cTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTa                        <  1:371130/143‑1 (MQ=255)
                                                                                                                    tGGGGTGCCGTGGCCTTGTTCTCCGCCTGTGCTAAATGAGTATTTAGCGTTTCATGCGGTGTTGATTGATAGCGATGAATTTGAAGCAGTGTGTTTTGCACTGGATGCCCCGGATTTTCAGGTGCAGTGTGTGTAGTTGTTaa                       <  2:203461/143‑1 (MQ=255)
                                                                                                                            ctttgcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGGGCGATGAAGTGGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGtgc               >  1:213508/1‑141 (MQ=255)
                                                                                                                              ttgcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTacac             >  2:534810/1‑143 (MQ=255)
                                                                                                                               ttcatttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACAct            >  1:186826/5‑143 (MQ=255)
                                                                                                                                 cctttATCTCTGACAGCGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACActct          <  2:513878/143‑1 (MQ=255)
                                                                                                                                 cctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACActct          <  2:217939/143‑1 (MQ=255)
                                                                                                                                  ctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTCGTTGTTAATTACGTACACTCTg         >  2:197475/1‑143 (MQ=255)
                                                                                                                                  atttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTg         >  2:335911/2‑143 (MQ=255)
                                                                                                                                   tttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGt        >  1:156197/1‑143 (MQ=255)
                                                                                                                                     tagccctGACAGAGATAAATAGATATTTAGCCACGCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTa      <  1:165019/138‑1 (MQ=255)
                                                                                                                                         tctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACaga  <  2:433560/143‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                          
AAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGA  >  NC_000913/1426037‑1426316

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