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,3210AG98.7% 514.8 / ‑3.1 146intergenic (+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/1);  new base G (78/66);  total (79/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

Only 100 of 149 total aligned reads displayed.

GATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAA  >  NC_000913/2071225‑2071416
                                                                                                |                                                                                                
gATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGata                                                                                               <  1:1293058/100‑1 (MQ=255)
  tGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATc                                                                                             <  2:2996665/100‑1 (MQ=255)
   gCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCt                                                                                            >  1:1031453/1‑100 (MQ=255)
    ccctccctGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTa                                                                                           >  2:2547119/1‑100 (MQ=255)
    ccctccctGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTa                                                                                           >  1:582849/1‑100 (MQ=255)
     cctccctGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTAt                                                                                          >  2:2152222/1‑100 (MQ=255)
       tccctGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTAttt                                                                                        <  2:2500153/100‑1 (MQ=255)
        ccctGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTAtttt                                                                                       >  1:255501/1‑100 (MQ=255)
           tGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTatc                                                                                    <  1:2662844/100‑1 (MQ=255)
           tGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTatc                                                                                    <  2:2748730/100‑1 (MQ=255)
           tGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTatc                                                                                    <  2:438041/100‑1 (MQ=255)
           tGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTatc                                                                                    <  1:1156208/100‑1 (MQ=255)
           tGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTatc                                                                                    <  2:2855791/100‑1 (MQ=255)
            gACTCTGAGTCTGGTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTatcg                                                                                   >  2:2277409/1‑100 (MQ=255)
              ctctGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTTTTGTCCGTGCAATAGCTCATTAATAGAATAAAACGATCGATATCTATTTTatcgat                                                                                 >  2:2377441/1‑100 (MQ=255)
                   aGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGttt                                                                            <  2:2826279/100‑1 (MQ=255)
                    gTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTta                                                                           <  1:3201832/100‑1 (MQ=255)
                     tccgatctcaAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTtat                                                                          <  1:1516011/92‑1 (MQ=255)
                     tCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTtat                                                                          <  1:2731989/100‑1 (MQ=255)
                      cTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATTATAGAATAAAACTATCGATATCTATTTTATCGATCGttt                                                                            >  1:110636/1‑97 (MQ=255)
                      cTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTtata                                                                        >  2:2894050/1‑100 (MQ=255)
                      cTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTtata                                                                        >  1:914682/1‑100 (MQ=255)
                       tGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTatat                                                                       >  1:1210891/1‑100 (MQ=255)
                       tGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTatat                                                                       <  2:38230/100‑1 (MQ=255)
                       gccTCACAAAAGCACTGTTTTCATTACTGTCCCTCTAGTCCGTGCATTAGCTTAATAATAGAATAAAACGATCGATACCTATTTTatcgatccgtta                                                                           <  1:2024720/95‑6 (MQ=255)
                        gCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatc                                                                      >  1:1671681/1‑100 (MQ=255)
                         cTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatcg                                                                     >  1:1061872/1‑100 (MQ=255)
                         cTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatcg                                                                     >  2:3336077/1‑100 (MQ=255)
                         cTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatcg                                                                     <  1:2431729/100‑1 (MQ=255)
                           cacaAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatcgat                                                                   >  1:695013/1‑100 (MQ=255)
                           cacaAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatcgat                                                                   >  1:1428402/1‑100 (MQ=255)
                            acaAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatcgatc                                                                  <  2:1794594/100‑1 (MQ=255)
                            acaAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATatcgatc                                                                  <  1:2982580/100‑1 (MQ=255)
                             caAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATTGATCGTTTATatcgatcg                                                                 >  1:2541265/1‑100 (MQ=255)
                             caAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTtatagatcgga                                                                 >  2:1516018/1‑93 (MQ=255)
                               aaaGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATcgatcgat                                                               <  2:1123564/100‑1 (MQ=255)
                                 aGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTAAATAATAAAATAAAACGATCGATATCTATTTTATCGATCGTTTATatcaatc                                                                  >  1:630245/1‑95 (MQ=255)
                                  gCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAg                                                            >  2:2909459/1‑100 (MQ=255)
                                    aCTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCt                                                          <  1:3206269/100‑1 (MQ=255)
                                     cTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCta                                                         >  2:3364664/1‑100 (MQ=255)
                                      tGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCtaa                                                        <  1:1521381/100‑1 (MQ=255)
                                      tGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGATATAAACGATCGATaaaatag                                                        <  2:1471963/100‑6 (MQ=25)
                                         tttCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCtaataa                                                     <  1:2985092/100‑1 (MQ=255)
                                         tttCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCtaataa                                                     >  2:2533379/1‑100 (MQ=255)
                                           tCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCtaataata                                                   <  1:1648290/100‑1 (MQ=255)
                                           tCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCtaataata                                                   >  1:580734/1‑100 (MQ=255)
                                           tCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCtaataata                                                   >  1:2679868/1‑100 (MQ=255)
                                             gTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCtaataataa                                                  >  2:3057389/1‑99 (MQ=255)
                                             gTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAAc                                                 >  1:3357151/1‑100 (MQ=255)
                                              ttACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAAcc                                                <  1:2268176/100‑1 (MQ=255)
                                              ccgatctctctcTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAAcc                                                <  1:2957592/94‑1 (MQ=255)
                                                aCTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCtt                                              >  2:2122534/1‑100 (MQ=255)
                                                    tctctcTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTc                                          <  1:2932521/100‑1 (MQ=255)
                                                    tctctcTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTc                                          <  2:2993650/100‑1 (MQ=255)
                                                    tctctcGTGTCCGTGCTATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTc                                          <  1:2833550/100‑1 (MQ=255)
                                                     ctctctTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCa                                         >  2:1468756/1‑100 (MQ=255)
                                                     ctctctTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCa                                         >  1:2437945/1‑100 (MQ=255)
                                                     ctctctTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCa                                         >  2:1880743/1‑100 (MQ=255)
                                                      tctctTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAg                                        <  1:495447/100‑1 (MQ=255)
                                                      tctctTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAg                                        <  2:263178/100‑1 (MQ=255)
                                                      tctctTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAg                                        <  1:2911147/100‑1 (MQ=255)
                                                       ctctTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGt                                       >  2:1269048/1‑100 (MQ=255)
                                                        tctTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTa                                      >  2:2915050/1‑100 (MQ=255)
                                                        tctTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTa                                      <  2:3080867/100‑1 (MQ=255)
                                                        tctTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTa                                      <  2:1753256/100‑1 (MQ=255)
                                                           tGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATGAGCTAATAATAACCTTTGTCAGTACCa                                   <  1:1918166/100‑1 (MQ=255)
                                                           tGTCCGTGCAATAGCTAAATAATAGAATAAAACGATTGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACa                                   <  1:373529/100‑1 (MQ=255)
                                                            gTCCGTGCAATAGCTCAATAGTAGAATAAAACGATCGATATCTATTTTATCGATCGTTGATATCGATCGATAAGCTAATAATAACCTTTGTCAGt                                       <  1:1006892/95‑1 (MQ=255)
                                                             tCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATg                                 <  1:1746265/100‑1 (MQ=255)
                                                             tCCGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTATCAGTAACATg                                 >  2:759348/1‑100 (MQ=255)
                                                              ccGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGc                                <  2:772997/100‑1 (MQ=255)
                                                              ccGTGCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGc                                <  1:1438259/100‑1 (MQ=255)
                                                               cGTGCAATAGCTCAATACTAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAAACTTTGTCAGTAACATg                                 >  1:2852672/1‑98 (MQ=255)
                                                               cGTGCAATAGCTCAATAATAGAATAAAACGATTGATATCTATTTTATCGATCGTTTATATCGATCCATAAGCTAATAATAACCTTCGTCAGTAACATGc                                <  1:257363/99‑1 (MQ=255)
                                                                  gCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAg                            >  2:1982803/1‑100 (MQ=255)
                                                                  gCAATAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAg                            >  1:2126286/1‑100 (MQ=255)
                                                                      tAGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATAc                        >  2:1726745/1‑100 (MQ=255)
                                                                       aGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACg                       <  2:582852/100‑1 (MQ=255)
                                                                       aGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACg                       >  2:132904/1‑100 (MQ=255)
                                                                       aGCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACg                       <  1:740242/100‑1 (MQ=255)
                                                                        gCTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATATGCTAATAATAACCTTTGTCAGTAACATGCACAGCTACGt                      >  2:1997540/1‑100 (MQ=255)
                                                                         cTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAAGTTTTGTCAGTAACATGCACAGATAc                        >  1:1994206/1‑97 (MQ=255)
                                                                         cTCAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTa                     >  1:2712376/1‑100 (MQ=255)
                                                                           cAATAATAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACa                   >  1:2180492/1‑100 (MQ=255)
                                                                            aataatAGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACa                   <  1:50012/99‑1 (MQ=255)
                                                                                ataGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAg              >  2:1489731/1‑100 (MQ=255)
                                                                                 taGAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGa             >  2:1904017/1‑100 (MQ=255)
                                                                                   gAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACa           >  2:1144083/1‑100 (MQ=255)
                                                                                   gAATAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACa           >  2:1806768/1‑100 (MQ=255)
                                                                                    aaTAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACAt          >  2:2627608/1‑100 (MQ=255)
                                                                                    aaTAAAACGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACAt          <  2:466424/100‑1 (MQ=255)
                                                                                        aaaCGATCGATATCTATTTTATCGATCGTTTATAACGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCa       >  1:922991/1‑99 (MQ=255)
                                                                                         aaCGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAgg     <  1:673055/100‑1 (MQ=255)
                                                                                          aCGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAggg    >  1:2949077/1‑100 (MQ=255)
                                                                                          aCGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAggg    <  1:2531733/100‑1 (MQ=255)
                                                                                          aCGATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAggg    <  2:2109808/100‑1 (MQ=255)
                                                                                          aCGATCGATATCTATTTTATCGATCGTTGATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAggg    <  1:760221/100‑1 (MQ=255)
                                                                                           cgatcgatATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAggg    >  1:2260383/1‑99 (MQ=255)
                                                                                           cgatcgatATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGa   <  1:1904004/100‑1 (MQ=255)
                                                                                            gATCGATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGaa  >  1:1734964/1‑100 (MQ=255)
                                                                                                |                                                                                                
GATGCCCTCCCTGACTCTGAGTCTGCTCACAAAAGCACTGTTTTCGTTACTGTCTCTCTTGTCCGTGCAATAGCTCAATAATAGAATAAAACGATCAATATCTATTTTATCGATCGTTTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAA  >  NC_000913/2071225‑2071416

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: