Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,725,151 A→C 8.9% S20R (AGT→AGG)  PP_1518 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,725,1510AC 8.9% 173.1 / 7.6 67S20R (AGT→AGG) PP_1518hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (31/30);  new base C (2/4);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

GGCAGCTCGATGCGCACGTCCGGGTTGTTGCTGAAGCCGCCTGGGGTGCTCAGTTGCTTGACGGCCAGTTGCGCACCTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCTTGTAAACGGATGAGTACGGGTATCAGCGAACCTGGTCTACCTTGATGCGTACCGATTGCGGGTCACTGCCATCGA  >  NC_002947/1725002‑1725296
                                                                                                                                                     |                                                                                                                                                 
ggCAGCTCGATGCGCACGTCCGGGTTGTTGCTGAAGCCGCCTGGGGTGCTCAGTTGCTTGACGGCCAGTTGCGCACCTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCa                                                                                                                                                   <  1:850445/150‑1 (MQ=255)
   aGCTCGATGCGCACGTCCGGGTTGTTGCTGAAGCCGCCTGGGGTGCTCAGTTGCTTGACGGCCAGTTGCGCACCTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTg                                                                                                                                                >  1:422890/1‑150 (MQ=255)
    gCTCGATGCGCACGTCCGGGTTGTTGCTGAAGCCGCCTGGGGTGCTCAGTTGCTTGACGGCCAGTTGCGCACCTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTgg                                                                                                                                               >  1:463173/1‑150 (MQ=255)
     cTCGATGCGCACGTCCGGGTTGTTGCTGAAGCCGCCTGGGGTGCTCAGTTGCTTGACGGCCAGTTGCGCACCTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGc                                                                                                                                              <  2:625107/150‑1 (MQ=255)
        gATCCGCATTTCCGGGTTGTTGCTGAAGCCGCCTGGGGTGCTCAGTTGCTTGACGGCCAGTTGCGCACCTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGa                                                                                                                                           >  2:1349173/1‑150 (MQ=255)
          tGCGCACGTCCGGGTTGTTGCTGAAGCCGCCTGGGGTGCTCAGTTGCTTGACGGCCAGTTGCGCACCTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGTGTCCTGGccctgt                                                                                                                                         <  1:16384/150‑5 (MQ=255)
          tGCGCACGTCCGGGTTGTTGCTGAAGCCGCCTGGGGTGCTCAGTTGCTTGACGGCCAGTTGCGCACCTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGTGCCCTGGccctgt                                                                                                                                         <  2:770813/150‑5 (MQ=255)
          tGCGCACGTCCGGGTTGTTGCTGAAGCCGCCTGGGGTGCTCAGTTGCTTGACGGCCAGTTGCGCACCTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGTGCCCTGGccctgt                                                                                                                                         <  1:213657/150‑5 (MQ=255)
                 gTCCGGGTTGTTGCTGAAGCCGCCTGGGGTGCTCAGTTGCTTGACGGCCAGTTGCGCACCTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGAcagcagcag                                                                                                                                  >  1:284553/1‑150 (MQ=255)
                    cGGGTTGTTGCTGAAGCCGCCTGGGGTGCTCAGTTGCTTGACGGCCAGTTGCGCACCTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGcc                                                                                                                               <  2:284553/150‑1 (MQ=255)
                      ggTTGTTGCTGAAGCCGCCTGGGGTGCTCAGTTGCTTGACGGCCAGTTGCGCACCTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCgg                                                                                                                             >  2:385635/1‑150 (MQ=255)
                              cTGAAGCCGCCTGGGGTGCTCAGTTGCTTGACGGCCAGTTGCGCACCTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAgg                                                                                                                     <  1:191065/150‑1 (MQ=255)
                                   gccgccTGGGGTGCTCAGTTGCTTGACGGCCAGTTGCGCACCTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGTGCCCTGGccct                                                                                                                                           <  2:1275110/123‑3 (MQ=255)
                                   gccgccTGGGGTGCTCAGTTGCTTGACGGCCAGTTGCGCACCTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGTGCCCTGGccct                                                                                                                                           >  1:1275110/1‑121 (MQ=255)
                                   gccgccTGGGGTGCTCAGTTGCTTGACGGCCAGTTGCGCACCTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGggtggt                                                                                                                >  2:392236/1‑150 (MQ=255)
                                      gccTGGGGTGCTCAGTTGCTTGACGGCCAGTTGCGCACCTCGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAAGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGCCAGCAGCAGGCCGGCGCACAGGGTGGGGaa                                                                                                             >  2:509634/1‑150 (MQ=255)
                                        cTGGGGTGCTCAGTTGCTTGACGGCCAGTTGCGCACCTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCACCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAgc                                                                                                           <  2:460925/150‑1 (MQ=255)
                                           gggTGCTCAGTTGCTTGACGGCCAGTTGCGCACCTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCa                                                                                                        <  1:1044173/150‑1 (MQ=255)
                                            ggTGCTCAGTTGCTTGACGGCCAGTTGCGCACCTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAg                                                                                                       <  2:463173/150‑1 (MQ=255)
                                                cTCAGTTGCTTGACGGCCAGTTGCGCACCTTGGGTCAGGGNATCTTTNAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGaa                                                                                                   >  1:1093273/1‑150 (MQ=255)
                                                 tCAGTTGCTTGACGGCCAGTTGCGCACCTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAg                                                                                                  >  1:1356860/1‑150 (MQ=255)
                                                  cAGTTGCTTGACGGCCAGTTGCGCACCTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGt                                                                                                 <  2:1356860/150‑1 (MQ=255)
                                                        cTTGACGGCCAGTTGCGCACCTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGAt                                                                                           <  1:925226/150‑1 (MQ=255)
                                                          tGACGGCCAGTTGCGCACCTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCa                                                                                         >  2:910732/1‑150 (MQ=255)
                                                          ggACGGCCAGTTGCGCACCTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACANCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCa                                                                                         <  1:509634/149‑1 (MQ=255)
                                                               gCCAGTTGCGCACCTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCATGAAg                                                                                    >  1:721197/1‑150 (MQ=255)
                                                                      gcgcACCTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCt                                                                             >  2:1351444/1‑150 (MQ=255)
                                                                       cgcACCTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCtt                                                                            <  1:1349173/150‑1 (MQ=255)
                                                                           ccTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCTTGTaa                                                                        <  2:422890/150‑1 (MQ=255)
                                                                           ccTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCTTGTaa                                                                        >  1:935429/1‑150 (MQ=255)
                                                                            cTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGaa                                                                                                             <  1:687427/112‑1 (MQ=255)
                                                                            cTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGaa                                                                                                             >  2:687427/1‑112 (MQ=255)
                                                                            cTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGaa                                                                                                   <  1:458807/122‑1 (MQ=255)
                                                                            cTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGaa                                                                                                   >  2:458807/1‑122 (MQ=255)
                                                                                 gTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCTTGTAAACGGAt                                                                  >  1:468854/1‑150 (MQ=255)
                                                                                     gggCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCTTGTAAACGGATGAGt                                                              >  2:1008004/1‑150 (MQ=255)
                                                                                         aTCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCTTGTAAACGGATGAGTACgg                                                          <  1:1008004/150‑1 (MQ=255)
                                                                                         aTCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCTTGTAAACGGATGAGTACgg                                                          <  2:1102993/150‑1 (MQ=255)
                                                                                         aTCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCTTGTAAACGGATGAGTACgg                                                          <  2:1273039/150‑1 (MQ=255)
                                                                                            tttCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCTTGTAAACGGATGAGTACGGGTa                                                       >  2:1071951/1‑150 (MQ=255)
                                                                                                aGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCTTGTAAACGGATGAGTACGGGTATCAg                                                   <  1:910732/150‑1 (MQ=255)
                                                                                                 ggCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCTTGTAAACGGATGAGTACGGGTATCAGc                                                  <  1:385635/150‑1 (MQ=255)
                                                                                                        ggTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCCTGAAGCTTTCCTTGTAAACGGATGAGTACGGGTATCAGCGAACATg                                           <  1:392236/150‑1 (MQ=255)
                                                                                                        ggTGGCATCATTCGGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCTTGTAAACGGATGAGTACGGGTATCAGCGAACCTg                                           <  2:1105500/150‑1 (MQ=255)
                                                                                                            gCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGggtggt                                                                                                                <  1:516410/77‑1 (MQ=255)
                                                                                                            gCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGggtggt                                                                                                                >  2:516410/1‑77 (MQ=255)
                                                                                                                 tttCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCTTGTAAACGGATGAGTACGGGTATCAGCGAACCTGGTCTACCtt                                  >  2:560632/1‑150 (MQ=255)
                                                                                                                   tCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGAGGTGAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCTTGTAAACGGATGAGTACGGGTATCAGCGAACCTGGTCTACCTTGa                                >  1:406116/1‑150 (MQ=255)
                                                                                                                    cTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGACGGCCAGGGAAGTGCGGATCATGAAGCTTTCCTTGTAAACGGATGAGTACGGGTATCAGCGAACCTGGTCTACCTTGAt                               <  1:93628/150‑1 (MQ=255)
                                                                                                                     tGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCTTGTAAACGGATGAGTACGGGTATCAGCGAACCTGGTCTACCTTGATg                              >  2:512664/1‑150 (MQ=255)
                                                                                                                       gCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCTTGTAAAc                                                                      >  2:597062/1‑108 (MQ=255)
                                                                                                                       gCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCTTGTAAAc                                                                      <  1:597062/108‑1 (MQ=255)
                                                                                                                         tGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGAt                                                                                           >  2:1024644/1‑85 (MQ=255)
                                                                                                                         tGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGAt                                                                                           <  1:1024644/85‑1 (MQ=255)
                                                                                                                            ggTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCTTGTAAACGGATGAGTACGGGTATCAGCGAACCTGGTCTACCTTGATGCgagatc                       >  1:155211/1‑145 (MQ=255)
                                                                                                                                gCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCTTGTAAACGGATGAGTACGGGTATCAGCGAACCTGGTCTACCTTGATGCGTACCGATTg                   >  2:365722/1‑150 (MQ=255)
                                                                                                                                  ccAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCTTGTAAACGGATGAGTACGGGTATCAGCGAACCTGGTCTACCTTGATGCGTACCGATTGCg                 >  1:589592/1‑150 (MQ=255)
                                                                                                                                        aCAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCAGAGGGTGGTGAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCTTGTAAACGGATGAGTACGGGTATCAGCGAACCTGGTCTACCTTGATGCGTACCGATTGCGGGTCAc           >  1:1362328/1‑150 (MQ=255)
                                                                                                                                         cagggccagggcACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCTTGTAAACGGATGAGTACGGGTATCAGCGAACCTGGTCTACCTTGATGCGTACCGATTGCGGGTCACt          <  2:589592/150‑1 (MQ=255)
                                                                                                                                           gggccagggcACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCTTGTAAACGGATGAGTACGGGTATCAGCGAACCTGGTc                                        <  1:223999/118‑1 (MQ=255)
                                                                                                                                           gggccagggcACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCTTGTAAACGGATGAGTACGGGTATCAGCGAACCTGGTc                                        >  2:223999/1‑118 (MQ=255)
                                                                                                                                            aaccagtgcCCTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCTTGTAAACGGATGAGTACGGGTATCAGCGAACCTGGTCTACCTTGATGCGTACCGATTGCGGGTCACTGcc       >  1:1119194/3‑150 (MQ=255)
                                                                                                                                               cagggcACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATc                                                                                          <  1:1003740/64‑1 (MQ=255)
                                                                                                                                               cagggcACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATc                                                                                          >  2:1003740/1‑64 (MQ=255)
                                                                                                                                               cagggcACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCTTGTAAACGGATGAGTACGGGTATCAGCGAACCTGGTCTACCTTGATGCGTACCGATTGCGGGTCACTGCCATc    <  2:58494/150‑1 (MQ=255)
                                                                                                                                               cagggcACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTCAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCTTGTAAACGGATGAGTACGGGTATCAGCGAACCTGGTCTACCTTGATGCGTACCGATTGCGGGTCACTGCCATc    <  1:512664/150‑1 (MQ=255)
                                                                                                                                                agggcACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCTTGTAAACGGATGAGTACGGGTATCAGCGAACCTGGTCTACCTTGATGCGTACCGATTGCGGGTCACTGCCATCg   >  1:1060703/1‑150 (MQ=255)
                                                                                                                                                 gggcACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCTTGTAAACGGATGAGTACGGGTATCAGCGAACCTGGTCTACCTTGATGCGTACCGATTGCGGGTCACTGCCATCGa  <  1:365722/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GGCAGCTCGATGCGCACGTCCGGGTTGTTGCTGAAGCCGCCTGGGGTGCTCAGTTGCTTGACGGCCAGTTGCGCACCTTGGGTCAGGGCATCTTTCAGGCCGCCGGTGGCATCTTTCTGGCTGAGGTCGCCCAGCGACAGGGCCAGGGCACTGGCCGACAGCAGCAGGCCGGCGCACAGGGTGGTGAAGCGCAGGGAAGTGCGGATCATGAAGCTTTCCTTGTAAACGGATGAGTACGGGTATCAGCGAACCTGGTCTACCTTGATGCGTACCGATTGCGGGTCACTGCCATCGA  >  NC_002947/1725002‑1725296

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