Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,938,222 T→A 5.4% F75Y (TTT→TAT) ‡ PP_3471 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,938,2220TA 5.4% 224.7 / 2.5 74F75Y (TTT→TAT) ‡PP_3471hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/40);  new base A (2/2);  total (32/42)
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.30e-01

ATAAAACATCCGATGGATATTTACTCGGGAATGGT‑CGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAA  >  NC_002947/3938075‑3938369
                                                                                                                                                    |                                                                                                                                                   
aTAAATCATCCGATGGATATTTACTCGGGAATGGT‑CGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACttt                                                                                                                                                    <  5:240552/149‑1 (MQ=255)
aTAAAACATCCGATGGATATTTACTCGGGAATGGTACGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACtt                                                                                                                                                     >  1:55196/1‑149 (MQ=255)
aTAAAACATCCGATGGATATTTACTCGGGAATGGT‑CGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACttt                                                                                                                                                    >  8:86290/1‑149 (MQ=255)
aTAAAACATCCGATGGATATTTACTCGGGAATGGT‑CGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACttt                                                                                                                                                    <  5:364420/149‑1 (MQ=255)
 tAAAACATCCGATGGATATTTACTCGGGAATGGT‑CGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTttg                                                                                                                                                   <  4:61436/149‑1 (MQ=255)
   aaaCATCCGATGGATATTTACTCGGGAATGGT‑CGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTttgt                                                                                                                                                  <  3:208128/148‑1 (MQ=255)
       aTCCGATGGATATTTACTCGGGAATGGT‑CGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATc                                                                                                                                             >  2:185413/1‑149 (MQ=255)
            aTGGATATTTAGTCGGGTATGGT‑CGGCGCGGATACAAGCCGCGCTTGTTGCGTTTTTGTTCGCCTGACAATTTAAGTCCATTTCAATCAGGCGTCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCATCTTTGTTGATCCAAGt                                                                                                                                        <  8:299219/149‑1 (MQ=255)
            aTGGATATTTACTCGGGAATGGT‑CGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGt                                                                                                                                        <  6:223677/149‑1 (MQ=255)
            aTGGATATTTACTCGGGAATGGT‑CGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACAAAGTTGATcggatc                                                                                                                                        <  6:13338/149‑6 (MQ=255)
               gATATTTACTCGGGAATGGT‑CGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGc                                                                                                                                     <  8:252142/149‑1 (MQ=255)
                   tttACTCGGGAATGGT‑CGGCGCTGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCAtt                                                                                                                                 >  7:55951/1‑149 (MQ=255)
                     tACTCGGGAATGGT‑CGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTc                                                                                                                                <  1:185413/148‑1 (MQ=255)
                       cTCGGGAATGGT‑CGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCg                                                                                                                              <  7:86290/148‑1 (MQ=255)
                       cTCGGGAATGGT‑CGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCg                                                                                                                              <  2:215938/148‑1 (MQ=255)
                        tCGGGAATGGT‑CGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGcc                                                                                                                                    >  6:2650/1‑141 (MQ=255)
                        tCGGGAATGGT‑CGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGcc                                                                                                                                    <  5:2650/141‑1 (MQ=255)
                        tCGGGAATGGT‑CGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGc                                                                                                                            >  4:11988/1‑149 (MQ=255)
                        tCGGGAATGGT‑CGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGc                                                                                                                            >  8:88741/1‑149 (MQ=255)
                          gggAATGGT‑CGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCggn                                                                                                                            <  2:309782/147‑2 (MQ=255)
                               tGGT‑CGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGc                                                                                                                     <  4:104053/149‑1 (MQ=255)
                                  t‑CGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCAc                                                                                                                  >  7:82046/1‑149 (MQ=255)
                                     ggCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATcc                                                                                                                                            >  8:156111/1‑121 (MQ=255)
                                       cgcgGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTAc                                                                                                               >  8:283241/1‑148 (MQ=255)
                                        gcgGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGAGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACaa                                                                                                             <  5:369383/149‑1 (MQ=255)
                                           gATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAAt                                                                                                           >  8:251969/1‑148 (MQ=255)
                                           gATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACAAAGTTGATcgg                                                                                                                                           <  8:249719/116‑3 (MQ=255)
                                           gATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACAAAGTTGATcgg                                                                                                                                           >  7:249719/1‑114 (MQ=255)
                                             tACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCt                                                                                                        <  4:295707/149‑1 (MQ=255)
                                               cAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTgg                                                                                                      <  7:88741/149‑1 (MQ=255)
                                               cAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTACATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTgg                                                                                                      <  8:100852/149‑1 (MQ=255)
                                                     gcgcTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTgg                                                                                                                                      <  2:97977/111‑1 (MQ=255)
                                                     gcgcTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTgg                                                                                                                                      >  1:97977/1‑111 (MQ=255)
                                                         ttGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCaaa                                                                                             >  4:156507/1‑148 (MQ=255)
                                                            aTGCGTTTTAGTTCGCCGGACAATTTAAGTCCATTTAAATCAGGCGGCATCAACAGCTACGCCTACTGCAAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGCCGATCTGCCAATACAAATTCTGACCGCGCAAAAAt                                                                                          >  1:265773/1‑148 (MQ=255)
                                                             tGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTAGCCGCGCAAAAATTa                                                                                        >  1:361847/1‑149 (MQ=255)
                                                             tGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTAGCCGCGCAAAAATTa                                                                                        <  2:361847/149‑1 (MQ=255)
                                                              gCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCcgc                                                                                                  >  6:167258/1‑138 (MQ=255)
                                                              gCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCcgc                                                                                                  <  5:167258/138‑1 (MQ=255)
                                                               cggttttGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCGAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTAcc                                                                                      <  7:291318/146‑1 (MQ=255)
                                                               cGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTAcc                                                                                      <  4:92916/149‑1 (MQ=255)
                                                               cGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTAcc                                                                                      >  1:208467/1‑149 (MQ=255)
                                                                      gTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATGAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTc                                                                                <  5:320221/148‑1 (MQ=255)
                                                                           ccGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAg                                                                          <  4:13880/149‑1 (MQ=255)
                                                                             ggAAAATTTAAGTCCATTTCAATCAGGGGGCATCAACATCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGc                                                                        <  2:265773/149‑1 (MQ=255)
                                                                               aCAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGc                                                                        >  1:166894/1‑147 (MQ=255)
                                                                                  aTTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGa                                                                   <  8:82046/149‑1 (MQ=255)
                                                                                   tttAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGa                                                                   >  5:374883/1‑148 (MQ=255)
                                                                                       aGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAAtt                                                               <  7:251969/148‑1 (MQ=255)
                                                                                        gTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAAttt                                                              >  4:274949/1‑148 (MQ=255)
                                                                                         tCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTaa                                                            <  6:105282/149‑1 (MQ=255)
                                                                                           cATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCACAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAag                                                           <  7:270449/148‑1 (MQ=255)
                                                                                             tttCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAagag                                                         >  8:59704/1‑148 (MQ=255)
                                                                                                 aaTCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATcc                                                    <  2:55196/149‑1 (MQ=255)
                                                                                                   tCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCaa                                                  <  7:283241/149‑1 (MQ=255)
                                                                                                     aGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAAc                                                 <  8:55951/148‑1 (MQ=255)
                                                                                                          gcATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCgg                                           >  2:137371/1‑149 (MQ=255)
                                                                                                               aaCAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATc                                      >  5:157082/1‑149 (MQ=255)
                                                                                                                aCAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATcc                                     <  7:16386/149‑1 (MQ=255)
                                                                                                                aCAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATcc                                     <  7:22007/149‑1 (MQ=255)
                                                                                                                aCAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATcc                                     >  2:129974/1‑149 (MQ=255)
                                                                                                                       cGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTaa                              >  5:184638/1‑149 (MQ=255)
                                                                                                                          cTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGt                           >  4:61323/1‑149 (MQ=255)
                                                                                                                             cTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGttt                         <  1:137371/148‑1 (MQ=255)
                                                                                                                             cTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGtt                          <  1:129974/147‑1 (MQ=255)
                                                                                                                                ccAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCTGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTtgat                      <  7:312229/148‑1 (MQ=255)
                                                                                                                                ccAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTtga                       <  3:274949/147‑1 (MQ=255)
                                                                                                                                ccAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCAAATTAATCTTtgat                      >  1:67297/1‑148 (MQ=255)
                                                                                                                                    ggAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGAtt                  >  3:292942/1‑148 (MQ=255)
                                                                                                                                       gaTCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACa               >  2:171610/1‑148 (MQ=255)
                                                                                                                                        cttggATCAACAAAGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACAt              >  6:86999/5‑148 (MQ=255)
                                                                                                                                           cGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTtgatga                    >  7:4587/1‑139 (MQ=255)
                                                                                                                                           cGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTtgatga                    <  8:4587/139‑1 (MQ=255)
                                                                                                                                           cGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATaaaa          <  6:374883/149‑1 (MQ=255)
                                                                                                                                           cGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATaaaa          <  2:208467/149‑1 (MQ=255)
                                                                                                                                            gATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTTGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAaga                                                          <  3:266619/100‑1 (MQ=255)
                                                                                                                                            gATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTTGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAaga                                                          >  4:266619/1‑100 (MQ=255)
                                                                                                                                              tCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTa       <  7:59704/149‑1 (MQ=255)
                                                                                                                                                 aCTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTAtt     >  4:8182/1‑148 (MQ=255)
                                                                                                                                                  cTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCTCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTaa   <  5:327583/149‑1 (MQ=255)
                                                                                                                                                   tttGTTGATCCAAGTGGCGATTACGGCGTTCTGCCACTACAAATTCTGGCAGCGCAAAAATTACCGACTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAAACTGCGGTCATCCCAATTAATGATGGATGATTACATAAAACtaataa   >  4:314697/1‑148 (MQ=255)
                                                                                                                                                   tttGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTaaa  <  3:11988/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATAAAACATCCGATGGATATTTACTCGGGAATGGT‑CGGCGCGGATACAAGCCGCGCTTGATGCGTTTTTGTTCGCCGGACAATTTAAGTCCATTTCAATCAGGCGGCATCAACAGCTACGCCTACTGCCAAGGAGATCCGATCAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAA  >  NC_002947/3938075‑3938369

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