Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,563,518:1 (A)5→6 100% intergenic (+401/‑74) yffP → / → yffQ CPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,563,5131.A94.4% 400.1 / 16.7 124intergenic (+396/‑79)yffP/yffQCPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base . (3/4);  new base A (62/54);  total (66/58)
Fisher's exact test for biased strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 135 total aligned reads displayed.

TGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGT‑AA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAATCG  >  NC_000913/2563417‑2563607
                                                                                                  |                                                                                               
tGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAAAGAAAGGGAGTTCCCAACAACAAAATAAAGATGGGA‑GGG‑aa                                                                                               >  1:2704925/1‑99 (MQ=255)
   aTTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGT‑aa‑aaaa                                                                                          >  1:3342675/1‑99 (MQ=255)
   aTTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGT‑aa‑aaaa                                                                                          <  2:4078917/100‑2 (MQ=255)
   aTTTTGTGATGCAGATCGATTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAAAAAAAGGAAGTTCCCAACAACAAGATAAAGATGGGA‑GGT‑aa‑gaaa                                                                                          >  2:1863587/1‑99 (MQ=255)
     tttGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGT‑aa‑aaaatg                                                                                        <  2:113590/100‑4 (MQ=255)
     tttGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGT‑aa‑aaaatg                                                                                        <  1:2968985/100‑4 (MQ=255)
     tttGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGT‑aa‑aaaatg                                                                                        <  1:1982230/100‑4 (MQ=255)
       tgtgATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGT‑aa‑aaaatgaa                                                                                      <  1:1695034/100‑6 (MQ=255)
       tgtgATGCAGATCGCTTTTTCCCATCAGGCATTTTTGTTTAATTAAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGT‑aa‑aaaatgaa                                                                                      <  2:726328/100‑6 (MQ=255)
        gtgATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGaaa                                                                                      <  1:3346365/100‑1 (MQ=255)
         ttaTGCAGTTCGCTTTTTCCCTTCAAGCATTTTTAATTAAGTCAATCAACGAAAGGCAGTTCCCAACAACAAATTAAAGATTGGA‑GGT‑aa‑aaaatgagac                                                                                    <  2:2002594/98‑8 (MQ=255)
          gATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACg                                                                                    >  1:3166429/1‑100 (MQ=255)
            tGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGt                                                                                  >  1:3871438/1‑100 (MQ=255)
             gCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTg                                                                                 <  1:2308846/100‑1 (MQ=255)
              cAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAATAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTgg                                                                                >  2:1391215/1‑100 (MQ=255)
              cAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTgg                                                                                <  2:1168555/100‑1 (MQ=255)
              cAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTgg                                                                                >  1:2578695/1‑100 (MQ=255)
               aGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTggg                                                                               >  1:654837/1‑100 (MQ=255)
               aGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAACTCAATCAACGAAAGGAAGTTCCCACCAACAAAATAAAGATGGGAGGGA‑AA‑AAATGAAACGGTggg                                                                               >  1:1344805/1‑100 (MQ=255)
                gATCGCTTTTTCCTTTTATGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGa                                                                              >  2:142931/1‑100 (MQ=255)
                gATCGCTTTTTCCCTTGAGGCATTTTTGTTTAATGCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGa                                                                              <  2:3897102/100‑1 (MQ=255)
                gATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGa                                                                              <  2:3978409/100‑1 (MQ=255)
                gATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGa                                                                              >  1:1472012/1‑100 (MQ=255)
                    gCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTg                                                                          >  1:3897406/1‑100 (MQ=255)
                     cTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGT‑AA‑AAATGAAACGGTGGGAAGTGGc                                                                        <  2:270619/100‑1 (MQ=255)
                      tttttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGc                                                                        <  1:7321/100‑1 (MQ=255)
                          tCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTa                                                                    >  1:1935426/1‑100 (MQ=255)
                            ccTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACt                                                                  >  1:2343153/1‑100 (MQ=255)
                            ccTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACt                                                                  >  1:1329650/1‑100 (MQ=255)
                             cTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTg                                                                 <  1:1801320/100‑1 (MQ=255)
                                cAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTgg                                                                >  2:4105196/1‑98 (MQ=255)
                                cAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTg                                                              >  2:454220/1‑100 (MQ=255)
                                cAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTg                                                              >  1:384618/1‑100 (MQ=255)
                                 aGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGa                                                             >  1:2897166/1‑100 (MQ=255)
                                    cATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATgg                                                          >  1:3993237/1‑100 (MQ=255)
                                    cATTTTTGAGTAATTCAATCATCGAACGGCAGTTCCCAACAACAATATAAAAACAAGA‑GGT‑AACAAATCAAACCTTCAGAAGTAGCTTTGCCGGttaacg                                                          >  2:3361527/1‑95 (MQ=255)
                                        tttGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAAt                                                      <  1:1460278/100‑1 (MQ=255)
                                           gTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGtt                                                   >  1:879474/1‑100 (MQ=255)
                                            tttAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTg                                                  >  2:2459229/1‑100 (MQ=255)
                                            tttAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGCGATGGTAATGTTg                                                  <  1:2785775/100‑1 (MQ=255)
                                             ttAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTgg                                                 >  2:880560/1‑100 (MQ=255)
                                             ttAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTgg                                                 <  1:592021/100‑1 (MQ=255)
                                              tAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGt                                                <  1:4133276/100‑1 (MQ=255)
                                               aaTTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGtt                                               <  1:1837264/100‑1 (MQ=255)
                                                 ttCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTg                                             <  1:2984480/100‑1 (MQ=255)
                                                 ttCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTg                                             <  1:840824/100‑1 (MQ=255)
                                                  tcaatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGc                                            <  1:2459243/100‑1 (MQ=255)
                                                  tcaatcaaCGAAAGGAAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGc                                            <  2:3337614/100‑1 (MQ=255)
                                                  tcaagcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGc                                            <  1:2017051/100‑1 (MQ=255)
                                                   caatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCt                                           >  1:1439307/1‑100 (MQ=255)
                                                    aatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTg                                          >  1:3363241/1‑100 (MQ=255)
                                                    aatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTg                                          >  2:1685333/1‑100 (MQ=255)
                                                    aatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTg                                          <  2:3582335/100‑1 (MQ=255)
                                                    aatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTg                                          <  1:3775992/100‑1 (MQ=255)
                                                    aatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTg                                          >  2:588430/1‑100 (MQ=255)
                                                     atcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAAAGGAGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGt                                         >  2:3838696/1‑100 (MQ=255)
                                                      tcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTa                                        >  2:1768067/1‑100 (MQ=255)
                                                       caaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTAt                                       >  2:2195035/1‑100 (MQ=255)
                                                       caaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTAt                                       >  2:2139641/1‑100 (MQ=255)
                                                       caaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTAt                                       <  1:92686/100‑1 (MQ=255)
                                                       caaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTAt                                       >  2:635379/1‑100 (MQ=255)
                                                        aaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTAtt                                      >  2:2229790/1‑100 (MQ=255)
                                                         aCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTg                                     <  1:3130730/100‑1 (MQ=255)
                                                           gAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGa                                    <  1:1232105/99‑1 (MQ=255)
                                                            aaaGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAAc                                  >  1:2705764/1‑100 (MQ=255)
                                                                gCAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACtgtg                              <  2:1189132/100‑1 (MQ=255)
                                                                 cAGTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGa                             >  2:276795/1‑100 (MQ=255)
                                                                   gTTCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAag                           <  1:276798/100‑1 (MQ=255)
                                                                    ttaccAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAag                           <  2:2468856/96‑1 (MQ=255)
                                                                    ttCCCAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAagg                          <  2:1732543/100‑1 (MQ=255)
                                                                      cccAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAaggag                        <  2:2428562/100‑1 (MQ=255)
                                                                      cccAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAaggag                        >  1:2820030/1‑100 (MQ=255)
                                                                      cccAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGAGGGAAGTGGCTTTACTGGTAATAGAAGAGTTGGATTGCTGTATTGAACTGTGAagg                          >  2:3460751/1‑98 (MQ=255)
                                                                       ccAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAaggagg                       <  1:2504773/100‑1 (MQ=255)
                                                                       ccAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAaggagg                       >  2:1546311/1‑100 (MQ=255)
                                                                       ccAACAACAAAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAaggagg                       >  2:3105122/1‑100 (MQ=255)
                                                                        caacaacaaAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTACTGTTGGGGTGGTGTATTGAAGGGTGCaggtgg                       >  2:3320442/1‑99 (MQ=255)
                                                                        caacaacaaAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAggagga                      >  2:4108043/1‑100 (MQ=255)
                                                                           caacaaAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGc                   >  1:1460174/1‑100 (MQ=255)
                                                                           caacaaAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGc                   <  1:3223046/100‑1 (MQ=255)
                                                                             acaaAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCa                 >  2:2835697/1‑100 (MQ=255)
                                                                             acaaAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCa                 >  2:1315002/1‑100 (MQ=255)
                                                                              caaAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCAt                <  2:3985119/100‑1 (MQ=255)
                                                                              caaAATAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCAt                >  1:1053578/1‑100 (MQ=255)
                                                                                aaaTAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGa              >  2:3971829/1‑100 (MQ=255)
                                                                                 aaTAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGaa             <  1:2197845/100‑1 (MQ=255)
                                                                                   tAAAGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAAc           <  1:738966/100‑1 (MQ=255)
                                                                                    aaaGATGGGG‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAAc           <  1:541337/99‑1 (MQ=255)
                                                                                     aaGATGGGA‑GCTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACtt         <  2:2705748/100‑1 (MQ=255)
                                                                                      aGATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACttt        >  1:452936/1‑100 (MQ=255)
                                                                                       gATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACtttt       >  2:151269/1‑100 (MQ=255)
                                                                                       gATGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACtttt       <  1:1184211/100‑1 (MQ=255)
                                                                                        aTGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTa      <  2:1988273/100‑1 (MQ=255)
                                                                                        aTGGGA‑GGT‑AA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTaa     <  1:710627/100‑1 (MQ=255)
                                                                                        aTGGGA‑GGT‑AA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTaa     <  1:3943726/100‑1 (MQ=255)
                                                                                         ttggA‑GGTAAA‑AAATGAAACGGGGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTAGCTGTATTGAACTGTGAAGGAGCACGCCATGAAACTTTTa      <  2:3055527/97‑1 (MQ=255)
                                                                                         tGGGA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTaa     <  2:2801089/100‑1 (MQ=255)
                                                                                           ggA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAATc   >  2:2804176/1‑100 (MQ=255)
                                                                                            gA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAATCg  <  2:3661325/100‑1 (MQ=255)
                                                                                            gA‑GGTAAA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAATCg  <  2:3878605/100‑1 (MQ=255)
                                                                                                  |                                                                                               
TGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGA‑GGT‑AA‑AAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAATCG  >  NC_000913/2563417‑2563607

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