Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,071,321 A→G 100% intergenic (+110/‑218) yeeP → / → flu pseudogene, CP4‑44 prophage; 50S ribosome‑binding GTPase family;Phage or Prophage Related; putative histone/CP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,071,3210AG95.6% 330.6 / 0.7 104intergenic (+110/‑218)yeeP/flupseudogene, CP4‑44 prophage; 50S ribosome‑binding GTPase family;Phage or Prophage Related; putative histone/CP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand):  ref base A (1/4);  new base G (50/49);  total (51/53)
Fisher's exact test for biased strand distribution p-value = 3.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

GATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACA  >  NC_000913/2071225‑2071418
                                                                                                |                                                                                                 
gATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGata                                                                                                <  1:2142862/100‑1 (MQ=255)
gATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGata                                                                                                <  1:1990092/100‑1 (MQ=255)
  tGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATc                                                                                              >  1:2537601/1‑100 (MQ=255)
    ccctccctGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTa                                                                                            >  1:1159372/1‑100 (MQ=255)
     cctccctGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTAt                                                                                           <  1:1329133/100‑1 (MQ=255)
      ctccctGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTAtt                                                                                          <  1:2359794/100‑1 (MQ=255)
         cctGACTCTCAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTa                                                                                       <  1:2302466/100‑1 (MQ=255)
          ctGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTat                                                                                      >  1:801564/1‑100 (MQ=255)
           tGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTatc                                                                                     >  2:1081764/1‑100 (MQ=255)
             aCTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTatcga                                                                                   <  2:1404703/100‑1 (MQ=255)
              ctctGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATTGAATAAAACGATCGATATCTATTTTatcgat                                                                                  >  2:495808/1‑100 (MQ=255)
               tctGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTatcgatc                                                                                 <  1:490732/100‑1 (MQ=255)
                   aGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGttt                                                                             >  1:2551976/1‑100 (MQ=255)
                    gTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTta                                                                            <  2:1730638/100‑1 (MQ=255)
                    gTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCGATTTTatcgatcg                                                                                <  2:1586604/96‑1 (MQ=255)
                     tCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTtat                                                                           <  1:1061753/100‑1 (MQ=255)
                       tGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTatat                                                                         >  1:783279/1‑100 (MQ=255)
                       tGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTatat                                                                         >  2:702002/1‑100 (MQ=255)
                       tGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTatat                                                                         >  1:1391791/1‑100 (MQ=255)
                        gttCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatc                                                                        >  1:1733842/3‑100 (MQ=255)
                        gCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatc                                                                        >  2:1595180/1‑100 (MQ=255)
                        gCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatc                                                                        <  2:1733830/100‑1 (MQ=255)
                        gCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatc                                                                        >  1:1324154/1‑100 (MQ=255)
                        gCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatc                                                                        <  2:2333010/100‑1 (MQ=255)
                           cacaAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCCATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatccat                                                                     >  2:687802/1‑100 (MQ=255)
                           cacaAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATACCACGATCGATATCTATTTTATCGATCGTTTATatc                                                                        >  2:1200440/1‑97 (MQ=255)
                           cacaAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatcgat                                                                     >  1:353911/1‑100 (MQ=255)
                           cacaAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatcgat                                                                     >  1:767916/1‑100 (MQ=255)
                           cacaAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatcgat                                                                     >  2:80981/1‑100 (MQ=255)
                           cacaAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatcgat                                                                     <  1:17621/100‑1 (MQ=255)
                            acaacaGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATACTAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATatcgatc                                                                    <  2:865121/100‑1 (MQ=255)
                             caAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatcgatcg                                                                   >  1:886938/1‑100 (MQ=255)
                             caAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatcgatcg                                                                   >  2:1640050/1‑100 (MQ=255)
                              aaaaGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatcgatc                                                                    <  2:1212738/98‑1 (MQ=255)
                                aaGCTCTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGATATAAACGATCGATa                                                                <  1:2343043/100‑1 (MQ=38)
                                   cACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGc                                                             >  1:1441784/1‑100 (MQ=255)
                                   cACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGc                                                             <  2:717281/100‑1 (MQ=255)
                                   cACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGc                                                             >  2:1666140/1‑100 (MQ=255)
                                    aCTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCt                                                            >  2:1291560/1‑100 (MQ=255)
                                    aCTGTTTTCGTTACTGTCGCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGATATAAACGATCGATaaaat                                                            <  1:2326808/100‑4 (MQ=25)
                                     cTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCta                                                           >  1:556481/1‑100 (MQ=255)
                                      tGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCtaa                                                          >  2:17334/1‑100 (MQ=255)
                                         tttCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCtaataa                                                       <  2:585273/100‑1 (MQ=255)
                                         tttCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCtaataa                                                       >  1:1576736/1‑100 (MQ=255)
                                          ttCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCtaataat                                                      >  2:1195275/1‑100 (MQ=255)
                                             gTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCtaataata                                                     <  2:2141301/98‑1 (MQ=255)
                                             gTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAAc                                                   >  1:299787/1‑100 (MQ=255)
                                              ttACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAAcc                                                  <  1:2016596/100‑1 (MQ=255)
                                              ttACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAAcc                                                  <  1:1291568/100‑1 (MQ=255)
                                               tACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCt                                                 >  2:1225870/1‑100 (MQ=255)
                                                 cTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCttt                                               >  1:1889906/1‑100 (MQ=255)
                                                  tGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTg                                              <  1:2034485/100‑1 (MQ=255)
                                                  tGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTg                                              >  2:128095/1‑100 (MQ=255)
                                                  tGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTg                                              <  1:1527160/100‑1 (MQ=255)
                                                  tGTCTCTCTTCTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTCTATCGATCGATAAGCTAATAATAACCTTTg                                              <  2:1689168/100‑1 (MQ=255)
                                                  tGTCTCTCTGGTCCGAGCACTAGCTCCATAATCGACTAACACGATCGCTATCTATTTTATCGATCCTATATAGCGATCGATAAGCTAGTATTAACAttt                                               >  2:128452/1‑99 (MQ=255)
                                                    tctctcTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTc                                            <  1:934462/100‑1 (MQ=255)
                                                    tctctcTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTc                                            <  2:1979931/100‑1 (MQ=255)
                                                    tctctcTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTc                                            <  2:1815563/100‑1 (MQ=255)
                                                    tctctcTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGATCATATCGATCTATAACCTAATAATACCCTTTGTc                                            <  2:88034/100‑1 (MQ=255)
                                                     ctctctTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCa                                           >  1:1824176/1‑100 (MQ=255)
                                                        tctTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTa                                        <  2:1740087/100‑1 (MQ=255)
                                                           tGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACa                                     >  1:1681583/1‑100 (MQ=255)
                                                            gTCCGTGCAATAGCTCAATAATAGAATAAAACGATTGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACAt                                    <  1:1762780/100‑1 (MQ=255)
                                                              ccGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGc                                  <  2:2184131/100‑1 (MQ=255)
                                                              ccGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGc                                  >  2:507791/1‑100 (MQ=255)
                                                              ccGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGc                                  <  1:194769/100‑1 (MQ=255)
                                                               cGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGca                                 >  2:641445/1‑100 (MQ=255)
                                                                gTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGcac                                >  1:1200112/1‑100 (MQ=255)
                                                                 tGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGcaca                               <  1:2048289/100‑1 (MQ=255)
                                                                  gCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAg                              >  1:380773/1‑100 (MQ=255)
                                                                    aaTAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGAt                            >  2:2528621/1‑100 (MQ=255)
                                                                     aTAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATa                           >  1:624334/1‑100 (MQ=255)
                                                                     aTAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTGTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATa                           <  2:247625/100‑1 (MQ=255)
                                                                        gCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGt                        >  2:52026/1‑100 (MQ=255)
                                                                        gCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGt                        >  1:1467945/1‑100 (MQ=255)
                                                                          tCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTAc                      <  1:1007071/100‑1 (MQ=255)
                                                                          tCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTAc                      >  1:2064104/1‑100 (MQ=255)
                                                                           cAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACa                     <  1:2006853/100‑1 (MQ=255)
                                                                             ataataGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGa                   <  1:1990370/100‑1 (MQ=255)
                                                                              taataGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGaa                  >  2:2213322/1‑100 (MQ=255)
                                                                                ataGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAg                <  2:1174975/100‑1 (MQ=255)
                                                                                 taGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGa               <  2:2095440/100‑1 (MQ=255)
                                                                                  aGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGa               <  2:2016429/99‑1 (MQ=255)
                                                                                  aGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGAc              >  2:2363985/1‑100 (MQ=255)
                                                                                   gAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACa             >  1:195744/1‑100 (MQ=255)
                                                                                    aaTAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACAt            <  2:1475362/100‑1 (MQ=255)
                                                                                    aaTAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACAt            >  1:330053/1‑100 (MQ=255)
                                                                                    aaTAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACAt            <  2:1230608/100‑1 (MQ=255)
                                                                                    aaTAAAACGATCGATATCTATTTTATCGATCGTGTATATCGATCGATAAGCTAATAATAACCTTTGTCATTAACATGCACAGATATGTACAGAAAGAc              >  2:1614906/1‑98 (MQ=255)
                                                                                     aTAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAAACTTTGTCAGTAACATGCACAGATACGTACAGAAAGAGAtt           >  2:2059789/1‑100 (MQ=255)
                                                                                      tAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTc          <  2:649362/100‑1 (MQ=255)
                                                                                       aaaaCGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCa         <  1:2388606/100‑1 (MQ=255)
                                                                                        aaaCGATGGATATCTATTTTATCGATCGTGTATATCGATCGATAACCTAATTATTACCGTTGTCAGTTACATGCACAGATAAGTACAGAAAGACAtgc          >  2:1527962/1‑96 (MQ=255)
                                                                                        aaaCGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAg        <  2:298414/100‑1 (MQ=255)
                                                                                        aaaCGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAg        >  2:1725858/1‑100 (MQ=255)
                                                                                          aCGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAggg      <  1:1631939/100‑1 (MQ=255)
                                                                                           cgatcgatATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGa     >  1:688264/1‑100 (MQ=255)
                                                                                            gATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGaa    >  1:175122/1‑100 (MQ=255)
                                                                                              tCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTGACATGCACAGATACGTACAGAAAGACATTCAGGGaaca  <  2:1576725/100‑1 (MQ=255)
                                                                                                |                                                                                                 
GATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACA  >  NC_000913/2071225‑2071418

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: