Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,344,346 C→G 10.2% Q261H (CAG→CAC)  PP_2941 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,344,3460CG10.2% 196.7 / 10.5 67Q261H (CAG→CAC) PP_2941hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (26/34);  new base G (2/5);  total (28/39)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

CATTCTCGAAGGGAAGGTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATCCA  >  NC_002947/3344199‑3344492
                                                                                                                                                   |                                                                                                                                                  
cATTCTCGAAGGGAAGGTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACt                                                                                                                                                   >  5:348620/1‑149 (MQ=255)
    ctcGAAGGGAAGGTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCAGTGAAg                                                                                                                                               >  7:212217/1‑149 (MQ=255)
       gaagggaaggTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCgg                                                                                                                                            <  6:348620/149‑1 (MQ=255)
       gaagggaaggTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCgg                                                                                                                                            <  3:337419/149‑1 (MQ=255)
        aagggaaggTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCAGTGAAGCGGc                                                                                                                                           <  4:129503/149‑1 (MQ=255)
          gggaaggTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCAGTGAAGCGGctg                                                                                                                                         <  3:202676/149‑2 (MQ=255)
             aaggGTAGGGCTGAACAGTACGTGGAAAGGCAGAAACGCGGAGCTTGTGCTGCCAGGTGAGTGCGCCGAGGTCTGAGCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTAt                                                                                                                                      <  6:240105/149‑1 (MQ=255)
                gTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCAGTGAAGCGGctgcacaga                                                                                                                                   <  8:212217/149‑8 (MQ=255)
                gTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTAtcag                                                                                                                                   >  4:72605/1‑149 (MQ=255)
                gTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTAtcag                                                                                                                                   <  3:234736/149‑1 (MQ=255)
                           aCAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCtt                                                                                                                        <  3:55942/149‑1 (MQ=255)
                              gTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGAt                                                                                                                     >  7:262169/1‑149 (MQ=255)
                               tACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGAtt                                                                                                                    <  7:135637/149‑1 (MQ=255)
                                 cGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTg                                                                                                                   <  2:121923/148‑1 (MQ=255)
                                         gTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTg                                                                                                          <  1:94739/149‑1 (MQ=255)
                                          tCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTgg                                                                                                         <  8:295698/149‑1 (MQ=255)
                                          tCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTgg                                                                                                         <  2:203305/149‑1 (MQ=255)
                                           cAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTgg                                                                                                         <  7:350514/148‑1 (MQ=255)
                                            aGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTAtc                                                                                                                                     >  4:37546/1‑119 (MQ=255)
                                            aGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTAtc                                                                                                                                     <  3:37546/119‑1 (MQ=255)
                                            aGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGc                                                                                                        >  6:360744/1‑148 (MQ=255)
                                             gAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCaa                                                                                                      >  8:160990/1‑149 (MQ=255)
                                              aaaaGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAAt                                                                                                     >  7:15654/1‑149 (MQ=255)
                                                aaGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATgg                                                                                                   <  7:160990/149‑1 (MQ=255)
                                                    ggAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCtgt                                                                                               <  6:18050/149‑1 (MQ=255)
                                                       gCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCt                                                                                            <  8:234273/149‑1 (MQ=255)
                                                           gtgCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACt                                                                                        >  1:230963/1‑149 (MQ=255)
                                                              ctgcCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGa                                                                                     <  8:15654/149‑1 (MQ=255)
                                                                 ccAGGTGAGTTCGCCGCGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATAc                                                                                  >  4:91409/1‑149 (MQ=255)
                                                                 ccAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATa                                                                                   <  2:126996/148‑1 (MQ=255)
                                                                 ccAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATAc                                                                                  <  5:360744/149‑1 (MQ=255)
                                                                 ccAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAAACACGGGCAGGTCGACAACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCAATCAGCTCAGATCCTAGATTGCGAACATAGGCAATGGCTGTGCTCACTCGATAc                                                                                  >  4:298597/1‑149 (MQ=255)
                                                                 ccAGGTGAGTTCGCCGAGGCCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATa                                                                                   <  4:158543/148‑1 (MQ=255)
                                                                   aGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCa                                                                                >  3:153700/1‑149 (MQ=255)
                                                                         gTTCGCCGAGGTCTGATCGGACCCCCCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCg                                                                            >  1:29393/1‑147 (MQ=255)
                                                                                    tCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCggctg                                                               <  2:230963/149‑1 (MQ=255)
                                                                                        aTCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCggctggc                                                             <  3:72605/147‑1 (MQ=255)
                                                                                         tCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAg                                                                         <  4:31914/134‑1 (MQ=255)
                                                                                         tCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAg                                                                         <  4:27042/134‑1 (MQ=255)
                                                                                         tCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAg                                                                         >  3:27042/1‑134 (MQ=255)
                                                                                         tCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAg                                                                         >  3:31914/1‑134 (MQ=255)
                                                                                            gACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCt                                                                                                 >  7:163856/1‑107 (MQ=255)
                                                                                            gACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCt                                                                                                 <  8:163856/107‑1 (MQ=255)
                                                                                               ccACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCt                                                                                            >  5:382328/1‑109 (MQ=255)
                                                                                               ccACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCt                                                                                            <  6:382328/109‑1 (MQ=255)
                                                                                                 aCCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGtt                                                  >  1:142370/1‑149 (MQ=255)
                                                                                                   cGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGtt                                                  >  3:242515/1‑147 (MQ=255)
                                                                                                   cGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGa                                                <  7:79037/149‑1 (MQ=255)
                                                                                                        ccACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAgcgc                                            >  2:205863/1‑148 (MQ=255)
                                                                                                            tcgCAGGTGGTCGACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCt                                       >  5:15700/3‑149 (MQ=255)
                                                                                                                     gACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCtt                              >  7:53364/1‑149 (MQ=255)
                                                                                                                             gCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGt                            >  8:353759/1‑143 (MQ=255)
                                                                                                                             gCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGt                            <  7:353759/143‑1 (MQ=255)
                                                                                                                              cGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCg                     >  5:194147/1‑149 (MQ=255)
                                                                                                                                      tgCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCa                                          <  3:250324/120‑1 (MQ=255)
                                                                                                                                      tgCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTAGATACCACTCGAAGGCCTCGGCTGGCTCGCAAAGGTTGAGCGCCa                                          >  4:250324/1‑120 (MQ=255)
                                                                                                                                        cAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCgg            <  8:262169/148‑1 (MQ=255)
                                                                                                                                          gCCGCTTCAGTGTAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCTATGGCTGTGCTCACTCGATACCTCTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGtt         <  8:163208/149‑1 (MQ=255)
                                                                                                                                          gCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCAGTCGCAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGtt         <  2:29393/149‑1 (MQ=255)
                                                                                                                                          gCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGtt         <  6:194147/149‑1 (MQ=255)
                                                                                                                                             gCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCgg            >  3:58160/1‑143 (MQ=255)
                                                                                                                                             gCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCgg            <  4:58160/143‑1 (MQ=255)
                                                                                                                                             gCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCa      >  8:320446/1‑149 (MQ=255)
                                                                                                                                              cTTCAGTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCAt     >  2:108110/1‑149 (MQ=255)
                                                                                                                                              cTTCAGTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCAt     <  1:108110/149‑1 (MQ=255)
                                                                                                                                                 cACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATc    <  2:142370/147‑1 (MQ=255)
                                                                                                                                                  aCTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATCCa  <  6:15700/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CATTCTCGAAGGGAAGGTTATGGCTGAACAGTACGTGGAAAGTCAGAAAAGCGGAGCTTGTGCTGCCAGGTGAGTTCGCCGAGGTCTGATCGGACCCACCGCAACCACGGGCAGGTCGACCACCTGCGAATCTGTGCAGCCGCTTCACTGAAGCGGCTCTATCAGCTCAGCGCCTTGATTGCGAACATTGGCAATGGCTGTGCTCACTCGATACCACTCGAAGGCCTCGGCTGGCTCGCCCAGGTTGAGCGCCATCTGCTCGGCTTGTTCGCCCGGCGTCGCGGGTTGCATCCA  >  NC_002947/3344199‑3344492

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: ‑