Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,699,528 G→A 5.1% V439M (GTG→ATG)  bcsG → DUF3260 family cellulose production inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,699,5280GA 5.1% 266.1 / 2.4 78V439M (GTG→ATG) bcsGDUF3260 family cellulose production inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (34/40);  new base A (2/2);  total (36/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 = 4.80e-01

ACAACACGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTG  >  NC_000913/3699382‑3699672
                                                                                                                                                  |                                                                                                                                                
acaacaCGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCgt                                                                                                                                                 >  6:216554/1‑148 (MQ=255)
 caacaCGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCgtg                                                                                                                                                <  5:322527/148‑1 (MQ=255)
 caacaCGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCg                                                                                                                                                  <  7:124250/146‑1 (MQ=255)
 caacaCGCTTCCACTGCATGACGGCAACCATTATCCGGGGGGCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCg                                                                                                                                                  <  2:26215/146‑1 (MQ=255)
  aacaCGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCgtgg                                                                                                                                               >  7:43008/1‑148 (MQ=255)
        cTTCCACTGCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCgg                                                                                                                                         <  6:31097/148‑1 (MQ=255)
              cTGCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGGCGTGGTGCCGGAACAcg                                                                                                                                   >  8:257972/1‑148 (MQ=255)
               tGCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACAcgg                                                                                                                                  >  5:140656/1‑148 (MQ=255)
                gCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACAcggc                                                                                                                                 >  5:286568/1‑148 (MQ=255)
                    taCGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGgcggcg                                                                                                                             >  7:242661/2‑148 (MQ=255)
                    gACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCgtg                                                                                                                                                >  2:43864/1‑129 (MQ=255)
                    gACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCgtg                                                                                                                                                <  1:43864/129‑1 (MQ=255)
                    gACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGgcggcg                                                                                                                             >  7:63069/1‑148 (MQ=255)
                         cAACCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTg                                                                                                                        <  8:63069/148‑1 (MQ=255)
                         cAACCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGGGGTGCCGGAACACGGCGGCGCGCTg                                                                                                                        <  5:369651/148‑1 (MQ=255)
                           aCCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGaa                                                                                                                      <  3:1341/148‑1 (MQ=255)
                            ccATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAg                                                                                                                     >  7:376465/1‑148 (MQ=255)
                                 aTCCGGGGGTCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTg                                                                                                                        <  4:157788/140‑1 (MQ=255)
                                 aTCCGGGGGTCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTg                                                                                                                        >  3:157788/1‑140 (MQ=255)
                                   ccGGGGGTCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGa                                                                                                                >  1:235470/1‑146 (MQ=255)
                                     gggggTCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGa                                                                                                            <  1:237074/148‑1 (MQ=255)
                                          tCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATg                                                                                                         >  2:141350/1‑146 (MQ=255)
                                          tCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTCAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCa                                                                                                       <  1:175933/148‑1 (MQ=255)
                                          tCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTCAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCa                                                                                                       >  2:175933/1‑148 (MQ=255)
                                           cAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTg                                                                                                                        <  2:312744/130‑1 (MQ=255)
                                           cAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTg                                                                                                                        >  1:312744/1‑130 (MQ=255)
                                           cAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGa                                                                                                                       <  1:286707/131‑1 (MQ=255)
                                           cAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGa                                                                                                                       >  2:286707/1‑131 (MQ=255)
                                             gCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAgg                                                                                                     >  7:61920/1‑147 (MQ=255)
                                              cAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGa                                                                                                                                        <  3:268222/111‑1 (MQ=255)
                                              cAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGa                                                                                                                                        >  4:268222/1‑111 (MQ=255)
                                              cAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCa                                                                                                       >  5:170018/1‑144 (MQ=255)
                                              cAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCa                                                                                                       <  6:170018/144‑1 (MQ=255)
                                              cAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTa                                                                                                   >  4:165476/1‑148 (MQ=255)
                                              cAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTa                                                                                                   <  2:345167/148‑1 (MQ=255)
                                               aaaaCAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTa                                                                                                   <  6:140656/147‑1 (MQ=255)
                                                         aTTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACg                                                                                        <  1:350546/148‑1 (MQ=255)
                                                          ttACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGt                                                                                       >  8:145502/1‑148 (MQ=255)
                                                           tACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTg                                                                                      <  3:165476/148‑1 (MQ=255)
                                                            aCAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTg                                                                                      <  2:235470/147‑1 (MQ=255)
                                                             cAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGat                                                                                    <  6:286568/148‑1 (MQ=255)
                                                             cAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGat                                                                                    <  8:43008/148‑1 (MQ=255)
                                                             cAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGGCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGat                                                                                    <  5:216554/148‑1 (MQ=255)
                                                              aaaGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGata                                                                                   <  2:10716/148‑1 (MQ=255)
                                                                aGCGCTGGCGCAGAAATTCTTTGATTAACTGGACGCATTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCATGGTGCCGGAACACGGCGGCGCGCTGAAGTGCGACAGAATGCAGGTATCTGGCCTACGTg                                                                                      <  6:360520/143‑1 (MQ=255)
                                                                aGCGCGGGCGAAGAATTTCTTTGATGAACTTGACGACTTCTTTACTGAACTTGAGAAATCGGTTCGTAAAGTGATGGTGGTCATGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTg                                                                                      <  6:366314/143‑1 (MQ=255)
                                                                 gcgcGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTg                                                                                                                        >  1:364189/1‑108 (MQ=255)
                                                                 gcgcGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTg                                                                                                                        <  2:364189/108‑1 (MQ=255)
                                                                 gcgcGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATcc                                                                                <  8:376465/148‑1 (MQ=255)
                                                                 gcgcGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCATGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGTCCTACGTGa                                                                                     >  5:360520/1‑143 (MQ=255)
                                                                 gcgcGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCATGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGa                                                                                     >  5:366314/1‑143 (MQ=255)
                                                                  cgcgGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGGCGTAACGTGATGGTGGGCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTCGCTGGCCTACGTGATATccc                                                                               <  2:272046/148‑1 (MQ=255)
                                                                      ggCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGc                                                                           <  1:141350/148‑1 (MQ=255)
                                                                            gAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTc                                                                      >  3:22413/1‑147 (MQ=255)
                                                                                ttCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCa                                                                 >  8:324043/1‑148 (MQ=255)
                                                                                     tgatgaACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCg                                                              <  8:242661/146‑1 (MQ=255)
                                                                                     tgatgaACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGa                                                             <  8:61920/147‑1 (MQ=255)
                                                                                         gaACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACATcc                                                        >  4:341863/1‑148 (MQ=255)
                                                                                           aCTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTcccc                                                      >  3:143944/1‑148 (MQ=255)
                                                                                                aCGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTgg                                                 >  7:63586/1‑148 (MQ=255)
                                                                                                aCGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTgg                                                 <  4:22413/148‑1 (MQ=255)
                                                                                                   ccTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGt                                              >  6:367273/1‑148 (MQ=255)
                                                                                                   ccTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGGCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGt                                              <  7:145502/148‑1 (MQ=255)
                                                                                                    cttcttTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTg                                             <  7:324043/148‑1 (MQ=255)
                                                                                                       cttTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGaaa                                          >  6:103/1‑148 (MQ=255)
                                                                                                       cttTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGaa                                           <  2:231675/147‑1 (MQ=255)
                                                                                                        tttACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAAt                                         >  6:364477/1‑148 (MQ=255)
                                                                                                         ttACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGACCACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAAtt                                        >  6:277160/1‑148 (MQ=255)
                                                                                                              gAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATc                                                                                                 <  6:126999/86‑1 (MQ=255)
                                                                                                              gAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATc                                                                                                 >  5:126999/1‑86 (MQ=255)
                                                                                                              gAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCg                                   >  7:312393/1‑148 (MQ=255)
                                                                                                                  ttGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCAt                               >  6:18885/1‑148 (MQ=255)
                                                                                                                    gagaAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATg                              <  4:143944/147‑1 (MQ=255)
                                                                                                                    gagaAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATg                              <  5:103/147‑1 (MQ=255)
                                                                                                                    gagaAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATg                              <  5:18885/147‑1 (MQ=255)
                                                                                                                      gaAATCGGGTCGTAAAGTGATGGTGGGCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACATCCCCGTTGGGGTGAAATTCTTCGGCATGAAg                           <  3:341863/148‑1 (MQ=255)
                                                                                                                        aaTCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAAt                                         <  1:364478/132‑1 (MQ=255)
                                                                                                                        aaTCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAAt                                         >  2:364478/1‑132 (MQ=255)
                                                                                                                          tCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGCAc                       <  8:63586/148‑1 (MQ=255)
                                                                                                                                cGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGCACCGCATc                 <  5:277160/148‑1 (MQ=255)
                                                                                                                                   aaaGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCgg                                  >  3:189244/1‑128 (MQ=255)
                                                                                                                                   aaaGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCgg                                  <  4:189244/128‑1 (MQ=255)
                                                                                                                                      gTGATGGTGGTCGTGGGGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGCACCGCATCAggggg           <  8:312393/148‑1 (MQ=255)
                                                                                                                                               gtCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTg  >  1:40772/1‑148 (MQ=255)
                                                                                                                                               gtCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTg  >  1:79342/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
ACAACACGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCTTTGATGAACTGGACGCCTTCTTTACTGAACTTGAGAAATCGGGTCGTAAAGTGATGGTGGTCGTGGTGCCGGAACACGGCGGCGCGCTGAAGGGCGACAGAATGCAGGTATCTGGCCTACGTGATATCCCTAGCCCGTCTATCACCGACGTCCCCGTTGGGGTGAAATTCTTCGGCATGAAGGCACCGCATCAGGGGGCACCGATTG  >  NC_000913/3699382‑3699672

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: