Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,472,408 A→G intergenic (+117/‑28) b4251 → / → b4252 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,472,4080AG99.9% 2269.4 / NA 690intergenic (+117/‑28)b4251/b4252hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (338/351);  total (339/351)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 702 total aligned reads displayed.

AAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGAT  >  NC_000913/4472318‑4472497
                                                                                          |                                                                                         
aaGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCg                                                                                         <  1:4232550/93‑1 (MQ=255)
aaGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCg                                                                                         <  2:1329868/93‑1 (MQ=255)
  gcgcCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCgtt                                                                                       >  1:9237664/1‑93 (MQ=255)
    cccTGCCGATAAGTTGTGATGAGACAATAATTCACGCGTTAAGGCTAGCGGAATGGATTATCTTTTCGTATAACGATAGAAATGAAGCgttgt                                                                                     <  2:7770706/92‑1 (MQ=255)
     ccTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCgttgtt                                                                                    <  1:3988927/93‑1 (MQ=255)
      cTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGttt                                                                                   >  1:9711251/1‑93 (MQ=255)
       tGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGtttt                                                                                  <  2:3888178/93‑1 (MQ=255)
       tGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGtttt                                                                                  >  1:2479231/1‑93 (MQ=255)
       tGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGtttt                                                                                  >  2:9133873/1‑93 (MQ=255)
       tGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGtttt                                                                                  >  2:8551142/1‑93 (MQ=255)
        gCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTtta                                                                                 >  1:5935049/1‑93 (MQ=255)
          cGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTttaat                                                                               >  1:7063819/1‑93 (MQ=255)
            aTAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTttaatta                                                                             <  2:411211/93‑1 (MQ=255)
            aTAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTttaatta                                                                             >  2:7436716/1‑93 (MQ=255)
             tAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTttaattaa                                                                            <  1:1261060/93‑1 (MQ=255)
               aGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAgg                                                                          <  2:4979139/93‑1 (MQ=255)
               aGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAgg                                                                          >  1:9647283/1‑93 (MQ=255)
               aGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAgg                                                                          >  2:3062314/1‑93 (MQ=255)
                gTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGa                                                                         <  1:3426621/93‑1 (MQ=255)
                gTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGa                                                                         <  1:9069708/93‑1 (MQ=255)
                  tgtgATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGt                                                                       <  1:5227338/93‑1 (MQ=255)
                     gATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGa                                                                    >  2:3637778/1‑93 (MQ=255)
                       tAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGaaa                                                                  >  2:1708306/1‑93 (MQ=255)
                       tAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGaaa                                                                  >  1:2052162/1‑93 (MQ=255)
                         aGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAg                                                                >  1:8864475/1‑93 (MQ=255)
                         aGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAg                                                                >  2:6588639/1‑93 (MQ=255)
                            cAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAgatg                                                             <  2:6774567/93‑1 (MQ=255)
                             aataatTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAgatga                                                            >  1:9922492/1‑93 (MQ=255)
                              ataatTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAgatgat                                                           <  2:8729206/93‑1 (MQ=255)
                              ataatTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAgatgat                                                           <  1:6820296/93‑1 (MQ=255)
                               taatTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGatc                                                          >  1:9849858/1‑93 (MQ=255)
                               taatTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGatc                                                          <  1:8126671/93‑1 (MQ=255)
                                aatTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGatca                                                         <  1:4692210/93‑1 (MQ=255)
                                 atTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGatcat                                                        <  2:2692980/93‑1 (MQ=255)
                                 atTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGatcat                                                        <  1:8471636/93‑1 (MQ=255)
                                 atTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGatcat                                                        <  2:2121318/93‑1 (MQ=255)
                                  ttCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGatcatc                                                       <  2:7213281/93‑1 (MQ=255)
                                  ttCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGatcatc                                                       >  2:6886063/1‑93 (MQ=255)
                                   tCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCg                                                      <  1:6104075/93‑1 (MQ=255)
                                    cACGCGTTAATGCTAGCGGAATTGATTATCTGTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCgg                                                     <  1:1751142/93‑1 (MQ=255)
                                     aCGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGa                                                    <  1:3280359/93‑1 (MQ=255)
                                       gcgTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGaaa                                                  <  2:1178727/93‑1 (MQ=255)
                                       gcgTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGaaa                                                  <  1:6368361/93‑1 (MQ=255)
                                       gcgTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGaaa                                                  <  2:3251022/93‑1 (MQ=255)
                                        cgTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGaact                                                 >  1:8636473/1‑91 (MQ=255)
                                        cgTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAat                                                 >  2:1951896/1‑93 (MQ=255)
                                          ttAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAatat                                               >  2:2819013/1‑93 (MQ=255)
                                          ttAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAatat                                               >  1:7609208/1‑93 (MQ=255)
                                          ttAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAatat                                               <  2:3785495/93‑1 (MQ=255)
                                          ttAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAatat                                               >  2:4444155/1‑93 (MQ=255)
                                          ttAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAatat                                               <  2:9451975/93‑1 (MQ=255)
                                            aaGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTc                                             <  1:10005519/93‑1 (MQ=255)
                                            aaGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTc                                             >  1:2917602/1‑93 (MQ=255)
                                             aGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCa                                            >  2:5765879/1‑93 (MQ=255)
                                             aGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCa                                            <  1:435382/93‑1 (MQ=255)
                                             aGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCa                                            >  1:2721128/1‑93 (MQ=255)
                                               gCTAGCGGAATTGATTATCTTTTCGGATAACGATAGAAATGAAGCGTTGTTTTCATTAAGTAGTGGAAAAGATGATCATCGGAAATATTCATa                                          <  2:4297129/93‑1 (MQ=255)
                                                cTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATaa                                         <  1:9357964/93‑1 (MQ=255)
                                                cTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATaa                                         <  2:5138174/93‑1 (MQ=255)
                                                 tAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAAt                                        <  1:7430717/93‑1 (MQ=255)
                                                 tAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAAt                                        >  1:5577488/1‑93 (MQ=255)
                                                 tAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAAt                                        <  1:5366378/93‑1 (MQ=255)
                                                 tAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAAt                                        >  1:2300842/1‑93 (MQ=255)
                                                   gCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCt                                      >  2:8391309/1‑93 (MQ=255)
                                                    cGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCtt                                     <  1:1105282/93‑1 (MQ=255)
                                                    cGGAATTGATTATCTTGTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCtt                                     >  2:9557343/1‑93 (MQ=255)
                                                     ggAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTc                                    >  1:4117891/1‑93 (MQ=255)
                                                     ggAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTc                                    <  2:5370700/93‑1 (MQ=255)
                                                      gAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTCAGGAGTGGAAAAGACGATCATCGGAAATATTCATAATCTTCa                                   >  1:4928174/1‑93 (MQ=255)
                                                      gAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCa                                   <  1:3651134/93‑1 (MQ=255)
                                                      gAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCa                                   >  1:3851650/1‑93 (MQ=255)
                                                        aTTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGc                                 >  1:8277878/1‑93 (MQ=255)
                                                        aTTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGc                                 <  2:4711531/93‑1 (MQ=255)
                                                        aTTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGc                                 <  2:5600897/93‑1 (MQ=255)
                                                          tGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCg                               >  1:7929799/1‑93 (MQ=255)
                                                           gATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGt                              >  1:2623978/1‑93 (MQ=255)
                                                           gATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGt                              >  1:6919838/1‑93 (MQ=255)
                                                           gATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGt                              <  2:2397442/93‑1 (MQ=255)
                                                             ttATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTgg                            <  1:6524151/93‑1 (MQ=255)
                                                              tATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGc                           <  2:64610/93‑1 (MQ=255)
                                                              tATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGc                           <  2:7684519/93‑1 (MQ=255)
                                                              tATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGc                           <  2:5249902/93‑1 (MQ=255)
                                                              tATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGc                           <  2:1013108/93‑1 (MQ=255)
                                                               aTCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCt                          <  2:4541617/93‑1 (MQ=255)
                                                                  ttttCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTAcc                       <  1:1641742/93‑1 (MQ=255)
                                                                  ttttCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTAcc                       >  2:1542294/1‑93 (MQ=255)
                                                                  ttttCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTAcc                       >  1:3106323/1‑93 (MQ=255)
                                                                   tttCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTAccc                      >  1:1401173/1‑93 (MQ=255)
                                                                      cGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAg                   <  2:6051752/93‑1 (MQ=255)
                                                                       gTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAgg                  >  1:6638711/1‑93 (MQ=255)
                                                                        tataACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGa                 <  1:7191845/93‑1 (MQ=255)
                                                                              gATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACg           >  2:3435975/1‑93 (MQ=255)
                                                                                aaGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGcc         <  2:7588681/92‑1 (MQ=255)
                                                                                 aGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCa        <  1:5687739/93‑1 (MQ=255)
                                                                                  gAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAg       <  1:7801979/93‑1 (MQ=255)
                                                                                   aaaTGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAgg      <  1:4827817/93‑1 (MQ=255)
                                                                                    aaTGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGc     <  2:89217/93‑1 (MQ=255)
                                                                                    aaTGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGc     <  1:5057976/93‑1 (MQ=255)
                                                                                      tGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGa   >  2:8062359/1‑93 (MQ=255)
                                                                                       gAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGAt  <  2:8252722/93‑1 (MQ=255)
                                                                                          |                                                                                         
AAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGAT  >  NC_000913/4472318‑4472497

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