Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1435,530,4210AC31.3% 97.1 / 10.7 67intergenic (‑53/+102)livH_5/braC_2High‑affinity branched‑chain amino acidtransport system permease protein LivH/Leucine‑, isoleucine‑, valine‑, threonine‑, and alanine‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (19/27);  new base C (21/0);  total (40/27)
Fisher's exact test for biased strand distribution p-value = 1.38e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTGTAGCCGATGGCGATCAGGCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTACCGCGGGGCTTGGGCCCGAGCAACCGCCCGGGCCGCGGGCCGTTGTTATTTTTGTTTTTCCAGCTGGTGGTATT  >  ACB143/5530289‑5530551
                                                                                                                                     |                                                                                                                                  
gtgtAGCCGATGGCGATCAGGCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCtg                                                                                                                                <  1:36445/137‑1 (MQ=255)
 tgtAGCCGATGGCGATCAGGCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCCCCtgt                                                                                                                               >  5:40268/1‑137 (MQ=255)
 tgtAGCCGATGGCGATCAGGCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCtgt                                                                                                                               >  1:313218/1‑137 (MQ=255)
 tgtAGCCGATGGCGATCAGGCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACAATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCtgt                                                                                                                               >  4:16659/1‑137 (MQ=255)
  gtAGCCGATGGCGATCAGGCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGGAAAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCCACtgt                                                                                                                               >  4:326637/1‑137 (MQ=255)
     gCCGATGGCGATCAGGCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCTGTGCgg                                                                                                                           >  1:377308/1‑137 (MQ=255)
      ccGATGGCGATCAGGCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCCCCTGTGgggg                                                                                                                          >  1:209294/1‑137 (MQ=255)
         aTGGCGATCAGGCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCTGTGCGGGGcc                                                                                                                       >  8:89675/1‑137 (MQ=255)
          tGGCGATCAGGCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGT‑AAAATACCATCCATACCGCAATCTCACCTGATTGAGATTGCACGAGCCCCCACCACGGCGCGGGCCACCTGTGCGGGGccc                                                                                                                      <  4:419024/137‑1 (MQ=255)
          tGGCGATCAGGCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCTGTGCGGGGccc                                                                                                                      <  4:153279/137‑1 (MQ=255)
          tGGCGATCAGGCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCTGTGCGGGGccc                                                                                                                      <  7:127238/137‑1 (MQ=255)
          gggCGATCCGGCCATAGACCGCCCCGCGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATACCGCAAGCTCCCCTGATTGAGATTGCACGAGCCCCCCCCCCGGCGCGGGCCACCTGTGCGGGGccc                                                                                                                      <  2:160880/136‑1 (MQ=255)
          gggCGATCAGCCCATAGCCCCACCCGCGGGACCGCCCGTTCACCATTTCCGCCAGC‑AAAAAACCACCCAAACCCCAATCGCACCTGATTCAGATTGCACGCGCGCCCACCACGGCGCGGGCCACCTGTGCGGGGccc                                                                                                                      <  2:456129/136‑1 (MQ=255)
            gCGATCAGGCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCtgtg                                                                                                                              >  1:58667/1‑127 (MQ=255)
            gCGATCAGGCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCtgtg                                                                                                                              <  2:58667/127‑1 (MQ=255)
            gCGATCAGGCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCTGTGCGGGGCCCTc                                                                                                                    >  7:232190/1‑137 (MQ=255)
            gCGATCAGGCCATAGACCGACCAGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCCCCTGTGCGGGGCCCTc                                                                                                                    >  4:28190/1‑137 (MQ=255)
              gATCAGGCCATAGCACGACCCGAGGGTCAGGCCGTTGCCCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTCAGATTGCACGAGCGCCGACCACGGCGCGGGGCACCTGGGGGGGGGCCgcgc                                                                                                                  >  3:75312/1‑137 (MQ=255)
              gATCAGGCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCCCCTGTGCGGGGCCCTcgc                                                                                                                  >  7:347982/1‑137 (MQ=255)
              gATCAGGCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCTGTGCGGGGCCCTcgc                                                                                                                  >  6:241806/1‑137 (MQ=255)
               aTCCGGCCATACCCCCACCCTCGGGTCAGGCCGTTCACCAGTTGCTCCAGA‑AAAAAACCACCCATAACGCAATCTCACCTGATTGATATTGCACGAGCGCCCCCCCCGGCGCGGGCCACCTGTGCGGGGCCCTcgca                                                                                                                 <  2:446258/137‑1 (MQ=255)
                    gCCATAGACCGCCCCGCGGGGCAGGCCGTTGACCAGTTGCGGCAGG‑AAAATCCCAACCAAACCGCAATCGACCCTGATTGAGATTGCACGAGCGCCGACCACGCCGCGGGCCACCTGTGCGGGGCCCTcgcagcgca                                                                                                            <  8:273671/137‑1 (MQ=255)
                    gCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCTGTGCGGGGCCCTcgcagcgca                                                                                                            <  5:469561/137‑1 (MQ=255)
                    gCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCTGTGCGGGGCCCTcgcagcgca                                                                                                            <  6:468830/137‑1 (MQ=255)
                                cccGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGa                                                                                                >  5:377126/1‑137 (MQ=255)
                                cccGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCTGTGCGGGGCCCGCGCAGCGCAGAACTGTGCAGa                                                                                                >  2:274388/1‑137 (MQ=255)
                                 ccGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCCCCGGTGCGGGGCCCTAGCAGCGCAGAACTGTGccgct                                                                                               >  7:301453/1‑133 (MQ=255)
                                 ccGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCCCCGGGGCGGGGCCCCCGCCGCGCAGAAATGTGCGGag                                                                                               >  5:289171/1‑136 (MQ=255)
                                 ccGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCCCCAGTGGGGGGGCCTCGCAGCGCAGCACTGTGCAGAt                                                                                               >  5:54782/1‑137 (MQ=255)
                                 ccGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCTGTGCGGGGCCCGCGCAGCGCAGAACTGTGCAGAt                                                                                               >  1:316609/1‑137 (MQ=255)
                                 ccGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCGGTGCGGGGCCCTCGCAGCGCAGAACTGTGGAGAt                                                                                               >  4:315625/1‑137 (MQ=255)
                                 ccGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCGGGGCGGGGCCCTCGCCGCGCAGAACAGTGCAGAt                                                                                               >  8:169327/1‑137 (MQ=255)
                                  cGAGGTTCAGGCCGTTTACCAGTTTCTGAAGG‑AAAATACCATCCATAAAGCAATCTCACCTTATTGATATTGCACGAGCGACGACAACGGCGCGGGCACCCTGTGGTGGGACCTCGCAGCGCAGAACTGTGCAgctt                                                                                              >  5:336071/1‑134 (MQ=255)
                                  cGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGCTTGCACGAGCGCCGCCCACGGCGCGGGGCCCCTGTGCGGGGCCCTCGCAGCGCAGAACGGTGCAgcta                                                                                              >  2:327240/1‑134 (MQ=255)
                                  cGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCCCCTGTGCGGGGCCCTCGCAGCGCCGAACTGTGCAGAta                                                                                              >  7:101801/1‑136 (MQ=255)
                                  cGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCCCCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCCGAtg                                                                                              >  5:80687/1‑136 (MQ=255)
                                  cGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCCCCTGTGCGGGGCCCGCTCAGCGCAGAACTGTGGAGATc                                                                                              >  1:43472/1‑137 (MQ=255)
                                  cGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCCCCTGTGCGGGGCCCGCGCAGCGCAGAACTGTGCAGATc                                                                                              >  7:39785/1‑137 (MQ=255)
                                  cGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCTGTGCGGGGCCCTCGCAGCGCAGCACTGGGCAGagg                                                                                              >  5:343307/1‑135 (MQ=255)
                                   gAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCCCCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCt                                                                                             >  4:365853/1‑137 (MQ=255)
                                       gTCAGGCCGTTGACCTGTTGCTGAAGT‑AAAATACCATCCATAACGCAATCTCACGTTATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGa                                                                                         <  8:94278/137‑1 (MQ=255)
                                            gCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATaaa                                                                                    <  2:377308/137‑1 (MQ=255)
                                             ccGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCCCCTGTGAGGGGCCCTCGCCGCGCAGAAATCTGCAGCTCTAAGCATAAAg                                                                                   >  4:274018/1‑137 (MQ=255)
                                                    ccAGTTGCTGCAGG‑AAAATACCATCCATAAGGCAATCTCCCCTGATTGTGATTGCACGACCGCCGACCACGGCGCGGGCCACCTTTGCGGGGCCCGCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTa                                                                            <  4:75312/137‑1 (MQ=255)
                                                    ccAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCTGTGCGGGGCTCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTa                                                                            >  2:290539/1‑137 (MQ=255)
                                                      atttGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTAcc                                                                          <  6:442651/135‑1 (MQ=255)
                                                           ctgcAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCCCCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTACCGCggg                                                                     >  3:74261/1‑137 (MQ=255)
                                                           ctgcAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCCCCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTACCGCggg                                                                     >  7:261284/1‑137 (MQ=255)
                                                                gg‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTACCGCGGGGCTTg                                                                <  5:391252/137‑1 (MQ=255)
                                                                gg‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCTGTGCGGGGCCCGCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTACCGCGGGGCTTg                                                                >  7:31988/1‑137 (MQ=255)
                                                                    aaaTACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCCCCTGTGCGGGGCCCGCGCAGCGCCGAACTGTGCAGATCGACGAATAAAGACAGGTACCGCGGGGCTTGGGc                                                             >  6:342228/1‑137 (MQ=255)
                                                                        aCCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGCCCACGGCGCGGGCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTACCGCGGGGCTTGGGCCCGa                                                         <  7:423515/137‑1 (MQ=255)
                                                                              aaaaaCCCAAACTCCCATGATTGAGATGGCACGAGCCCCGACCCGGGCGCGGGCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTACCGCGGGGCTTGGGCCCGAGCAAcc                                                   <  6:40268/134‑1 (MQ=255)
                                                                                         ctcACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTACCGCGGGGCTTGGGCCCGAGCAACCGCCcgggccgc                                        >  4:470807/1‑137 (MQ=255)
                                                                                         ctcACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTACCGCGGGGCTTGGGCCCCAGCAACCGCCcgggccgc                                        >  4:307427/1‑137 (MQ=255)
                                                                                           cACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCCCCAGTGCCGGGCCCGCGCAGCGCAGAACTGTGGAGATCGAACAATAAAGACAGGTAACGCGGGGCTTGGGCCCCAGCAACCCCCccggccgcgg                                      >  8:22967/1‑137 (MQ=255)
                                                                                                gATTTACATTGCCCGCGCCCCCCCCCCGGCGCGGGCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTACCGCGGGGCTTGGGCCCGAGCAACCGCCCGGGCCGCGGGCCgt                                 <  1:327240/137‑1 (MQ=255)
                                                                                                gATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCCCCTGTGCGGGGCCCTCGCAGCGCAGAACAGTGCAGATCTACGAATAAAGACAGGTACCGCGGGGCTTGGGCCCGAGCAACCCCCccggccgcggggccg                                 >  3:395235/1‑133 (MQ=255)
                                                                                                   tGAGATTGCACGAGCGCCGACCACGGCGCGGGGCACCTGTGCGGGGCCCAAGCAGCGCAGAACAGTGCAGATCTACGAAAAAAGACCCGTACCGCGGGGCTTGGGCCCGCGCAACCCCCccggccgccgggcggggg                              >  6:44423/1‑133 (MQ=255)
                                                                                                               cgcgcCGCCCACGGCGCGGGCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTACCGCGGGGCTTGGGCCCGAGCAACCGCCCGGGCCGCGGGCCGTTGTTATTTTTGtttt                  <  6:51975/136‑1 (MQ=255)
                                                                                                                 cgcCGACCACGGCGCGGGCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTACCGCGGGGCTTGGGCCCGAGCAACCGCCCGGGCCGCGGGCCGTTGTTATTTTTGTTTTTc                <  8:124376/137‑1 (MQ=255)
                                                                                                                  gcCGACCACGGCGCGGGCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTACCGCGGGGCTTGGGCCCGAGCAACCGCCCGGGCCGCGGGCCGTTGTTATTTTTGTTTTTcc               <  2:209294/137‑1 (MQ=255)
                                                                                                                    cGACCCGGGCGCGGGCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTACCGCGGGGCTTGGGCCCGAGCAACCGCCCGGGCCGCGGGCCGTTGTTATTTTTGTTTTTCCAg             <  3:16659/137‑1 (MQ=255)
                                                                                                                    cGACCACGGCGCGGGCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTACCGCGGGGCTTGGGCCCGAGCAACCGCCCGGGCCGCGGGCCGTTGTTATTTTTGTTTTTCCAg             <  4:457783/137‑1 (MQ=255)
                                                                                                                     gACCACGGGGCGGGCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTACCGCGGGGCTTGGGCCCGAGCAACCGCCCGGGCCGCGGGCCGTTGTTATTTTTGTTTTTCCAGc            <  3:219121/137‑1 (MQ=255)
                                                                                                                          gggCGCGGGCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTACCGCGGGGCTTGGGCCCGAGCAACCGCCCGGGCCGCGGGCCGTTGTTATTTTTGTTTTTCCAGCtggtg       <  4:395235/136‑1 (MQ=255)
                                                                                                                               cgGGCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTACCGCGGGGCTTGGGCCCGAGCAACCGCCCGGGCCGCGGGCCGTTGTTATTTTTGTTTTTCCAGCTGGTGGTAtt  <  1:188144/137‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                  
GTGTAGCCGATGGCGATCAGGCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGG‑AAAATACCATCCATAACGCAATCTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTACCGCGGGGCTTGGGCCCGAGCAACCGCCCGGGCCGCGGGCCGTTGTTATTTTTGTTTTTCCAGCTGGTGGTATT  >  ACB143/5530289‑5530551

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