Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 34,120 (T)9→8 coding (146/180 nt) yaaV ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091334,1120T.99.6% 1011.9 / ‑3.8 220coding (154/180 nt)yaaVhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (96/123);  minor base C (1/0);  total (97/123)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 250 total aligned reads displayed.

GCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAG  >  NC_000913/34020‑34200
                                                                                            |                                                                                        
gcacgcacAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGt                                                                                          <  1:3642108/93‑1 (MQ=255)
gcacgcacAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGt                                                                                          >  2:2187854/1‑93 (MQ=255)
 cacgcacAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGtt                                                                                         >  1:1534259/1‑93 (MQ=255)
  acgcacAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGttt                                                                                        >  2:3292870/1‑93 (MQ=255)
      acaGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGttttttt                                                                                    >  2:5877307/1‑93 (MQ=255)
        aGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGttttttttg                                                                                  <  1:936813/93‑2 (MQ=255)
         gATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGttttttttgc                                                                                 <  2:4891677/93‑3 (MQ=255)
          aTCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGttttttttgcc                                                                                >  1:5321666/1‑90 (MQ=255)
           tCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGttttttttgccc                                                                               >  1:1285971/1‑89 (MQ=255)
             aaaTAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCAt                                                                            <  1:3852127/93‑1 (MQ=255)
              aataatAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCAtt                                                                           <  1:2476545/93‑1 (MQ=255)
                taataGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTaa                                                                         >  1:1717306/1‑93 (MQ=255)
                taataGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTaa                                                                         >  1:958502/1‑93 (MQ=255)
                taataGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTaa                                                                         <  2:5496010/93‑1 (MQ=255)
                 aataGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAg                                                                        >  1:2172721/1‑93 (MQ=255)
                  ataGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGt                                                                       <  2:46788/93‑1 (MQ=255)
                   taGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTa                                                                      <  1:2328342/93‑1 (MQ=255)
                   taGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTa                                                                      >  2:2009223/1‑93 (MQ=255)
                    aGCATGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTaa                                                                     <  2:5446583/93‑1 (MQ=255)
                     gCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTaaa                                                                    <  2:4739077/93‑1 (MQ=255)
                       gtgtCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATc                                                                  >  2:3701542/1‑93 (MQ=255)
                          tCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCttt                                                               >  1:5172819/1‑93 (MQ=255)
                          tCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCttt                                                               <  2:6234489/93‑1 (MQ=255)
                           cATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCtttt                                                              <  2:3926609/93‑1 (MQ=255)
                            aTGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGTCCATTAAGTAAATCTTTTg                                                             >  2:2297119/1‑93 (MQ=255)
                            aTGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTg                                                             >  2:2127084/1‑93 (MQ=255)
                            aTGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTg                                                             >  2:498081/1‑93 (MQ=255)
                            aTGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTg                                                             >  2:4218803/1‑93 (MQ=255)
                             tGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTgg                                                            <  1:3763838/93‑1 (MQ=255)
                             tGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTgg                                                            <  1:5628946/93‑1 (MQ=255)
                              ggCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTggg                                                           >  2:731107/1‑93 (MQ=255)
                               gCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTgggg                                                          <  2:5332096/93‑1 (MQ=255)
                               gCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTgggg                                                          <  1:4632598/93‑1 (MQ=255)
                                 aGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGaa                                                        >  1:718925/1‑93 (MQ=255)
                                 aGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGaa                                                        >  2:3028479/1‑93 (MQ=255)
                                 aGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGaa                                                        <  1:2627968/93‑1 (MQ=255)
                                 aGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGaa                                                        <  2:4880641/93‑1 (MQ=255)
                                   atatTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATc                                                      <  2:4983993/93‑1 (MQ=255)
                                   atatTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATc                                                      >  2:4746481/1‑93 (MQ=255)
                                   atatTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATc                                                      <  1:5924757/93‑1 (MQ=255)
                                    tatTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCg                                                     >  1:3016194/1‑93 (MQ=255)
                                    tatTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCg                                                     >  2:1135315/1‑93 (MQ=255)
                                     atTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGa                                                    <  1:3215923/93‑1 (MQ=255)
                                      tttttCATCCGCTAATTTGATCCAATAAATAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGat                                                   <  1:4157655/93‑1 (MQ=255)
                                       ttttCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGata                                                  <  2:2511535/93‑1 (MQ=255)
                                       ttttCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGata                                                  <  2:3429005/93‑1 (MQ=255)
                                       ttttCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGata                                                  >  2:1354025/1‑93 (MQ=255)
                                        tttCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGatat                                                 <  1:3499878/93‑1 (MQ=255)
                                            aTCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGTACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATAttttt                                             <  1:3945949/93‑1 (MQ=255)
                                            aTCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATAttttt                                             >  2:5764386/1‑93 (MQ=255)
                                            aTCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATAttttt                                             <  2:837918/93‑1 (MQ=255)
                                              ccGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTtga                                           <  2:2422443/93‑1 (MQ=255)
                                               cGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTtgat                                          >  2:1562112/1‑93 (MQ=255)
                                               cGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTtgat                                          >  1:5050389/1‑93 (MQ=255)
                                                gCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTtgatg                                         >  1:3311678/1‑93 (MQ=255)
                                                gCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTtgatg                                         <  1:5748190/93‑1 (MQ=255)
                                                 cTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACC‑CTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTtgatga                                        >  1:4429116/1‑93 (MQ=255)
                                                  tAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGAc                                       >  2:4578472/1‑93 (MQ=255)
                                                   aaTTTGATCGATTAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACa                                      <  2:5643626/93‑1 (MQ=255)
                                                   aaTTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACa                                      >  1:1916036/1‑93 (MQ=255)
                                                    aTTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGGATCGATATTTTTGATGACAt                                     >  2:2617586/1‑93 (MQ=255)
                                                     tttGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATa                                    >  2:4275115/1‑93 (MQ=255)
                                                      ttGATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATaa                                   <  1:1637521/93‑1 (MQ=255)
                                                        gATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGc                                 <  2:3139101/93‑1 (MQ=255)
                                                        gATCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGc                                 >  1:341890/1‑93 (MQ=255)
                                                          tCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAg                               <  1:6142045/93‑1 (MQ=255)
                                                          tCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAg                               <  1:2618050/93‑1 (MQ=255)
                                                          tCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAg                               >  2:139339/1‑93 (MQ=255)
                                                          tCGAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAg                               >  1:1135408/1‑93 (MQ=255)
                                                            gAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAggt                             >  1:5257422/1‑92 (MQ=255)
                                                            gAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGa                             >  1:4826104/1‑93 (MQ=255)
                                                            gAATAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGa                             >  2:4200450/1‑93 (MQ=255)
                                                              aTAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGAtt                           <  2:1617372/93‑1 (MQ=255)
                                                              aTAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGAtt                           >  1:1032690/1‑93 (MQ=255)
                                                               tAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGAttt                          >  2:3985832/1‑93 (MQ=255)
                                                               tAACTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGAttt                          <  2:462934/93‑1 (MQ=255)
                                                                  cTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGc                       >  2:1112313/1‑93 (MQ=255)
                                                                  cTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGc                       <  2:5064858/93‑1 (MQ=255)
                                                                  cTAATACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGc                       <  2:6076345/93‑1 (MQ=255)
                                                                    aaTACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGGATCGATATTTTTGATGACATAAGCAGGATTTAGCTc                     >  1:1313994/1‑93 (MQ=255)
                                                                     agccGGTACTCTGATGAGGAGCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGcaca                    <  2:172030/90‑1 (MQ=255)
                                                                     aTACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTca                    >  2:3664083/1‑93 (MQ=255)
                                                                     aTACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTca                    <  2:5724028/93‑1 (MQ=255)
                                                                      tACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTcac                   >  1:6086120/1‑93 (MQ=255)
                                                                      tACGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTcac                   >  1:3569520/1‑93 (MQ=255)
                                                                       aCGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTcaca                  <  2:5818055/93‑1 (MQ=255)
                                                                        cGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCacac                 >  2:2306366/1‑93 (MQ=255)
                                                                        cGGTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCacac                 <  1:1897189/93‑1 (MQ=255)
                                                                         ggTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACt                <  1:2270517/93‑1 (MQ=255)
                                                                         ggTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACt                <  1:5769355/93‑1 (MQ=255)
                                                                         ggTTCTCTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACt                >  1:4733535/1‑93 (MQ=255)
                                                                            tctcTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTAt             >  2:5377454/1‑93 (MQ=255)
                                                                            tctcTGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTAt             >  1:4496161/1‑93 (MQ=255)
                                                                               ctGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGa          <  1:1582636/93‑1 (MQ=255)
                                                                               ctGATGAGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGa          <  2:4534795/93‑1 (MQ=255)
                                                                                tgatgaGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGAc         <  2:4436364/93‑1 (MQ=255)
                                                                                 gatgaGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACg        <  2:2788038/93‑1 (MQ=255)
                                                                                   tgaGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGt      >  1:734962/1‑93 (MQ=255)
                                                                                     aGGACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGa    <  2:6268330/93‑1 (MQ=255)
                                                                                       gACCG‑TTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAg  >  2:5732209/1‑93 (MQ=255)
                                                                                            |                                                                                        
GCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAG  >  NC_000913/34020‑34200

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