Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,411,512 G→A 20.8% A146V (GCG→GTG)  ttcA ← tRNA s(2)C32 thioltransferase, iron sulfur cluster protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,411,5120GA20.8% 206.1 / 45.0 96A146V (GCG→GTG) ttcAtRNA s(2)C32 thioltransferase, iron sulfur cluster protein
Reads supporting (aligned to +/- strand):  ref base G (47/29);  new base A (11/9);  total (58/38)
Fisher's exact test for biased strand distribution p-value = 6.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

CATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTC  >  NC_000913/1411414‑1411610
                                                                                                  |                                                                                                  
catcGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAAcgc                                                                                                   >  2:2301319/1‑100 (MQ=255)
catcGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAAcgc                                                                                                   >  1:389748/1‑100 (MQ=255)
  tcGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGa                                                                                                 >  1:1363165/1‑100 (MQ=255)
  tcGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGa                                                                                                 <  1:1593290/100‑1 (MQ=255)
   ccgatctcaGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGAt                                                                                                <  1:1983531/94‑1 (MQ=255)
    gCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATc                                                                                               >  2:2343129/1‑100 (MQ=255)
    gCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACACGATc                                                                                               >  2:531303/1‑100 (MQ=255)
     cTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCt                                                                                              <  1:2176417/100‑1 (MQ=255)
     cTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCt                                                                                              >  1:1989481/1‑100 (MQ=255)
     cTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCt                                                                                              <  2:1817151/100‑1 (MQ=255)
      tcatcaGTTTCGTAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACACGATCt                                                                                              <  1:1886167/99‑1 (MQ=255)
       catcaGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACACGATCTtc                                                                                            >  1:1472385/1‑100 (MQ=255)
       catcaGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGCCACGATGGTGACCCAACGCGATCTtc                                                                                            >  1:2258041/1‑100 (MQ=255)
        atcaGTTTCGGAGTCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTtcg                                                                                           <  1:2534063/100‑1 (MQ=255)
        atcaGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTtcg                                                                                           >  1:7438/1‑100 (MQ=255)
         tcaGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTtcgt                                                                                          >  1:890634/1‑100 (MQ=255)
          caGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAATAACAACGTTTGCAGTATATCGTCACGATGGTGACCCAACGCGATCTtcgtc                                                                                         <  2:623429/100‑1 (MQ=255)
          caGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGAtagatcgg                                                                                         >  2:1983512/1‑93 (MQ=255)
          caGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACACGATCTtcgtc                                                                                         >  2:1311213/1‑100 (MQ=255)
           aGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTtcgtcg                                                                                        >  2:587057/1‑100 (MQ=255)
           aGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACACGATCTtcgtcg                                                                                        <  2:36650/100‑1 (MQ=255)
           aGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACACGATCTtcgtcg                                                                                        <  1:1339659/100‑1 (MQ=255)
            gTTTCAGAGGCATACCTTTCATCTTACAGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGAGCTtcgtcg                                                                                        >  2:1859993/1‑99 (MQ=255)
                cGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGccccc                                                                                   >  1:986174/1‑100 (MQ=255)
                cGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGccccc                                                                                   >  1:962954/1‑100 (MQ=255)
                      cATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACACGATCTTCGTCGCCCCCAGTTcc                                                                             >  2:392385/1‑100 (MQ=255)
                      cATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGAATATAATCACGATAGTGACCCAACGCGATCTTCGTCGCCCCCAGTTcc                                                                             >  1:27010/1‑100 (MQ=255)
                       aTACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCg                                                                            <  2:2218418/100‑1 (MQ=255)
                        tACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGt                                                                           >  1:2149560/1‑100 (MQ=255)
                           cTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGc                                                                        >  2:317449/1‑100 (MQ=255)
                            tttCATCTTACCGTCGTAGTACCTATTTAAGTCCATCGTTTTCAGGTTATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCg                                                                       <  1:47976/100‑1 (MQ=255)
                             ttCATCTTACTTCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCgg                                                                      <  1:2461186/100‑1 (MQ=255)
                             ttCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACACGATCTTCGTCGCCCCCAGTTCCGTTGCgg                                                                      <  1:2643262/100‑1 (MQ=255)
                               cATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTa                                                                    >  2:876818/1‑100 (MQ=255)
                               cATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGCTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTa                                                                    >  2:1861765/1‑100 (MQ=255)
                                aTCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACACGATCTTCGTCGCCCCCAGTTCCGTTGCGGTAc                                                                   <  2:1472371/100‑1 (MQ=255)
                                 tCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACg                                                                  <  2:123692/100‑1 (MQ=255)
                                    tACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATa                                                               >  1:277814/1‑100 (MQ=255)
                                    tACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATa                                                               >  2:2547710/1‑100 (MQ=255)
                                    tACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTACCGTTGCGGTACGATa                                                               >  2:1545313/1‑100 (MQ=255)
                                      ccgccgTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATaaa                                                             >  2:358504/1‑100 (MQ=255)
                                      ccgccgTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATaaa                                                             >  2:885907/1‑100 (MQ=255)
                                        gccgTAGCACACATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAgg                                                           >  2:1385301/1‑100 (MQ=255)
                                        gccgTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAgg                                                           >  1:903329/1‑100 (MQ=255)
                                        gccgTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAgg                                                           <  2:1033078/100‑1 (MQ=255)
                                         ccgTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGa                                                          <  2:820873/100‑1 (MQ=255)
                                          cgTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGGGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGAt                                                         <  2:102371/100‑1 (MQ=255)
                                           gTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATa                                                        <  2:1875770/100‑1 (MQ=255)
                                           gTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACACGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATa                                                        <  2:1241832/100‑1 (MQ=255)
                                             aGAACATATTTAAGAACAACGTTTGCTGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATAcc                                                      <  1:1369953/100‑1 (MQ=255)
                                             aGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATAcc                                                      <  2:1383134/100‑1 (MQ=255)
                                              gAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACcg                                                     >  1:172329/1‑100 (MQ=255)
                                               aaCATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACcgc                                                    <  1:1902097/100‑1 (MQ=255)
                                               aaCATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACcgc                                                    <  1:198167/100‑1 (MQ=255)
                                               aaCATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACcgc                                                    <  1:1839128/100‑1 (MQ=255)
                                                aCATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACcgcg                                                   >  2:1134079/1‑100 (MQ=255)
                                                 cATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGcga                                                  <  2:389746/100‑1 (MQ=255)
                                                   tatTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGcgacg                                                >  1:180989/1‑100 (MQ=255)
                                                   tatTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACACGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGcgacg                                                >  1:2401476/1‑100 (MQ=255)
                                                    atTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGcgacga                                               >  2:2279949/1‑100 (MQ=255)
                                                     tttAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGaa                                              >  2:283678/1‑100 (MQ=255)
                                                     tttAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGaa                                              >  1:1326394/1‑100 (MQ=255)
                                                     tttAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGaa                                              >  1:1326287/1‑100 (MQ=255)
                                                       tAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAgg                                            >  1:2126295/1‑100 (MQ=255)
                                                       tAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAgg                                            >  1:900309/1‑100 (MQ=255)
                                                        aaGAAAAACGTTTGCAGGATATCGTCACGATGGTGACCCAATGCGAGCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGc                                           <  2:840918/100‑1 (MQ=255)
                                                         aGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCg                                          >  1:2282394/1‑100 (MQ=255)
                                                          gAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCga                                         >  1:1318511/1‑100 (MQ=255)
                                                          gAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCga                                         >  2:2232150/1‑100 (MQ=255)
                                                            acaacGTTTGCAGGATATCGTCACGATGGTGACCCAACACGATCTTCGTAGCCCACAGTAACGTTAAGATAAGATAAAGGATAAAGAGACGAAGGCgaga                                       >  1:1414442/1‑100 (MQ=255)
                                                             caacGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGaa                                      >  1:432909/1‑100 (MQ=255)
                                                               acGTTTGCAGGATATCGTCACGATGGTGACCCAACACGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAaca                                    >  1:2038812/1‑100 (MQ=255)
                                                                cGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAac                                     <  1:550812/98‑1 (MQ=255)
                                                                    tGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTg                               <  1:170548/100‑1 (MQ=255)
                                                                    tGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTg                               >  1:843052/1‑100 (MQ=255)
                                                                    tGCAGGATATCGTCACGATGGTGACCCAACACGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTg                               <  2:1478807/100‑1 (MQ=255)
                                                                     gCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTg                               >  2:2075495/1‑99 (MQ=255)
                                                                     gCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGa                              <  1:208888/100‑1 (MQ=255)
                                                                        gcgttgggTCACCATCGTGACCCAACGCGATCTTCGTCGCACCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCa                           >  1:1544116/8‑100 (MQ=255)
                                                                            atCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTg                       <  1:1701493/100‑1 (MQ=255)
                                                                              cGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGt                     >  1:76555/1‑100 (MQ=255)
                                                                              cGTCACGATGGTGACCCAACGCGATCTTCGTAGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAAAACAGTGAGCAAGTg                       >  2:2386397/1‑98 (MQ=255)
                                                                               ttCACGATGGTGACACAACACGATCTTCGTCGCCCCCAGTTCCGTTGCGGAACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGtt                    >  2:591692/2‑100 (MQ=255)
                                                                                   cGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGc                >  2:982270/1‑100 (MQ=255)
                                                                                   cGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGc                >  2:2017060/1‑100 (MQ=255)
                                                                                   cGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGc                >  2:138362/1‑100 (MQ=255)
                                                                                    gATGGTGACCCAACACGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGcc               <  1:1926144/100‑1 (MQ=255)
                                                                                     aTGGTGACCCAACACGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGccc              <  1:1237766/100‑1 (MQ=255)
                                                                                      tGGTGATCCAACGGTATATTTGTCGCCCCCAGTAACGTTGAGGTACGATAAAGGATACTGCGCCGAAGGCGCGATCACAATGAGCAAGTGGTTTTGCCct             >  2:1578399/1‑100 (MQ=255)
                                                                                        gTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCctct           <  2:2091880/100‑1 (MQ=255)
                                                                                          gACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTgg         <  1:1545328/100‑1 (MQ=255)
                                                                                          gACCCAACACGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTgg         >  1:1156777/1‑100 (MQ=255)
                                                                                           aCCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGa        <  2:463706/100‑1 (MQ=255)
                                                                                            cccAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGaa       >  1:436063/1‑100 (MQ=255)
                                                                                            cccAACACGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGaa       >  1:615868/1‑100 (MQ=255)
                                                                                            cccAACACGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGaa       >  1:1952410/1‑100 (MQ=255)
                                                                                               aaCGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAAtct    >  2:2496086/1‑100 (MQ=255)
                                                                                               aaCGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAAtct    <  1:2398663/100‑1 (MQ=255)
                                                                                                 cgcgATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAAtcttc  >  1:917083/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
CATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTC  >  NC_000913/1411414‑1411610

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: