Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,950,296:1 +C 17.6% 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.C17.6% 495.4 / 119.3 153coding (228/1773 nt)aspSaspartyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base . (58/68);  new base C (20/7);  total (78/75)
Fisher's exact test for biased strand distribution p-value = 1.04e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

Only 100 of 158 total aligned reads displayed.

CGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATA  >  NC_000913/1950196‑1950394
                                                                                                    |                                                                                                   
cGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCa                                                                                                      >  1:2333335/1‑100 (MQ=255)
cGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCa                                                                                                      <  1:1509351/100‑1 (MQ=255)
cGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCa                                                                                                      >  1:1554325/1‑100 (MQ=255)
 ggCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCAc                                                                                                    >  2:2064191/1‑100 (MQ=255)
    aGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTa                                                                                                 >  2:232626/1‑100 (MQ=255)
    aGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTa                                                                                                 >  1:3702149/1‑100 (MQ=255)
     gAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTAc                                                                                                >  2:519116/1‑100 (MQ=255)
     gAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTAc                                                                                                <  1:2571297/100‑1 (MQ=255)
      aaCATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACg                                                                                               <  1:3549053/100‑1 (MQ=255)
        cATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGt                                                                                             <  2:2137949/100‑1 (MQ=255)
           cTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGcc                                                                                          >  1:1824408/1‑100 (MQ=255)
            tGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACgtgcgc                                                                                          >  2:1933552/1‑95 (MQ=255)
            tGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGccc                                                                                         >  1:3230667/1‑100 (MQ=255)
            tGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCACACGTACGGTGccc                                                                                         <  1:514527/100‑1 (MQ=255)
              cgcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGt                                                                                       >  2:2048358/1‑100 (MQ=255)
               gcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTg                                                                                      >  1:1037867/1‑100 (MQ=255)
               gcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTg                                                                                      <  2:2903707/100‑1 (MQ=255)
                cgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCCCGCGCACGTACGGTGCCCGTGa                                                                                     >  1:3098593/1‑100 (MQ=255)
                cgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGa                                                                                     <  2:721523/100‑1 (MQ=255)
                 ggTTGATGATAGTCAGCGAGGCCGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGGTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGAc                                                                                    <  2:1760657/100‑1 (MQ=255)
                 ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGAc                                                                                    <  1:1904447/100‑1 (MQ=38)
                 ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGAc                                                                                    <  1:3255183/100‑1 (MQ=38)
                 ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGAc                                                                                    <  1:1933561/100‑1 (MQ=38)
                 ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTg                                                                                      <  2:2531808/98‑1 (MQ=255)
                 ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGAc                                                                                    <  1:2078547/100‑1 (MQ=255)
                 ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGAc                                                                                    <  2:596708/100‑1 (MQ=255)
                    tgatgaTAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTg                                                                                 <  1:2491714/100‑1 (MQ=255)
                     gttgatAGTCAGCGAGGACGCCGCCGCTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGAcgaaa                                                                                <  2:1742274/98‑5 (MQ=25)
                     gatgatAGTCAGCGAGGACGTCCGCACTTCGATTTTGCCGGTCGCCATAGCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGAcgaa                                                                                 <  2:2545776/99‑4 (MQ=25)
                      atgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGaa                                                                               >  1:1045292/1‑100 (MQ=255)
                      atgatAGTCAACGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGTGGTTAATATTTTTGTCGTCACGCGCACGTACGGTGCTCGTGACCTGaa                                                                               <  2:860050/100‑1 (MQ=255)
                         atAGTCGGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGc                                                                            <  2:3278303/100‑1 (MQ=255)
                          tAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCa                                                                           <  1:2027081/100‑1 (MQ=255)
                            gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGa                                                                         <  2:473431/100‑1 (MQ=255)
                            gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGa                                                                         >  2:2454837/1‑100 (MQ=255)
                            gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGa                                                                         <  1:2922119/100‑1 (MQ=255)
                            gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGa                                                                         >  2:3816176/1‑100 (MQ=255)
                            gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGa                                                                         >  2:3752629/1‑100 (MQ=255)
                              cAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                       <  2:3629495/100‑1 (MQ=255)
                              cAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                       <  2:2691568/100‑1 (MQ=255)
                              cAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCACACGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                       >  2:3941767/1‑100 (MQ=255)
                               aGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                       <  2:273059/99‑1 (MQ=255)
                                gCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTc                                                                     >  1:2954564/1‑100 (MQ=255)
                                gCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTc                                                                     <  2:1578915/100‑1 (MQ=255)
                                 cGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCa                                                                    <  1:1358589/100‑1 (MQ=255)
                                     gACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTAc                                                                <  1:2455492/100‑1 (MQ=255)
                                      acaccAGCAACTCGATCTCGCCGGTCGCCATATCGAGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTAc                                                                <  1:2555943/96‑1 (MQ=255)
                                      aCGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGGACGTACGGTGCCTGTGTCCTGAATGGAGAACTCATTACg                                                               <  2:1012886/100‑1 (MQ=255)
                                         ccAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAg                                                            >  2:3116657/1‑100 (MQ=255)
                                          cAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGt                                                           <  1:2637659/100‑1 (MQ=255)
                                               cTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTTGTTTTTTGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCaga                                                      >  2:3758738/1‑100 (MQ=255)
                                                ttCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCagag                                                     <  2:2232879/100‑1 (MQ=255)
                                                 tCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCagag                                                     <  2:3093620/99‑1 (MQ=255)
                                                   gATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGcc                                                  <  2:2009232/100‑1 (MQ=255)
                                                    aTTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCa                                                 >  1:2821569/1‑100 (MQ=255)
                                                    aTTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGGGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCa                                                 >  1:3762086/1‑100 (MQ=255)
                                                     tttCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAg                                                >  1:2401693/1‑100 (MQ=255)
                                                      ttCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGc                                               <  2:1113522/100‑1 (MQ=255)
                                                       tCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCt                                              >  1:2253771/1‑100 (MQ=255)
                                                            ggTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAc                                         <  2:3368156/100‑1 (MQ=255)
                                                              tCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGACCTCATTACGCAGTTCAGAGGCCAGCTTTAAcg                                        <  2:644389/99‑1 (MQ=255)
                                                              tCGCCATATCGCGGTTAATATTTTTTACGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGATCTCATTACGCAGTTCAGAGGCCAGCTTTAAc                                         <  2:1535272/98‑1 (MQ=255)
                                                              tCGCCATATCGCGGTTAATATCTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcgc                                       <  1:380982/100‑1 (MQ=255)
                                                               cGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcgcg                                      >  1:1120353/1‑100 (MQ=255)
                                                                gCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGGATGCAGAAATCATTACGAAGTTCCGAGGCCAGCTTTAAcgag                                      >  2:3272916/1‑97 (MQ=255)
                                                                gCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGt                                     >  1:199924/1‑100 (MQ=255)
                                                                   atatCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCg                                  <  2:682535/100‑1 (MQ=255)
                                                                     atCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCa                                <  1:2698016/100‑1 (MQ=255)
                                                                     atCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCa                                <  1:2096680/100‑1 (MQ=255)
                                                                       cgcgGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACg                              <  2:1037861/100‑1 (MQ=255)
                                                                        gcgGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGa                             >  1:801275/1‑100 (MQ=255)
                                                                        gcgGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGa                             >  1:861467/1‑100 (MQ=255)
                                                                         cgGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGAt                            <  1:2794422/100‑1 (MQ=255)
                                                                         cgGTTAATATTTTTTCCGTCACGCGCACGTACGGTGCCCGTGACCTGAATCCAGAACTCATTACGGAGTTCAGAGGCCAGCTTTAACACGTCAGCACGa                             >  2:3632530/1‑99 (MQ=255)
                                                                          ggTTAATATTTTTTTCGTCCCGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGGGGCCAGCTTTAACGCGTCCCCACGacc                           >  2:955989/1‑98 (MQ=255)
                                                                           gTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATcc                          >  1:3461380/1‑100 (MQ=255)
                                                                             tAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCgg                        <  2:1411997/100‑1 (MQ=255)
                                                                               atatTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGAt                      >  1:1237901/1‑100 (MQ=255)
                                                                                 atTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCg                    <  1:187308/100‑1 (MQ=255)
                                                                                  tttttttCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGa                   <  2:2306498/100‑1 (MQ=255)
                                                                                  tttttttCGCCACGCGCACTTCCGGTGCCCGTGCCCTGAATGCATAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATc                     <  2:445254/98‑1 (MQ=255)
                                                                                  tttgtttCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAGCTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATc                     <  2:847513/98‑1 (MQ=255)
                                                                                   ttttttCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaa                  <  1:3014589/100‑1 (MQ=255)
                                                                                        aggccACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCCGTTCAGAGGCCAGCTGTAACGCGTCCGCACGATCCGGATCGAAAAATa             >  1:1776018/5‑100 (MQ=37)
                                                                                         ggTCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc            >  2:1341458/2‑100 (MQ=38)
                                                                                         ggTCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc            >  2:1974390/2‑100 (MQ=38)
                                                                                         ggTCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc            >  2:1543523/2‑100 (MQ=38)
                                                                                         ggTCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCAGATCggaagagcg            >  2:2697374/2‑92 (MQ=25)
                                                                                         cGTCACGCGGACGTACGGTGCCCGTGTCCTGAGTGCAGAACTCATTTCGCAGTCCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc            <  2:1487070/100‑1 (MQ=255)
                                                                                          gTCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc           >  1:2731618/1‑100 (MQ=38)
                                                                                          gTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc           <  1:519119/100‑1 (MQ=255)
                                                                                           tCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt          >  2:3843539/1‑100 (MQ=38)
                                                                                           tCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt          >  2:2208452/1‑100 (MQ=38)
                                                                                           tCACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt          >  1:672363/1‑100 (MQ=38)
                                                                                            cACGGGCACCGTACGTGCGCGTGACCTGAGTGCAGAACCCATTACGCAGTGCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTg         >  1:2288822/1‑100 (MQ=255)
                                                                                            cACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTg         >  1:1204922/1‑100 (MQ=255)
                                                                                              cgcgCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCa       <  1:3929119/100‑1 (MQ=255)
                                                                                              cgcgCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCa       >  1:3717063/1‑100 (MQ=255)
                                                                                               gcgcACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCAc      <  1:3296111/100‑1 (MQ=255)
                                                                                                cgcACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACg     <  1:3541752/100‑1 (MQ=255)
                                                                                                    |                                                                                                   
CGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATA  >  NC_000913/1950196‑1950394

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: