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.

CGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATA  >  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)
 ggCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑c                                                                                                    >  2:2064191/1‑100 (MQ=255)
    aGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTa                                                                                                 >  2:232626/1‑100 (MQ=255)
    aGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTa                                                                                                 >  1:3702149/1‑100 (MQ=255)
     gAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTAc                                                                                                >  2:519116/1‑100 (MQ=255)
     gAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTAc                                                                                                <  1:2571297/100‑1 (MQ=255)
      aaCATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACg                                                                                               <  1:3549053/100‑1 (MQ=255)
        cATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGt                                                                                             <  2:2137949/100‑1 (MQ=255)
           cTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGcc                                                                                          >  1:1824408/1‑100 (MQ=255)
            tGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACgtgcgc                                                                                          >  2:1933552/1‑95 (MQ=255)
            tGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGccc                                                                                         >  1:3230667/1‑100 (MQ=255)
            tGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCACA‑CGTACGGTGccc                                                                                         <  1:514527/100‑1 (MQ=255)
              cgcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGt                                                                                       >  2:2048358/1‑100 (MQ=255)
               gcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTg                                                                                      >  1:1037867/1‑100 (MQ=255)
               gcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTg                                                                                      <  2:2903707/100‑1 (MQ=255)
                cgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCCCGCGCA‑CGTACGGTGCCCGTGa                                                                                     >  1:3098593/1‑100 (MQ=255)
                cgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGa                                                                                     <  2:721523/100‑1 (MQ=255)
                 ggTTGATGATAGTCAGCGAGGCCGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGGTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGAc                                                                                    <  2:1760657/100‑1 (MQ=255)
                 ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAc                                                                                    <  1:1904447/100‑1 (MQ=38)
                 ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAc                                                                                    <  1:3255183/100‑1 (MQ=38)
                 ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAc                                                                                    <  1:1933561/100‑1 (MQ=38)
                 ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTg                                                                                      <  2:2531808/98‑1 (MQ=255)
                 ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGAc                                                                                    <  1:2078547/100‑1 (MQ=255)
                 ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGAc                                                                                    <  2:596708/100‑1 (MQ=255)
                    tgatgaTAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTg                                                                                 <  1:2491714/100‑1 (MQ=255)
                     gttgatAGTCAGCGAGGACGCCGCCGCTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAcgaaa                                                                                <  2:1742274/98‑5 (MQ=25)
                     gatgatAGTCAGCGAGGACGTCCGCACTTCGATTTTGCCGGTCGCCATAGCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAcgaa                                                                                 <  2:2545776/99‑4 (MQ=25)
                      atgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGaa                                                                               >  1:1045292/1‑100 (MQ=255)
                      atgatAGTCAACGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGTGGTTAATATTTTTGTCGTCACGCGCA‑CGTACGGTGCTCGTGACCTGaa                                                                               <  2:860050/100‑1 (MQ=255)
                         atAGTCGGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGc                                                                            <  2:3278303/100‑1 (MQ=255)
                          tAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCa                                                                           <  1:2027081/100‑1 (MQ=255)
                            gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGa                                                                         <  2:473431/100‑1 (MQ=255)
                            gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGa                                                                         >  2:2454837/1‑100 (MQ=255)
                            gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGa                                                                         <  1:2922119/100‑1 (MQ=255)
                            gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGa                                                                         >  2:3816176/1‑100 (MQ=255)
                            gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGa                                                                         >  2:3752629/1‑100 (MQ=255)
                              cAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                       <  2:3629495/100‑1 (MQ=255)
                              cAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                       <  2:2691568/100‑1 (MQ=255)
                              cAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCACA‑CGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                       >  2:3941767/1‑100 (MQ=255)
                               aGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                       <  2:273059/99‑1 (MQ=255)
                                gCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTc                                                                     >  1:2954564/1‑100 (MQ=255)
                                gCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTc                                                                     <  2:1578915/100‑1 (MQ=255)
                                 cGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCa                                                                    <  1:1358589/100‑1 (MQ=255)
                                     gACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTAc                                                                <  1:2455492/100‑1 (MQ=255)
                                      acaccAGCAACTCGATCTCGCCGGTCGCCATATCGAGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTAc                                                                <  1:2555943/96‑1 (MQ=255)
                                      aCGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGGA‑CGTACGGTGCCTGTGTCCTGAATGGAGAACTCATTACg                                                               <  2:1012886/100‑1 (MQ=255)
                                         ccAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAg                                                            >  2:3116657/1‑100 (MQ=255)
                                          cAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGt                                                           <  1:2637659/100‑1 (MQ=255)
                                               cTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTTGTTTTTTGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCaga                                                      >  2:3758738/1‑100 (MQ=255)
                                                ttCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCagag                                                     <  2:2232879/100‑1 (MQ=255)
                                                 tCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCagag                                                     <  2:3093620/99‑1 (MQ=255)
                                                   gATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGcc                                                  <  2:2009232/100‑1 (MQ=255)
                                                    aTTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCa                                                 >  1:2821569/1‑100 (MQ=255)
                                                    aTTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGGGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCa                                                 >  1:3762086/1‑100 (MQ=255)
                                                     tttCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAg                                                >  1:2401693/1‑100 (MQ=255)
                                                      ttCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGc                                               <  2:1113522/100‑1 (MQ=255)
                                                       tCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCt                                              >  1:2253771/1‑100 (MQ=255)
                                                            ggTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAc                                         <  2:3368156/100‑1 (MQ=255)
                                                              tCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGACCTCATTACGCAGTTCAGAGGCCAGCTTTAAcg                                        <  2:644389/99‑1 (MQ=255)
                                                              tCGCCATATCGCGGTTAATATTTTTTACGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGATCTCATTACGCAGTTCAGAGGCCAGCTTTAAc                                         <  2:1535272/98‑1 (MQ=255)
                                                              tCGCCATATCGCGGTTAATATCTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcgc                                       <  1:380982/100‑1 (MQ=255)
                                                               cGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcgcg                                      >  1:1120353/1‑100 (MQ=255)
                                                                gCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGGATGCAGAAATCATTACGAAGTTCCGAGGCCAGCTTTAAcgag                                      >  2:3272916/1‑97 (MQ=255)
                                                                gCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGt                                     >  1:199924/1‑100 (MQ=255)
                                                                   atatCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCg                                  <  2:682535/100‑1 (MQ=255)
                                                                     atCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCa                                <  1:2698016/100‑1 (MQ=255)
                                                                     atCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCa                                <  1:2096680/100‑1 (MQ=255)
                                                                       cgcgGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACg                              <  2:1037861/100‑1 (MQ=255)
                                                                        gcgGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGa                             >  1:801275/1‑100 (MQ=255)
                                                                        gcgGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGa                             >  1:861467/1‑100 (MQ=255)
                                                                         cgGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGAt                            <  1:2794422/100‑1 (MQ=255)
                                                                         cgGTTAATATTTTTTCCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATCCAGAACTCATTACGGAGTTCAGAGGCCAGCTTTAACACGTCAGCACGa                             >  2:3632530/1‑99 (MQ=255)
                                                                          ggTTAATATTTTTTTCGTCCCGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGGGGCCAGCTTTAACGCGTCCCCACGacc                           >  2:955989/1‑98 (MQ=255)
                                                                           gTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATcc                          >  1:3461380/1‑100 (MQ=255)
                                                                             tAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCgg                        <  2:1411997/100‑1 (MQ=255)
                                                                               atatTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGAt                      >  1:1237901/1‑100 (MQ=255)
                                                                                 atTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCg                    <  1:187308/100‑1 (MQ=255)
                                                                                  tttttttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGa                   <  2:2306498/100‑1 (MQ=255)
                                                                                  tttttttCGCCACGCGCA‑CTTCCGGTGCCCGTGCCCTGAATGCATAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATc                     <  2:445254/98‑1 (MQ=255)
                                                                                  tttgtttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAGCTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATc                     <  2:847513/98‑1 (MQ=255)
                                                                                   ttttttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaa                  <  1:3014589/100‑1 (MQ=255)
                                                                                        aggccACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCCGTTCAGAGGCCAGCTGTAACGCGTCCGCACGATCCGGATCGAAAAATa             >  1:1776018/5‑100 (MQ=37)
                                                                                         ggTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc            >  2:1341458/2‑100 (MQ=38)
                                                                                         ggTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc            >  2:1974390/2‑100 (MQ=38)
                                                                                         ggTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc            >  2:1543523/2‑100 (MQ=38)
                                                                                         ggTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCAGATCggaagagcg            >  2:2697374/2‑92 (MQ=25)
                                                                                         cGTCACGCGGA‑CGTACGGTGCCCGTGTCCTGAGTGCAGAACTCATTTCGCAGTCCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc            <  2:1487070/100‑1 (MQ=255)
                                                                                          gTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc           >  1:2731618/1‑100 (MQ=38)
                                                                                          gTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAcc           <  1:519119/100‑1 (MQ=255)
                                                                                           tCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt          >  2:3843539/1‑100 (MQ=38)
                                                                                           tCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt          >  2:2208452/1‑100 (MQ=38)
                                                                                           tCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt          >  1:672363/1‑100 (MQ=38)
                                                                                            cACGGGCACCGTAC‑GTGCGCGTGACCTGAGTGCAGAACCCATTACGCAGTGCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTg         >  1:2288822/1‑100 (MQ=255)
                                                                                            cACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTg         >  1:1204922/1‑100 (MQ=255)
                                                                                              cgcgCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCa       <  1:3929119/100‑1 (MQ=255)
                                                                                              cgcgCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCa       >  1:3717063/1‑100 (MQ=255)
                                                                                               gcgcA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCAc      <  1:3296111/100‑1 (MQ=255)
                                                                                                cgcA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACg     <  1:3541752/100‑1 (MQ=255)
                                                                                                    |                                                                                                   
CGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATA  >  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: ‑