Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,059,768 A→G intergenic (+44/‑58) yihM → / → yihN hypothetical protein/putative resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,059,7680AG85.1% 755.1 / 125.2 322intergenic (+44/‑58)yihM/yihNhypothetical protein/putative resistance protein
Reads supporting (aligned to +/- strand):  ref base A (20/28);  new base G (150/124);  total (170/152)
Fisher's exact test for biased strand distribution p-value = 1.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 324 total aligned reads displayed.

CGTACTACAACAAATCAAACAAGAGGCAAACCATCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAATGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCT  >  NC_000913/4059676‑4059860
                                                                                            |                                                                                            
cGTACTACAACAAATCAAACAAGAGGCAAACCATCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGaa                                                                                              >  1:5907848/1‑93 (MQ=255)
  tactacAACAAATCAAACAAGAGGCAAACCATCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTg                                                                                            >  1:6580932/1‑93 (MQ=255)
    ctacAACAAATCAAACAAGAGGCAAACCATCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGt                                                                                          <  2:6874566/93‑1 (MQ=255)
     tacAACAAATCAAACAAGAGGCAAACCATCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTc                                                                                         >  1:5025057/1‑93 (MQ=255)
      acaacaAATCAAACAAGAGGCAAACCATCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCg                                                                                        <  2:3360713/93‑1 (MQ=255)
      acaacaAATCAAACAAGAGGCAAACCATCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCg                                                                                        >  1:2980671/1‑93 (MQ=255)
         acaaATCAAACAAGAGGCAAACCATCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGGGTGGTCGGTa                                                                                     >  2:7288810/1‑93 (MQ=255)
            aaTCAAACAAGAGGCAAACCATCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCt                                                                                  <  2:1395585/93‑1 (MQ=255)
            aaTCAAACAAGAGGCAAACCATCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCt                                                                                  <  2:3560062/93‑1 (MQ=255)
              tCAAACAAGAGGCAAACCATCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTa                                                                                <  2:6804126/93‑1 (MQ=255)
              tCAAACAAGAGGCAAACCATCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTa                                                                                >  2:1208565/1‑93 (MQ=255)
              tCAAACAAGAGGCAAACCATCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTa                                                                                <  2:7207554/93‑1 (MQ=255)
               cAAACAAGAGGCAAACCATCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTAt                                                                               <  2:3922413/93‑1 (MQ=255)
                aaaCAAGAGGCAAACCATCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTAtt                                                                              <  2:7419374/93‑1 (MQ=255)
                 aaCAAGAGGCAAACCATCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTg                                                                             <  2:472000/93‑1 (MQ=255)
                   cAAGAGGCAAACCATCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGcc                                                                           >  1:539985/1‑93 (MQ=255)
                   cAAGAGGCAAACCATCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGcc                                                                           >  1:7117223/1‑93 (MQ=255)
                    aaGAGGCAAACCATCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCt                                                                          >  1:1019428/1‑93 (MQ=255)
                    aaGAGGCAAACCATCTTCTGAACCACTAATAAGCTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCt                                                                          <  2:6525864/93‑1 (MQ=255)
                     agagGCAAACCATCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCtt                                                                         <  1:854950/93‑1 (MQ=255)
                     agagGCAAACCATCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCtt                                                                         >  2:1804904/1‑93 (MQ=255)
                      gagGCAAACCATCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCttt                                                                        >  1:2506009/1‑93 (MQ=255)
                           aaaCCATCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAAc                                                                   <  2:4079627/93‑1 (MQ=255)
                            aaCCATCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACa                                                                  >  1:1787151/1‑93 (MQ=255)
                            aaCCATCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACa                                                                  <  2:700235/93‑1 (MQ=255)
                             aCCATCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAg                                                                 <  1:3451791/93‑1 (MQ=255)
                              ccATCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGt                                                                <  1:3531002/93‑1 (MQ=255)
                               cATCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTg                                                               >  2:3669793/1‑93 (MQ=255)
                                aTCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGa                                                              <  1:6906764/93‑1 (MQ=255)
                                aTCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGa                                                              >  2:5168703/1‑93 (MQ=255)
                                aTCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTGAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGa                                                              >  1:6216609/1‑93 (MQ=255)
                                 tcttctGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGAt                                                             >  2:6708776/1‑93 (MQ=255)
                                 tcttctGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGAt                                                             >  1:305285/1‑93 (MQ=255)
                                  cttctGAACAACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAATGGTCGGTATCTTATTGCCTTTAAAACAGTGAtt                                                            <  2:2484335/93‑1 (MQ=255)
                                    tctGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAATGGTCGGTATCTTATTGCCTTTAAAACAGTGAtttt                                                          <  1:5783920/93‑1 (MQ=255)
                                      tGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAATGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTaa                                                        <  1:1874464/93‑1 (MQ=255)
                                         aCCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTaaaaa                                                     <  1:1238215/93‑1 (MQ=255)
                                         aCCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGCATCTTATTGCCTTTAAAACAGTGATTTTaaaaa                                                     >  1:7370266/1‑93 (MQ=255)
                                          ccACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAAt                                                    <  1:5062816/93‑1 (MQ=255)
                                          ccACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAAt                                                    <  2:4821499/93‑1 (MQ=255)
                                           cACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATc                                                   >  1:1135602/1‑93 (MQ=255)
                                            aCTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCg                                                  >  2:5667344/1‑93 (MQ=255)
                                            aCTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAATGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCg                                                  >  1:1360514/1‑93 (MQ=255)
                                             cTAATAAGTTTTAACTCATACGCAGTTCATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGc                                                 >  2:6420757/1‑93 (MQ=255)
                                             cTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCNc                                                 >  2:6834901/1‑93 (MQ=255)
                                             cTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGc                                                 >  2:15057/1‑93 (MQ=255)
                                             cTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGc                                                 >  2:752070/1‑93 (MQ=255)
                                              taataaGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCt                                                >  2:2584434/1‑93 (MQ=255)
                                              taataaGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCt                                                >  2:5886174/1‑93 (MQ=255)
                                               aataaGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCtt                                               >  2:7069447/1‑93 (MQ=255)
                                               aataaGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCtt                                               >  1:5837450/1‑93 (MQ=255)
                                                ataaGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTa                                              <  2:6905107/93‑1 (MQ=255)
                                                ataaGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTa                                              >  2:108728/1‑93 (MQ=255)
                                                 taaGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAg                                             <  2:3249017/93‑1 (MQ=255)
                                                 taaGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAg                                             >  2:5532489/1‑93 (MQ=255)
                                                 taaGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAg                                             >  1:5311470/1‑93 (MQ=255)
                                                     ttttAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGa                                         <  1:1103295/93‑1 (MQ=255)
                                                     ttttAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGa                                         <  1:1977368/93‑1 (MQ=255)
                                                      tttAACTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGAt                                        >  2:2776529/1‑93 (MQ=255)
                                                          aCTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTa                                    <  2:1335627/93‑1 (MQ=255)
                                                           cTCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTAt                                   <  2:4571835/93‑1 (MQ=255)
                                                            tCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATg                                  >  1:4802804/1‑93 (MQ=255)
                                                            tCATACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATg                                  >  2:2199574/1‑93 (MQ=255)
                                                              aTACGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCt                                >  2:4186877/1‑93 (MQ=255)
                                                              aTACGCAGTTAATAAAACCATTCCATATGAATGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCt                                <  2:4456672/93‑1 (MQ=255)
                                                                aCGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCa                              <  1:6328103/93‑1 (MQ=255)
                                                                aCGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCa                              <  1:2310826/93‑1 (MQ=255)
                                                                aCGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCa                              <  1:4979742/93‑1 (MQ=255)
                                                                aCGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCa                              <  2:6151107/93‑1 (MQ=255)
                                                                 cGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCAc                             >  1:6412636/1‑93 (MQ=255)
                                                                 cGCAGTTAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCAc                             <  1:939745/93‑1 (MQ=255)
                                                                   cAGTTAATAAAACCATTCCATATGAATGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCACGa                           >  1:4590003/1‑93 (MQ=255)
                                                                     gTTCATAAAACCATTCCATATGAATGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCACGaaa                         >  1:6977573/1‑93 (MQ=255)
                                                                      ttAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCACGaaaa                        >  1:4844832/1‑93 (MQ=255)
                                                                      ttAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCACGaaaa                        >  2:6033123/1‑93 (MQ=255)
                                                                      ttAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCACGaaaa                        >  1:6749781/1‑93 (MQ=255)
                                                                      ttAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCACGaaaa                        <  2:5691731/93‑1 (MQ=255)
                                                                      ttAATAAAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCACGaaaa                        <  1:5700734/93‑1 (MQ=255)
                                                                      ttAATAAAACCATTCCATATGAATGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCACGaaaa                        >  2:4916058/1‑93 (MQ=255)
                                                                      ttAATAAAACCATTCCATATGAATGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCACGaaaa                        <  1:4385410/93‑1 (MQ=255)
                                                                       taataaAACCATTCCATATGAATGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCACGaaaaa                       <  1:3148434/93‑1 (MQ=255)
                                                                         ataaAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCACGAAAAAGa                     <  2:3179066/93‑1 (MQ=255)
                                                                         ataaAACCATTCCATATGAATGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCACGAAAAAGa                     <  2:3002802/93‑1 (MQ=255)
                                                                          taaAACCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCACGAAAAAGaa                    >  1:5801024/1‑93 (MQ=255)
                                                                           aaaaCCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCACGAAAAAGaaa                   <  2:1563009/93‑1 (MQ=255)
                                                                           aaaaCCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCACGAAAAAGaaa                   <  1:5137049/93‑1 (MQ=255)
                                                                           aaaaCCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCACGAAAAAGaaa                   <  1:1647735/93‑1 (MQ=255)
                                                                            aaaCCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCACGAAAAAGAAAt                  >  1:7096263/1‑93 (MQ=255)
                                                                            aaaCCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCACGAAAAAGAAAt                  >  2:1912319/1‑93 (MQ=255)
                                                                             aaCCATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATg                 >  1:3813442/1‑93 (MQ=255)
                                                                             aaCCATTCCATATGAATGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATg                 <  1:3928106/93‑1 (MQ=255)
                                                                               ccATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATggg               <  2:1046447/93‑1 (MQ=255)
                                                                               ccATTCCATATGAATGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATggg               <  2:5323935/93‑1 (MQ=255)
                                                                               ccATTCCATATGAATGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATggg               <  1:3737700/93‑1 (MQ=255)
                                                                                cATTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGc              <  1:2644963/93‑1 (MQ=255)
                                                                                 aTTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCg             >  1:5530200/1‑93 (MQ=255)
                                                                                 aTTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCg             >  1:2191096/1‑93 (MQ=255)
                                                                                 aTTCCATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCg             >  1:1733681/1‑93 (MQ=255)
                                                                                    ccATATGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTa          <  1:4255377/93‑1 (MQ=255)
                                                                                        atGAGTGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTa      <  1:7192361/93‑1 (MQ=255)
                                                                                            |                                                                                            
CGTACTACAACAAATCAAACAAGAGGCAAACCATCTTCTGAACCACTAATAAGTTTTAACTCATACGCAGTTAATAAAACCATTCCATATGAATGGTCGGTATCTTATTGCCTTTAAAACAGTGATTTTAAAAATCGCTTAGGGGATGTTATGCTCACGAAAAAGAAATGGGCGTTATTTAGTCT  >  NC_000913/4059676‑4059860

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