Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 594,948 G→C 10.6% Q166E (CAG→GAG) ‡ cusR ← DNA‑binding transcriptional activator CusR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913594,9480GC10.6% 206.4 / 11.8 65Q166E (CAG→GAG) ‡cusRDNA‑binding transcriptional activator CusR
Reads supporting (aligned to +/- strand):  ref base G (31/27);  new base C (3/4);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCA  >  NC_000913/594808‑595097
                                                                                                                                            |                                                                                                                                                     
cggAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCgaag                                                                                                                                              <  2:1183891/148‑1 (MQ=255)
cTGAATTAGCTTCGTCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCgaag                                                                                                                                              <  2:1184318/150‑1 (MQ=255)
 tGAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGACCGGGGCAGCACTTCGCCCTGATGGCgaagg                                                                                                                                             >  1:1279069/1‑150 (MQ=255)
  gAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCgaagga                                                                                                                                            >  2:502790/1‑150 (MQ=255)
   aaTTATCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCATCAGGGCgaagtgc                                                                                                                                           <  2:614017/150‑4 (MQ=255)
   aaTTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCATCAGGGCgaagtgc                                                                                                                                           <  1:1010896/150‑4 (MQ=255)
   aaTTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCATCAGGGCgaagtgc                                                                                                                                           <  2:372411/150‑4 (MQ=255)
   aaTTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCATCAGGGCgaagtgc                                                                                                                                           <  2:488943/150‑4 (MQ=255)
    aTTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCgaaggaag                                                                                                                                          <  2:176166/150‑1 (MQ=255)
        gCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACt                                                                                                                                      <  2:1410759/150‑1 (MQ=255)
        gCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACt                                                                                                                                      >  2:1227739/1‑150 (MQ=255)
            cGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCcag                                                                                                                                  <  1:1400490/150‑1 (MQ=255)
            cGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCcag                                                                                                                                  <  1:48392/150‑1 (MQ=255)
                 cAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATAGCTGTCAAAATTCATGTCCCATACCTGAGAGGCGATAAGAGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCagag                                                                                                                             >  1:397631/1‑150 (MQ=255)
                     gTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTaaa                                                                                                                         >  2:1368566/1‑150 (MQ=255)
                        gTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAActc                                                                                                                      <  2:57522/150‑1 (MQ=255)
                           gTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGCTGGCGAAGGAAGAACTCCAGCAGAGTAAACTCttt                                                                                                                   >  2:821827/1‑150 (MQ=255)
                             cGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATTCCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTAc                                                                                                                 >  2:1128487/1‑150 (MQ=255)
                             cGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTAc                                                                                                                 >  2:1043552/1‑150 (MQ=255)
                               aTTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTg                                                                                                               <  1:120914/150‑1 (MQ=255)
                                  ttGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTc                                                                                                            <  1:1133175/150‑1 (MQ=255)
                                   tGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCa                                                                                                           <  2:263838/150‑1 (MQ=255)
                                   tGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCa                                                                                                           >  1:968566/1‑150 (MQ=255)
                                        cgcAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTg                                                                                                      >  2:1004462/1‑150 (MQ=255)
                                         gcAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGa                                                                                                     <  2:1146278/150‑1 (MQ=255)
                                           aGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATg                                                                                                   <  1:1227739/150‑1 (MQ=255)
                                                     ccGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGc                                                                                         >  2:433231/1‑150 (MQ=255)
                                                      cGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCt                                                                                        <  1:149593/150‑1 (MQ=255)
                                                       gCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTg                                                                                       >  2:736887/1‑150 (MQ=255)
                                                             tCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATg                                                                                                   >  1:782079/1‑132 (MQ=255)
                                                             tCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATg                                                                                                   <  2:782079/132‑1 (MQ=255)
                                                             tCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTg                                                                                 >  1:1359620/1‑150 (MQ=255)
                                                                   gCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTg                                                                           >  2:544589/1‑150 (MQ=255)
                                                                   gCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTg                                                                           >  1:513737/1‑150 (MQ=255)
                                                                    cATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTg                                                                                 >  1:461093/1‑143 (MQ=255)
                                                                    cATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTg                                                                                 <  2:461093/143‑1 (MQ=255)
                                                                     aTTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCg                                                                         <  1:137086/150‑1 (MQ=255)
                                                                        ggTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCt                                                                      >  1:548366/1‑150 (MQ=255)
                                                                          tATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGa                                                                    <  2:833341/150‑1 (MQ=255)
                                                                            tCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACg                                                                  <  2:446279/150‑1 (MQ=255)
                                                                            tCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGATGCGGGTGACTTTGCGGCTGACg                                                                  >  1:1202858/1‑150 (MQ=255)
                                                                                tGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGAt                                                              >  2:72099/1‑150 (MQ=255)
                                                                                 gTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATc                                                             >  2:590607/1‑150 (MQ=255)
                                                                                   cAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGa                                                           >  2:947857/1‑150 (MQ=255)
                                                                                   cAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGa                                                           >  2:312423/1‑150 (MQ=255)
                                                                                     aaaTTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGAcc                                                         <  1:502790/150‑1 (MQ=255)
                                                                                          cATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCaa                                                    >  2:1344341/1‑150 (MQ=255)
                                                                                             gTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATc                                                 >  2:1381641/1‑150 (MQ=255)
                                                                                               cccATACCTGCGAGGCTATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGAAGAGTAAACTCTTTACTGGTCAAAGTAATGCGCGTGCCGCTGCGGGTTACTTTGCGGCTGACGAGATCGACCATCAAATCgt                                               <  2:397631/150‑2 (MQ=255)
                                                                                               cccATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCgg                                               <  2:1359620/150‑1 (MQ=255)
                                                                                                 cATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCa                                             <  2:737056/150‑1 (MQ=255)
                                                                                                               gATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGAGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACt                               >  2:487608/1‑150 (MQ=255)
                                                                                                                   agcgagcgGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCg                           >  2:1319511/1‑150 (MQ=255)
                                                                                                                   agcgagcgGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCGGAAACTGACTTTCg                           >  1:623961/1‑150 (MQ=255)
                                                                                                                    gcgagcgGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGa                          <  1:72099/150‑1 (MQ=255)
                                                                                                                       agcgGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATaa                       <  1:1004462/150‑1 (MQ=255)
                                                                                                                          ggggCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCa                    <  1:736887/150‑1 (MQ=255)
                                                                                                                            gttcttccTTCGCCATCAGGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCAcc                  >  1:334363/8‑150 (MQ=255)
                                                                                                                             gcagcaCTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGCCCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCAccg                 >  1:1117238/1‑150 (MQ=255)
                                                                                                                               cttccTTCGCCATCAGGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCAccgcc               >  2:442535/5‑150 (MQ=255)
                                                                                                                               agcaCTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCAccgcc               >  2:219639/1‑150 (MQ=255)
                                                                                                                                  aCTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCgccgcc            <  1:1128487/150‑1 (MQ=255)
                                                                                                                                   cTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGccc           <  2:1307237/150‑1 (MQ=255)
                                                                                                                                   cTTCGCCATCAGGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGccc           >  2:320244/1‑150 (MQ=255)
                                                                                                                                      cGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCcgc        >  2:1235548/1‑150 (MQ=255)
                                                                                                                                            gATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCa  <  1:356887/150‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                     
CTGAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCA  >  NC_000913/594808‑595097

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: