Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,998,290 +G coding (379/486 nt) b4566 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,998,2891.G99.5% 715.9 / ‑3.5 207coding (380/486 nt)b4566hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (118/88);  total (119/88)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.34e-01

Only 100 of 226 total aligned reads displayed.

GACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCAT  >  NC_000913/3998197‑3998376
                                                                                             |                                                                                       
gACTGCGCCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCa                                                                                          <  1:5914669/93‑1 (MQ=255)
 aCTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGACCGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAg                                                                                        <  2:3446811/93‑1 (MQ=255)
 aCTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAg                                                                                        <  2:6345434/93‑1 (MQ=255)
  cTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAgg                                                                                       >  2:5022591/1‑92 (MQ=255)
  cTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAgg                                                                                       <  2:2380807/93‑2 (MQ=255)
    gCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAggac                                                                                     >  1:1960755/1‑90 (MQ=255)
       ccGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGa                                                                                   >  2:5544730/1‑93 (MQ=255)
       ccGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGa                                                                                   >  2:5350332/1‑93 (MQ=255)
        cGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGAt                                                                                  >  1:5678143/1‑93 (MQ=255)
        cGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGAt                                                                                  >  1:742253/1‑93 (MQ=255)
        cGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGAt                                                                                  >  2:3385915/1‑93 (MQ=255)
          cAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATaa                                                                                >  1:6341408/1‑93 (MQ=255)
            aaCTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAg                                                                              >  1:1019933/1‑93 (MQ=255)
            aaCTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAg                                                                              >  2:3472038/1‑93 (MQ=255)
            aaCTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAg                                                                              >  1:1198500/1‑93 (MQ=255)
              cTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGaa                                                                            >  1:5747215/1‑93 (MQ=255)
              cTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGaa                                                                            >  1:3135548/1‑93 (MQ=255)
              cTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGaa                                                                            >  2:5418133/1‑93 (MQ=255)
               tGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGaaa                                                                           >  1:4925159/1‑93 (MQ=255)
               tGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGaaa                                                                           >  2:1398793/1‑93 (MQ=255)
                 aaTCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAAtt                                                                         >  2:284285/1‑93 (MQ=255)
                 aaTCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAAtt                                                                         >  1:2189131/1‑93 (MQ=255)
                   tctcTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTcc                                                                       >  1:2278339/1‑93 (MQ=255)
                   tctcTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTcc                                                                       >  2:1955075/1‑93 (MQ=255)
                    ctctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTccc                                                                      <  2:3781805/93‑1 (MQ=255)
                    ctctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTccc                                                                      >  1:1157283/1‑93 (MQ=255)
                    ctctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTccc                                                                      >  1:2505795/1‑93 (MQ=255)
                     tctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACCTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCt                                                                     >  2:1975368/1‑93 (MQ=255)
                       tttACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTc                                                                   >  2:1286081/1‑93 (MQ=255)
                       tttACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTc                                                                   >  1:3601270/1‑93 (MQ=255)
                       tttACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTc                                                                   <  1:530546/93‑1 (MQ=255)
                        ttACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCa                                                                  <  1:2325337/93‑1 (MQ=255)
                            tGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTg                                                              <  2:4150579/93‑1 (MQ=255)
                            tGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTg                                                              >  2:825160/1‑93 (MQ=255)
                             ggTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGa                                                             <  1:4669444/93‑1 (MQ=255)
                             ggTGTTGAAGGCATCCATACGAACGATGCCGACTTTATCGCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGa                                                             <  1:3060390/93‑1 (MQ=255)
                              gTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGaa                                                            >  2:1017101/1‑93 (MQ=255)
                               ttttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAg                                                           <  2:6160382/93‑1 (MQ=255)
                                tttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGc                                                          <  1:53616/93‑1 (MQ=255)
                                 ttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCAACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCg                                                         >  2:5899333/1‑93 (MQ=255)
                                  tGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCga                                                        <  2:6418064/93‑1 (MQ=255)
                                   gAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCACTTCACCTGAAGCgag                                                       >  1:200513/1‑93 (MQ=255)
                                    aaGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCgaga                                                      >  1:5142601/1‑93 (MQ=255)
                                    aaGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCgaga                                                      <  2:3057534/93‑1 (MQ=255)
                                      ggCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGAtt                                                    >  2:3819759/1‑93 (MQ=255)
                                        cATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTa                                                  >  1:817140/1‑93 (MQ=255)
                                         aTCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTaa                                                 >  1:4020407/1‑93 (MQ=255)
                                         aTCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTaa                                                 <  2:5009267/93‑1 (MQ=255)
                                         aTCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTaa                                                 >  2:251031/1‑93 (MQ=255)
                                          tCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAAc                                                >  1:1634236/1‑93 (MQ=255)
                                             aGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGcc                                             >  2:2218777/1‑93 (MQ=255)
                                              gACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCa                                            >  2:4287622/1‑93 (MQ=255)
                                              gACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCa                                            >  2:6169608/1‑93 (MQ=255)
                                               acgaacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCAc                                           <  2:4118422/93‑1 (MQ=255)
                                               acgaacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCAc                                           >  1:895355/1‑93 (MQ=255)
                                               acgaacgaTGCCGACTTTACCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCAc                                           <  2:3730342/93‑1 (MQ=255)
                                                cgaacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCAcc                                          <  1:984724/93‑1 (MQ=255)
                                                 gaacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCt                                         >  1:5884950/1‑93 (MQ=255)
                                                 gaacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCt                                         <  2:4778685/93‑1 (MQ=255)
                                                 gaacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCt                                         <  2:4582895/93‑1 (MQ=255)
                                                 caccgaGGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCt                                         <  1:1394541/90‑1 (MQ=255)
                                                  aacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCtt                                        >  2:5199274/1‑93 (MQ=255)
                                                  aacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCtt                                        >  1:4758094/1‑93 (MQ=255)
                                                   acgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTc                                       >  1:5771825/1‑93 (MQ=255)
                                                   acgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTc                                       <  2:298539/93‑1 (MQ=255)
                                                     gaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctg                                     >  2:2097140/1‑93 (MQ=255)
                                                     gaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctg                                     >  1:4912267/1‑93 (MQ=255)
                                                      aTGCGGACTTTATCTCCCGCAACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctgc                                    <  1:4428613/93‑1 (MQ=255)
                                                      aTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctgc                                    >  1:1675652/1‑93 (MQ=255)
                                                      aTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctgc                                    <  2:5759152/93‑1 (MQ=255)
                                                           gACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGAcc                               <  2:734397/93‑1 (MQ=255)
                                                            aCTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCg                              >  1:3761892/1‑93 (MQ=255)
                                                              tttATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGtt                            <  1:4619821/93‑1 (MQ=255)
                                                              tttATCTCCCGCGACATAAAACACGCCGCCAGGCGGATAAAGAAATTCCCTTCAACTGAAGCGAGATTTAACGCCACCCTCTGCTGACCGttt                           >  2:6005468/1‑93 (MQ=255)
                                                               ttATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGttt                           <  1:5205993/93‑1 (MQ=255)
                                                                tATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGtttt                          <  1:3018697/93‑1 (MQ=255)
                                                                tATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGtttt                          <  2:3176299/93‑1 (MQ=255)
                                                                 aTCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTc                         >  1:6425678/1‑93 (MQ=255)
                                                                 aTCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTc                         <  1:1581494/93‑1 (MQ=255)
                                                                  tctcCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCg                        <  2:5694344/93‑1 (MQ=255)
                                                                  tctcCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCg                        <  2:2630739/93‑1 (MQ=255)
                                                                   ctcCCGCGACATAAAACACGCCGCCAGGACGGATAATGAAATTCACTTCGCCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGa                       <  2:6471779/93‑1 (MQ=255)
                                                                        gcgACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTAc                  <  1:2208004/93‑1 (MQ=255)
                                                                          gACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACaa                >  2:4498198/1‑93 (MQ=255)
                                                                          gACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACaa                >  1:3530907/1‑93 (MQ=255)
                                                                          gACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACaa                >  2:3078308/1‑93 (MQ=255)
                                                                           aCATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACaaa               >  1:1098111/1‑93 (MQ=255)
                                                                             aTAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAAcc             <  1:5446001/93‑1 (MQ=255)
                                                                              tAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCg            >  2:1757488/1‑93 (MQ=255)
                                                                              tAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCg            >  1:4294913/1‑93 (MQ=255)
                                                                                aaaCACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAg          <  2:1857985/93‑1 (MQ=255)
                                                                                  acacGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgcg        >  2:2800006/1‑93 (MQ=255)
                                                                                  acacGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgcg        >  1:3415268/1‑93 (MQ=255)
                                                                                    acGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgcggc      <  2:1641605/93‑1 (MQ=255)
                                                                                    acGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgcggc      <  2:1541971/93‑1 (MQ=255)
                                                                                    acGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgcggc      >  1:4083656/1‑93 (MQ=255)
                                                                                    acGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGCTTTCGATGTACAAACCGAgcggc      >  2:3647545/1‑93 (MQ=255)
                                                                                      gccgccAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGgcgc    <  1:3266985/93‑1 (MQ=255)
                                                                                       ccgccAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCa   >  1:2824798/1‑93 (MQ=255)
                                                                                        cgccAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCAt  <  2:2723683/93‑1 (MQ=255)
                                                                                             |                                                                                       
GACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCAT  >  NC_000913/3998197‑3998376

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