Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,436,556 C→T 6.2% A170T (GCC→ACC)  PP_1257 ← 1‑pyrroline‑4‑hydroxy‑2‑carboxylate deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,436,5560CT 6.2% 191.2 / 2.0 64A170T (GCC→ACC) PP_12571‑pyrroline‑4‑hydroxy‑2‑carboxylate deaminase
Reads supporting (aligned to +/- strand):  ref base C (30/30);  new base T (2/2);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

CGCGGCAGCGCATTGCCGATGCCGGTGATGGCGCCGGTGGCGTTGCAGTTGACGAAGCCATGCACCACCTGGGTGTCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCGTAGTAAGGGCTGTTGTAGATCACCGCCGGCAGTTTCGGTGCGGCAGCGAGGATGG  >  NC_002947/1436411‑1436696
                                                                                                                                                 |                                                                                                                                            
cgcgGCAGCGCATTGCCGATGCCGGTGATGGCGCCGGTGGCGTTGCAGTTGACGAAGCCATGCACCACCTGGGTGTCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGcc                                                                                                                                           <  1:151668/149‑1 (MQ=255)
        cgcATTGCCGATGCCGGTGATGGCGCCGGTGGCGTTGCAGTTGACGAAGCCATGCACCACCTGGGTGTCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCaa                                                                                                                                   <  4:46641/149‑1 (MQ=255)
         gcATTGCCGATGCCGGTGATGGCGCCGGTGGCGTTGCAGTTGACGAAGCCATGCACCACCTGGGTGTCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCaaa                                                                                                                                  >  8:297823/1‑149 (MQ=255)
          cATTGCCGATGCCGGTGATGGCGCCGGTGGCGTTGCAGTTGACGAAGCCATGCACCACCTGGGTGTCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGTGCCACCGCCAAAc                                                                                                                                 >  4:241438/1‑149 (MQ=255)
          cATTGCCGATGCCGGTGATGGCGCCGGTGGCGTTGCAGTTGACGAAGCCATGCACCACCTGGGTGTCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGTGCCACCGCCAAAc                                                                                                                                 >  4:245796/1‑149 (MQ=255)
              gCCGATGCCGGTGATGGCGCCGGTGGCGTTGCAGTTGACGAAGCCATGCACCACCTGGGTGTCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTcc                                                                                                                              <  8:127838/148‑1 (MQ=255)
               ccGATGCCGGTGATGGCGCCGGTGGCGTTGCAGTTGACGAAGCCATGCACCACCTGGGTGTCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCtt                                                                                                                            <  8:134021/149‑1 (MQ=255)
                cGATGCCGGTGATGGCGCCGGTGGCGTTGCAGTTGACGAAGCCATGCACCACCTGGGTGTCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCtt                                                                                                                            >  2:160853/1‑148 (MQ=255)
                 gATGCCGGTGATGGCGCCGGTGGCGTTGCAGTTGACGAAGCCATGCACCACCTTGGTGTCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGa                                                                                                                          <  8:104965/149‑1 (MQ=255)
                 gATGCCGGTGATGGCGCCGGTGGCGTTGCAGTTGACGAAGCCATGCACCACCTGGGTGTCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTg                                                                                                                           >  8:119976/1‑148 (MQ=255)
                      cGGTGATGGCGCCGGTGGCGTTGCAGTTGACGAAGCCATGCACCACCTGGGTGTCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGcc                                                                                                                      >  1:12231/1‑148 (MQ=255)
                         tGATGGCGCCGGTGGCGTTGCAGTTGACGAAGCCATGCACCACCTGGGTGTCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATc                                                                                                                  <  1:160853/149‑1 (MQ=255)
                                 ccGGTGGCGTTGCAGTTGACGAAGCCATGCACCACCTGGGTGTCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTggg                                                                                                          <  4:276953/149‑1 (MQ=255)
                                        cGTTGCAGTTGACGAAGCCATGCACCACCTGGGTGTCGAAACCGAGCATGAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTGGGCTCCACCGCCAAACTGCTTGAAGCGGTTCAGGTTGGGGAACTgg                                                                                                   <  2:135611/149‑1 (MQ=255)
                                          ttGCAGTTGACGAAGCCATGCACCACCTGGGTGTCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCt                                                                                                 >  2:173899/1‑149 (MQ=255)
                                               gTTGACGAAGCCATGCACCACCTGGGTGTCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGgcgc                                                                                                                                            >  6:166691/1‑101 (MQ=255)
                                               gTTGACGAAGCCATGCACCACCTGGGTGTCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGgcgc                                                                                                                                            <  5:166691/101‑1 (MQ=255)
                                                    cGAAGCCATGCACCACCTGGGTGTCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTc                                                                                        >  8:27899/1‑148 (MQ=255)
                                                           aTGCACCACCTGGGTGTCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCgaagaa                                                                                  >  4:223760/1‑147 (MQ=255)
                                                                ccaccTGGGTGTCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTggg                                                                                                          <  4:95742/118‑1 (MQ=255)
                                                                ccaccTGGGTGTCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTggg                                                                                                          >  3:95742/1‑118 (MQ=255)
                                                                      gggTGTCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGgcgc                                                                       >  4:305372/1‑147 (MQ=255)
                                                                       ggTGTCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGa                                                                     <  1:165931/148‑1 (MQ=255)
                                                                       ggTGTCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGa                                                                     <  1:173899/148‑1 (MQ=255)
                                                                         tgtCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTg                                                                  >  6:33471/1‑149 (MQ=255)
                                                                         ggtCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTg                                                                  <  8:190730/148‑1 (MQ=255)
                                                                          gtCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGTGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTgg                                                                 <  3:245796/149‑1 (MQ=255)
                                                                          gtCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTgcg                                                                                              <  1:167497/120‑1 (MQ=255)
                                                                          gtCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTgcg                                                                                              >  2:167497/1‑120 (MQ=255)
                                                                           tCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGTGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTgg                                                                 <  3:241438/148‑1 (MQ=255)
                                                                           tCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGc                                                                <  7:27899/149‑1 (MQ=255)
                                                                                acCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGc                                                           <  7:297823/149‑1 (MQ=255)
                                                                                  cGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCg                                                         >  6:284891/1‑149 (MQ=255)
                                                                                   gACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCgt                                                        <  7:119976/149‑1 (MQ=255)
                                                                                   gACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCgt                                                        <  5:284891/149‑1 (MQ=255)
                                                                                   gACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTACTTGAAGCCGATCAGGATGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCg                                                         >  3:90026/1‑148 (MQ=255)
                                                                                    aCCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCgta                                                       <  5:188137/149‑1 (MQ=255)
                                                                                     ccATCAGTGTGACTTCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCgta                                                       <  8:220130/148‑1 (MQ=255)
                                                                                      cATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCgtag                                                      >  3:108902/1‑148 (MQ=255)
                                                                                          aGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCa                                                                                            >  4:144421/1‑106 (MQ=255)
                                                                                          aGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCa                                                                                            <  3:144421/106‑1 (MQ=255)
                                                                                           gggTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCGTAGTAAggg                                                >  4:263271/1‑149 (MQ=255)
                                                                                            ggTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCGTAGTAAGGGc                                               >  3:38508/1‑149 (MQ=255)
                                                                                                  aTCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGCTCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCgta                                                       >  3:175284/1‑135 (MQ=255)
                                                                                                  aTCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCgtag                                                      >  2:269358/1‑136 (MQ=255)
                                                                                                  aTCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCgtag                                                      <  1:269358/136‑1 (MQ=255)
                                                                                                  aTCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCgta                                                       <  4:175284/135‑1 (MQ=255)
                                                                                                  aTCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCGTAGTAAGGGCtgttgt                                         <  5:127974/149‑1 (MQ=255)
                                                                                                  aTCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCGTAGTAAGGGCtgttg                                          >  6:127974/1‑148 (MQ=255)
                                                                                                  aTCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCGTAGTAAGGGCtgttg                                          >  6:98537/1‑148 (MQ=255)
                                                                                                  aTCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCGTAGTAAGGGCtgttg                                          >  7:165954/1‑148 (MQ=255)
                                                                                                   tCGTCATCCTTGGAGGTGATGTGCTCGGCGGAGTAGCGCAAGTCGGCGCAACCGCCAAACTCCTTGAAGCCGATCAGGTTGGTGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCGTAGTAAGGGCtgttg                                          >  5:279914/1‑147 (MQ=255)
                                                                                                               gAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCGTAGTAAGGGCTGTTGTAGATCACCGCCgg                            <  5:33471/149‑1 (MQ=255)
                                                                                                                   tGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCGTAGTAAGGGCTGTTGTAGATCACCGCCGGCAGt                        >  2:24862/1‑149 (MQ=255)
                                                                                                                      tgtgCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCGTAGTAAGGGCTGTTGTAGATCACCGCCGGCAGTTTc                     <  3:223760/149‑1 (MQ=255)
                                                                                                                       gtgCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCGTAGTAAGGGCTGTTGTAGATCACCGCCGGCAGTTTCg                    >  5:259691/1‑149 (MQ=255)
                                                                                                                       gtgCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCGTAGTAAGGGCTGTTGTAGATCACCGCCGGCAGTTTCg                    >  2:231343/1‑149 (MQ=255)
                                                                                                                        tgCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCGTAGTAAGGGCTGTTGTAGATCACCGCCGGCAGTTTCg                    <  3:263271/148‑1 (MQ=255)
                                                                                                                         gCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCGTAGTAAGGGCTGTTGTAGATCACCGCCGGCAGTTTCGGt                  >  7:219182/1‑149 (MQ=255)
                                                                                                                         gCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCGTAGTAAGGGCTGTTGTAGATCACCGCCGGCAGTTTCGGt                  >  3:157570/1‑149 (MQ=255)
                                                                                                                          cTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCGTAGTAAGGGCTGTTGTAGATCACCGCCGGCAGTTTCGGTg                 >  8:137436/1‑149 (MQ=255)
                                                                                                                              gcggcgTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCGTAGTAAGGGCTGTTGTAGATCACCGCCGGCAGTTTCGGTGCGGc             <  4:157570/149‑1 (MQ=255)
                                                                                                                                      gcgcAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCGTAGTAAGGGCTGTTGTAGATCACCGCCGGCAGTTTCGGTGCGGCAGCGAGGa     >  3:207080/1‑149 (MQ=255)
                                                                                                                                        gcAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCGTAGTAAGGGCTGTTGTAGATCACCGCCGGCAGTTTCGGTGCGGCAGCGAGGATg   <  6:259691/149‑1 (MQ=255)
                                                                                                                                         cAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCGTAGTAAGGGCTGTTGTAGATCACCGCCGGCAGTTTCGGTGCGGCAGCGAGGATgg  <  4:207080/149‑1 (MQ=255)
                                                                                                                                         cAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCGTAGTAAGGGCTGTTGTAGATCACCGCCGGCAGTTTCGGTGCGGCAGCGAGGATgg  <  8:219182/149‑1 (MQ=255)
                                                                                                                                                 cgcCACGTCAAAACTCCTTGATGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCGTAGTAAGGGCTGTTGTAGATCACCGCCGGCAGTTTCGGTGCGGCAGCGAGGa     <  6:228139/138‑1 (MQ=255)
                                                                                                                                                 cgcCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTAGAAGAATAGGTCGGCGCGAGTGGCAAAGCCGTAGTAAGGGCTGTTGTAGATCACCGCCGGCAGTTTCGGTGCGGCAGCGAGGa     >  5:228139/1‑138 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
CGCGGCAGCGCATTGCCGATGCCGGTGATGGCGCCGGTGGCGTTGCAGTTGACGAAGCCATGCACCACCTGGGTGTCGACACCGACCATCAGGGTGACATCGTCATCCTTGGAGGTGATGTGCTCGGCGGCGTAGCGCAAGTCGGCGCCACCGCCAAACTCCTTGAAGCCGATCAGGTTGGGGAACTGGCTGCGCAGTTCGAAGAATAGGTCGGCGCGAGTGGCAAAGCCGTAGTAAGGGCTGTTGTAGATCACCGCCGGCAGTTTCGGTGCGGCAGCGAGGATGG  >  NC_002947/1436411‑1436696

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