Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,950,296:1 +C 15.8% coding (227/1773 nt) aspS ← aspartyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,950,2951.C15.8% 485.9 / 97.2 139coding (228/1773 nt)aspSaspartyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base . (60/57);  new base C (13/9);  total (73/66)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

Only 100 of 141 total aligned reads displayed.

CGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCAC  >  NC_000913/1950196‑1950390
                                                                                                    |                                                                                               
cGGCAGAACATCTGCGCGGGTAATGATAGGCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCa                                                                                                  <  1:2294356/100‑1 (MQ=255)
   cAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGt                                                                                              >  2:2152339/1‑100 (MQ=255)
    aGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTa                                                                                             >  2:2524999/1‑100 (MQ=255)
       aCATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACg                                                                                           <  2:3122829/100‑1 (MQ=255)
          tCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGc                                                                                       <  1:622506/100‑1 (MQ=255)
           cTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGcc                                                                                      <  1:2768388/100‑1 (MQ=255)
           cTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGcc                                                                                      >  2:2859947/1‑100 (MQ=255)
            tGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGccc                                                                                     <  1:406110/100‑1 (MQ=255)
            tGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGccc                                                                                     >  2:1676055/1‑100 (MQ=255)
                cgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGa                                                                                 <  1:1572083/100‑1 (MQ=38)
                 ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGAc                                                                                <  1:227702/100‑1 (MQ=38)
                 ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTg                                                                                  <  1:2674960/98‑1 (MQ=255)
                 ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGAc                                                                                >  2:3540029/1‑100 (MQ=255)
                 ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGAc                                                                                <  2:3644544/100‑1 (MQ=255)
                 ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGAc                                                                                <  2:2944384/100‑1 (MQ=255)
                 ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGACCGTGa                                                                                 >  1:1332621/1‑99 (MQ=255)
                  gTTGATGATTGTCAGCGCGGACGCCAGCTCTTCGATTTCGCCGGTCGCCATATCGCGGTTGATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTgcc                                                                                <  1:2602981/99‑3 (MQ=255)
                  gTTGATGATAGTCAGCGAGGACGCCCGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGAcg                                                                               <  2:645759/100‑2 (MQ=37)
                  gTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGAcg                                                                               <  1:983416/100‑2 (MQ=38)
                  gTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGAcc                                                                               >  2:451763/1‑100 (MQ=255)
                    tgatgaTAGTCAGCGAGGACGCCAGCCCTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCTCGTGAcgaa                                                                             <  2:836261/100‑4 (MQ=37)
                    tgatgaTAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTg                                                                             <  2:2232549/100‑1 (MQ=255)
                       tgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTTTTTACTCTCACGTACGGTGCCCGTGACCTGAAt                                                                          >  1:2707293/1‑100 (MQ=255)
                       tgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGAcgaaaaa                                                                          <  2:3545802/100‑7 (MQ=25)
                       tgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAAt                                                                          >  1:2151737/1‑100 (MQ=255)
                        gatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATg                                                                         <  1:2864170/100‑1 (MQ=255)
                          tAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCa                                                                       <  2:2786555/100‑1 (MQ=255)
                          tAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTGCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCa                                                                       <  1:1199549/100‑1 (MQ=255)
                            gTCAGCGAGTACGCCAGCACTTCGTTTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCTTCACCCGCACGTACGGTGCCCGTGTCCTGAATGCAGa                                                                     <  1:2434712/100‑1 (MQ=255)
                            gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGCTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGa                                                                     <  2:615211/100‑1 (MQ=255)
                             tCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGaa                                                                    <  2:1671023/100‑1 (MQ=255)
                              cAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                   >  2:2958988/1‑100 (MQ=255)
                              cAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                   <  1:1611557/100‑1 (MQ=255)
                              cAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                   >  1:424946/1‑100 (MQ=255)
                              cAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                   >  1:1264691/1‑100 (MQ=255)
                              cAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                   <  2:1978231/100‑1 (MQ=255)
                               aGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACt                                                                  >  2:2634753/1‑100 (MQ=255)
                               aGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACt                                                                  >  1:1469474/1‑100 (MQ=255)
                                gCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTc                                                                 >  2:89845/1‑100 (MQ=255)
                                    ggACGCCAGCACTTCGATTTCGCCGGTCGCCCTATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTa                                                             >  1:913890/1‑100 (MQ=255)
                                    ggACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTa                                                             >  2:1061371/1‑100 (MQ=255)
                                         ccAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAg                                                        <  2:3593434/100‑1 (MQ=255)
                                         ccAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAg                                                        >  2:1430446/1‑100 (MQ=255)
                                           aGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGtt                                                      >  2:939987/1‑100 (MQ=255)
                                               cTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCaga                                                  >  2:1235072/1‑100 (MQ=255)
                                               cTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCaga                                                  >  2:556717/1‑100 (MQ=255)
                                                ttCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCagag                                                 <  1:3759021/100‑1 (MQ=255)
                                                ttCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTATAAGCAGAACTCATTACGCAGCTccaa                                                  >  1:2705242/1‑96 (MQ=255)
                                                 tCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAgg                                                >  1:1929744/1‑100 (MQ=255)
                                                 tCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAgg                                                >  1:1929759/1‑100 (MQ=255)
                                                 tCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAgg                                                >  2:3799814/1‑100 (MQ=255)
                                                 tCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGGCACGCGCACGTACGGTGCCCGTAACCTGAATACACAACTCACTCCGCAGTACAGAgg                                                >  1:2253582/1‑100 (MQ=255)
                                                    aTTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTAGGCAGTTCAGAGGcc                                              >  2:3811663/1‑99 (MQ=255)
                                                    aTTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCa                                             >  1:3665881/1‑100 (MQ=255)
                                                     tttCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAg                                            <  2:3062018/100‑1 (MQ=255)
                                                        cGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCtt                                         <  2:1264699/100‑1 (MQ=255)
                                                        cGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCtt                                         >  2:2810014/1‑100 (MQ=255)
                                                        cGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCtt                                         >  1:2378886/1‑100 (MQ=255)
                                                        cGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCtt                                         >  2:3574017/1‑100 (MQ=255)
                                                        cGCCGGTCGCCATATCGCGGTGAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCtt                                         <  1:3366903/100‑1 (MQ=255)
                                                         gCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCttt                                        >  2:970118/1‑100 (MQ=255)
                                                          ccGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTa                                       <  2:3418716/100‑1 (MQ=255)
                                                          ccGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTa                                       <  2:131710/100‑1 (MQ=255)
                                                          ccGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACCGTGCCCGTGACCTGAATGCATAACTCATTACCCAGCTCCGAGCCCACCttt                                        >  1:3040107/1‑99 (MQ=255)
                                                           cGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTaa                                      <  1:1110149/100‑1 (MQ=255)
                                                           cGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTaa                                      <  1:1110263/100‑1 (MQ=255)
                                                            ggTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAc                                     <  1:733470/100‑1 (MQ=255)
                                                              tCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcgc                                   >  2:2383120/1‑100 (MQ=255)
                                                               cGCCATATCGCGGTTAATATTTTTTTCTTCTTTCCCACGTACGGTGCCCGTCACCTGACTGCAGAACTCCTTACGCAGTTCAGCGGCCAGCTTTAAccc                                   >  1:1537770/1‑97 (MQ=255)
                                                                 ccATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTc                                >  2:1269036/1‑100 (MQ=255)
                                                                     atCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCa                            >  2:2197672/1‑100 (MQ=255)
                                                                      tCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCAc                           <  2:3330697/100‑1 (MQ=255)
                                                                        gcgGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGa                         >  2:851064/1‑100 (MQ=255)
                                                                         cgGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGAt                        <  2:582330/100‑1 (MQ=255)
                                                                         cgGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGAt                        >  1:2934203/1‑100 (MQ=255)
                                                                           gTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATcc                      <  2:3248375/100‑1 (MQ=255)
                                                                                tatTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATc                 >  1:1442187/1‑100 (MQ=255)
                                                                                 atTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCg                <  2:2049982/100‑1 (MQ=255)
                                                                                   nattcaggTCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaa              >  2:2364931/8‑100 (MQ=25)
                                                                                    tttttCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaa             <  2:2425119/100‑1 (MQ=255)
                                                                                     ttttCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaaa            >  2:26639/1‑100 (MQ=255)
                                                                                     ttttCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaaa            <  2:48165/100‑1 (MQ=255)
                                                                                      tttCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaaaa           <  2:643264/100‑1 (MQ=255)
                                                                                         ggTCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCTGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc        >  1:3574708/2‑100 (MQ=37)
                                                                                         ggTCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc        >  2:2193040/2‑100 (MQ=38)
                                                                                         ggTCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc        >  2:1256916/2‑100 (MQ=38)
                                                                                          gTCTCGCGCACGTACGGTGCCCGTGACTTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc       <  2:3764386/100‑1 (MQ=255)
                                                                                          gTCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc       >  2:3623442/1‑100 (MQ=38)
                                                                                          gTCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc       >  2:2514282/1‑100 (MQ=38)
                                                                                          gTCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc       >  2:951715/1‑100 (MQ=38)
                                                                                          gTCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc       >  2:35932/1‑100 (MQ=38)
                                                                                          gTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc       <  2:2530234/100‑1 (MQ=255)
                                                                                           tCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt      >  1:786311/1‑100 (MQ=38)
                                                                                            cACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTg     >  2:969892/1‑100 (MQ=255)
                                                                                             aCGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGGGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGc    <  1:2152324/100‑1 (MQ=255)
                                                                                             aCGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATGGAAAAATACCTGc    >  1:3389285/1‑100 (MQ=255)
                                                                                              cgcgCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCa   >  2:134575/1‑100 (MQ=255)
                                                                                               gcgcACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCAc  <  1:199639/100‑1 (MQ=255)
                                                                                                ngcACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTCAACGCGTCCGCACGATCCGGATCGAAAAATACCTGccc  >  1:3022361/2‑97 (MQ=255)
                                                                                                 gcACGTACGGTGCCCGTGCCCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTTCAc  <  1:2591600/98‑1 (MQ=255)
                                                                                                    |                                                                                               
CGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCAC  >  NC_000913/1950196‑1950390

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: