Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,950,302 Δ1 bp 15.7% coding (221/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,3010G.15.7% 475.2 / 86.7 159coding (222/1773 nt)aspSaspartyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (66/68);  new base . (18/7);  total (84/75)
Fisher's exact test for biased strand distribution p-value = 4.89e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-01

Only 100 of 161 total aligned reads displayed.

ACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCG  >  NC_000913/1950203‑1950400
                                                                                                    |                                                                                                   
aCATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACgg                                                                                                    <  2:3705783/100‑1 (MQ=255)
 cATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGt                                                                                                   >  2:300048/1‑100 (MQ=255)
   tCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGc                                                                                                 <  1:1895261/100‑1 (MQ=255)
      gcgcGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCg                                                                                              >  1:3977982/1‑100 (MQ=255)
      gcgcGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCg                                                                                              <  2:3438354/100‑1 (MQ=255)
        gcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTg                                                                                            <  2:253426/100‑1 (MQ=255)
        gcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTg                                                                                            <  1:1998328/100‑1 (MQ=255)
          ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGAc                                                                                          <  2:4019513/100‑1 (MQ=38)
          ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTGCGTCACGGGCACCGTACGTGCGCGTGAc                                                                                          <  2:2622848/100‑1 (MQ=38)
           gTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGAcg                                                                                         <  2:3911784/100‑2 (MQ=38)
           gTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGtta                                                                                           >  1:1719680/1‑96 (MQ=255)
           gTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGAcc                                                                                         <  2:1072081/100‑1 (MQ=255)
             tgattaTCGTCAGCGAGGACGCCGGTAATTCGATTACGCGGGTCGCCATATCGCGGTTAATTTTTTTTTCGTCACGGGCTCCGTACGTGCGCGTGAcgaa                                                                                       <  1:99127/100‑4 (MQ=16)
             tgatgaTAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGAcgaa                                                                                       <  2:588472/100‑4 (MQ=37)
             tgatgaTAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGAcgaa                                                                                       <  2:4326837/100‑4 (MQ=37)
             tgatgaTAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTg                                                                                       <  2:1550810/100‑1 (MQ=255)
             tgatgaTAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTg                                                                                       <  1:3605120/100‑1 (MQ=255)
               atgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGaa                                                                                     >  1:1536735/1‑100 (MQ=255)
                tgatAGTCAGCGAGGACGCCAGAACTTCGATTTAGTCGGTGGCCATATCGCGGTTAATATTTTTTTCGTTACGCGCACGTAAGGTGGCCGTGACCTGaa                                                                                     >  1:2651104/1‑99 (MQ=255)
                 gatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGATCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATg                                                                                   >  1:2345422/1‑100 (MQ=255)
                  atAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGc                                                                                  >  1:2289389/1‑100 (MQ=255)
                   tAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGc                                                                                  <  1:2246495/99‑1 (MQ=255)
                   tAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCa                                                                                 <  2:3954423/100‑1 (MQ=255)
                   tAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCa                                                                                 <  2:2132309/100‑1 (MQ=255)
                   tAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCa                                                                                 >  1:2958989/1‑100 (MQ=255)
                    aggcAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAg                                                                                <  1:1051647/97‑1 (MQ=255)
                       cAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                             >  2:1598245/1‑100 (MQ=255)
                        aGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                             >  1:70664/1‑99 (MQ=255)
                        aGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACt                                                                            >  1:2160162/1‑100 (MQ=255)
                         gCGAGGACGCGGGCACGTCGGGTTGGTCGGTCGCCATATTGCGGTGAGTATTTTTTTCGGCACGCGCAGGTACGGTGCCCGTGACCTGAATGCAg                                                                                <  1:3760444/95‑1 (MQ=255)
                          cGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCa                                                                          <  2:2734524/100‑1 (MQ=255)
                          cGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCa                                                                          <  1:1058035/100‑1 (MQ=255)
                          cGAGGACGCCAGCACTTCGATATCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATACAGAACTCa                                                                          >  1:4264138/1‑100 (MQ=255)
                          cGAGGACGCCAGCACTTCGATATCGCCGGTAGCCATATCGCGGTTAATATTTTTTTCTTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCa                                                                          >  1:498005/1‑100 (MQ=255)
                                cGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGc                                                                    >  1:2573541/1‑100 (MQ=255)
                                cGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGc                                                                    >  2:260646/1‑100 (MQ=255)
                                 gCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCa                                                                   >  1:4067131/1‑100 (MQ=255)
                                  ccAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAg                                                                  <  2:958163/100‑1 (MQ=255)
                                   cAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGt                                                                 <  2:231818/100‑1 (MQ=255)
                                    aGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTa                                                                       >  1:1034375/1‑93 (MQ=255)
                                     gCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTc                                                               >  2:3637710/1‑100 (MQ=255)
                                          tCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAgg                                                          <  2:524176/100‑1 (MQ=255)
                                           cGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGc                                                         >  1:3759312/1‑100 (MQ=255)
                                           cGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGc                                                         <  2:3571630/100‑1 (MQ=255)
                                            gATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGcc                                                        >  1:128693/1‑100 (MQ=255)
                                            gATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGcc                                                        <  2:948377/100‑1 (MQ=255)
                                             aTTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCa                                                       <  1:562576/100‑1 (MQ=255)
                                              tttCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAg                                                      >  2:2203726/1‑100 (MQ=255)
                                              tttCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAg                                                      <  1:1161912/100‑1 (MQ=255)
                                                tCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGc                                                     >  1:529914/1‑99 (MQ=255)
                                                tCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCt                                                    <  1:2104144/100‑1 (MQ=255)
                                                 cGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCtt                                                   >  1:1694897/1‑100 (MQ=255)
                                                  gCCGGTCGCCATATCGCGGTTAATATTTTTTTCTTCACGTGCACGTACGGTGCCCGTGACCTCAATGCAGACCTCATTACGAAGTTCAGAGGCCAGCttt                                                  >  1:677124/1‑100 (MQ=255)
                                                  gCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCttt                                                  <  2:4067160/100‑1 (MQ=255)
                                                   ccGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTa                                                 >  1:2484468/1‑100 (MQ=255)
                                                     ggTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAc                                               <  2:3399719/100‑1 (MQ=255)
                                                      gTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcg                                              <  2:1155793/100‑1 (MQ=255)
                                                      gTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcg                                              >  1:3465046/1‑100 (MQ=255)
                                                        cGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcgcg                                            >  2:2259018/1‑100 (MQ=255)
                                                          ccATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTGAAcgcg                                            >  2:491899/1‑98 (MQ=255)
                                                            atatCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCg                                        <  2:902835/100‑1 (MQ=255)
                                                             tatCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGc                                       <  1:2396273/100‑1 (MQ=255)
                                                              atCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCa                                      >  1:1558931/1‑100 (MQ=255)
                                                                cgcgGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGCCCTGAATGCAGAACTCATTAAGCAGTTCAAAGGCCCGCTTTAACGCGTCCGCACg                                    >  2:2355527/1‑100 (MQ=255)
                                                                 gcgGTTAATATTTTTTTCGTCACGCGCCCGTACGGTGCCCGTAACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGCCCGCACGa                                   >  2:2656956/1‑100 (MQ=255)
                                                                    gTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATcc                                <  1:3922494/100‑1 (MQ=255)
                                                                       aaTATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGa                             >  2:2024282/1‑100 (MQ=255)
                                                                        atatTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGAt                            >  2:3021393/1‑100 (MQ=255)
                                                                          atTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCg                          <  2:1278957/100‑1 (MQ=255)
                                                                          atTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCg                          <  2:1279701/100‑1 (MQ=255)
                                                                           tttttttCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCCGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGa                         >  2:1692606/1‑100 (MQ=255)
                                                                             tttttCGTCACGCGCACGTACGGTGCCCGTGTCCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaa                       <  1:102825/100‑1 (MQ=255)
                                                                             tttttCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaa                       <  1:300047/100‑1 (MQ=255)
                                                                             tgtttCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGGCCGCACGATCCGGATCGaaa                       <  1:113921/98‑1 (MQ=255)
                                                                              ttttCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaaa                      <  2:2850803/100‑1 (MQ=255)
                                                                               tttCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaaaa                     >  1:878684/1‑100 (MQ=255)
                                                                                 tCGTCACGCGCACGTACGGTGCCCGGGACCTGAACGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGAACCGGATCGAAAAATa                   <  1:3183759/100‑1 (MQ=255)
                                                                                  ggTCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc                  >  2:41852/2‑100 (MQ=38)
                                                                                  ggTCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc                  >  2:1642417/2‑100 (MQ=38)
                                                                                  ggTCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc                  >  1:3219154/2‑100 (MQ=38)
                                                                                  cGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc                  >  1:3335356/1‑100 (MQ=255)
                                                                                  cGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGAGTCCGGACGATCAGGATCGAAAAATa                   >  1:1142083/1‑99 (MQ=255)
                                                                                   gTCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc                 >  1:237733/1‑100 (MQ=38)
                                                                                   gTCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc                 >  1:2483560/1‑100 (MQ=38)
                                                                                    tCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt                >  2:4060655/1‑100 (MQ=38)
                                                                                    tCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGCCCGCACGATCCGGATCGAAAAATACCt                >  1:3833209/1‑100 (MQ=38)
                                                                                    tCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt                <  1:2216696/100‑1 (MQ=255)
                                                                                     cACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTg               >  2:2600675/1‑100 (MQ=255)
                                                                                     cACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTg               >  1:7800/1‑100 (MQ=255)
                                                                                     cACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTg               <  1:4092077/100‑1 (MQ=255)
                                                                                       cgcgCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCa             >  2:1413299/1‑100 (MQ=255)
                                                                                       acgggcaccGTACGTGCGCGTGACCTGAATGCAGCACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGc              >  2:2186332/9‑100 (MQ=255)
                                                                                       acgggcaccGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGc              >  2:3267905/9‑100 (MQ=255)
                                                                                         cgcACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACg           >  2:2855296/1‑100 (MQ=255)
                                                                                               cGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATAc       <  2:4351082/100‑1 (MQ=255)
                                                                                               cGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATAc       <  2:803741/100‑1 (MQ=255)
                                                                                                 tACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGGGTCCGCACGATCCGGATCGAAAAAGACCTGCACGAtgc       >  1:2677704/1‑96 (MQ=255)
                                                                                                 tACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATAcc      <  1:1833539/99‑1 (MQ=255)
                                                                                                 tACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCt     >  2:1960304/1‑100 (MQ=255)
                                                                                                    ggTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTcg  >  2:4222828/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
ACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCG  >  NC_000913/1950203‑1950400

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