Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,792,826 C→T *286* (TAG→TAA yibB ← YibB family protein, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,792,8260CT98.8% 273.0 / NA 86*286* (TAG→TAAyibBYibB family protein, function unknown
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (45/40);  total (45/41)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATA  >  NC_000913/3792692‑3792960
                                                                                                                                      |                                                                                                                                      
aaaaCTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTAtt                                                                                                                                    >  2:123288/1‑139 (MQ=255)
    cTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCt                                                                                                                                >  2:795252/1‑139 (MQ=255)
       tACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGa                                                                                                                             <  2:307305/139‑1 (MQ=255)
                    gggCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATtctc                                                                                                                >  1:785497/1‑139 (MQ=255)
                      gCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAg                                                                                                              >  1:1244215/1‑139 (MQ=255)
                                 agtgattTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGAGCCTGCAttg                                                                                                   >  2:1250571/3‑138 (MQ=255)
                                  ttgattTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCAtttt                                                                                                  <  1:123288/139‑1 (MQ=255)
                                    gattTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTg                                                                                                <  1:885637/139‑1 (MQ=255)
                                    gattTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTg                                                                                                <  2:1244215/139‑1 (MQ=255)
                                     attTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGc                                                                                               <  2:944862/139‑1 (MQ=255)
                                     attTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGc                                                                                               <  1:1452275/139‑1 (MQ=255)
                                      ttttACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGcc                                                                                              >  2:799230/1‑139 (MQ=255)
                                         tACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGa                                                                                                                             <  1:121558/105‑1 (MQ=255)
                                         tACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGa                                                                                                                             >  2:121558/1‑105 (MQ=255)
                                           cTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTaa                                                                                         >  1:1228759/1‑139 (MQ=255)
                                              cATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAgt                                                                                      >  2:781605/1‑139 (MQ=255)
                                                 tGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGtt                                                                                   <  2:423508/139‑1 (MQ=255)
                                                    tataAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACt                                                                                >  1:428653/1‑139 (MQ=255)
                                                        aaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAgtg                                                                                     >  2:379312/1‑130 (MQ=255)
                                                        aaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAgtg                                                                                     <  1:379312/130‑1 (MQ=255)
                                                           aGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGTAAGGAAAATTCTCAGCTCCTGCATTTTTGCCCCTAAAGTGTTACTCATGAag                                                                         >  1:1406434/1‑139 (MQ=255)
                                                             aaaaaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCAtttt                                                                                                  <  2:386547/112‑1 (MQ=255)
                                                             aaaaaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCAtttt                                                                                                  >  1:386547/1‑112 (MQ=255)
                                                             aaaaaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAagca                                                                       >  1:954491/1‑139 (MQ=255)
                                                             aaaaaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAagca                                                                       >  2:1175123/1‑139 (MQ=255)
                                                                            aGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGaa                                                        >  1:193375/1‑139 (MQ=255)
                                                                            aGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGaa                                                        <  2:428653/139‑1 (MQ=255)
                                                                                   tttttCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTg                                                                                                <  2:134666/92‑1 (MQ=255)
                                                                                   tttttCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTg                                                                                                >  1:134666/1‑92 (MQ=255)
                                                                                    ttttCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTc                                                <  1:799230/139‑1 (MQ=255)
                                                                                        ccTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGAGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCctct                                            >  2:853601/1‑139 (MQ=255)
                                                                                         cTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCtctc                                           <  2:1228759/139‑1 (MQ=255)
                                                                                         cTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCtctc                                           <  2:193375/139‑1 (MQ=255)
                                                                                          tATCTCACTTCGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTCTCTCGAAAGGAAAATTCTCAGCGCCTCCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTcc                                          <  2:1406434/139‑1 (MQ=255)
                                                                                           aTCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTcc                                               >  2:822175/1‑133 (MQ=255)
                                                                                           aTCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTcc                                               <  1:822175/133‑1 (MQ=255)
                                                                                             ctcAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCa                                       >  1:220015/1‑139 (MQ=255)
                                                                                             ctcAGTTGGTATCAAAACGCTACCCAATAATGGAGTGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCa                                       <  1:255096/139‑1 (MQ=255)
                                                                                                  ttGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGa                                                                                                                             <  1:349783/48‑1 (MQ=255)
                                                                                                  ttGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGa                                                                                                                             >  2:349783/1‑48 (MQ=255)
                                                                                                   tGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGt                                 <  1:1175123/139‑1 (MQ=255)
                                                                                                   tGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGt                                 <  2:220015/139‑1 (MQ=255)
                                                                                                      tATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTaa                                                                                         >  1:1022971/1‑80 (MQ=255)
                                                                                                      tATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTaa                                                                                         <  2:1022971/80‑1 (MQ=255)
                                                                                                      tATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTa                              >  1:133167/1‑139 (MQ=255)
                                                                                                      tATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTa                              >  1:115182/1‑139 (MQ=255)
                                                                                                        tCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAgg                            <  2:954491/139‑1 (MQ=255)
                                                                                                        tCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAgg                            >  1:1178813/1‑139 (MQ=255)
                                                                                                          aaaaCGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGtt                          >  1:351950/1‑139 (MQ=255)
                                                                                                          aaaaCGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGtt                          >  1:724675/1‑139 (MQ=255)
                                                                                                           aaaCGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGaaa                                                                 <  2:683864/99‑1 (MQ=255)
                                                                                                           aaaCGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGaaa                                                                 >  1:683864/1‑99 (MQ=255)
                                                                                                              cGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTCTCCTCTCCCCAGATAGTTTAGGTTGAag                      <  2:115182/139‑1 (MQ=255)
                                                                                                                cTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAg                                      >  1:760865/1‑121 (MQ=255)
                                                                                                                cTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAg                                      <  2:760865/121‑1 (MQ=255)
                                                                                                                cTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGt                                 >  1:969623/1‑126 (MQ=255)
                                                                                                                cTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGt                                 <  2:969623/126‑1 (MQ=255)
                                                                                                                cTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAagag                    >  1:1085851/1‑139 (MQ=255)
                                                                                                                      aCTAATTGAGGGAAACTTATTTTCTTGAAAGGAAAAGTCTCAGCGCCTGAATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCTTTTTCCTCGCCCCAGATAGTTTAGGTTGAAGAGATCggg              <  1:1240948/139‑1 (MQ=255)
                                                                                                                      aCTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCggg              <  2:228641/139‑1 (MQ=255)
                                                                                                                      aCTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCggg              <  2:785497/139‑1 (MQ=255)
                                                                                                                        tAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGTTCGGGGc            <  2:1085851/139‑1 (MQ=255)
                                                                                                                        tAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGc            <  1:795252/139‑1 (MQ=255)
                                                                                                                           tGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGtt         <  1:469517/139‑1 (MQ=255)
                                                                                                                             gAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGcccc                                                                                            <  2:1078073/54‑1 (MQ=255)
                                                                                                                             gAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGcccc                                                                                            >  1:1078073/1‑54 (MQ=255)
                                                                                                                             gAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTg       >  1:1176405/1‑139 (MQ=255)
                                                                                                                               gggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCg                                                          >  1:649048/1‑86 (MQ=255)
                                                                                                                               gggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCg                                                          <  2:649048/86‑1 (MQ=255)
                                                                                                                               gggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTa     >  2:1426691/1‑139 (MQ=255)
                                                                                                                               gggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTa     >  1:1262293/1‑139 (MQ=255)
                                                                                                                                ggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTaa                                                                                         >  2:202917/1‑54 (MQ=255)
                                                                                                                                ggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTaa                                                                                         >  1:541391/1‑54 (MQ=255)
                                                                                                                                ggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTaa                                                                                         <  2:541391/54‑1 (MQ=255)
                                                                                                                                ggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTaa                                                                                         <  1:202917/54‑1 (MQ=255)
                                                                                                                                ggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCctc                                             >  1:1205257/1‑98 (MQ=255)
                                                                                                                                ggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCctc                                             <  2:1205257/98‑1 (MQ=255)
                                                                                                                                ggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGtt         <  1:1393189/134‑1 (MQ=255)
                                                                                                                                ggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGtt         >  2:1393189/1‑134 (MQ=255)
                                                                                                                                ggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTaa    >  2:982364/1‑139 (MQ=255)
                                                                                                                                 gAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCctc                                             <  2:179213/97‑1 (MQ=255)
                                                                                                                                 gAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCctc                                             >  1:179213/1‑97 (MQ=255)
                                                                                                                                  aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGt          >  2:985495/1‑131 (MQ=255)
                                                                                                                                  aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGt          <  1:985495/131‑1 (MQ=255)
                                                                                                                                  aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATa  >  1:602849/1‑139 (MQ=255)
                                                                                                                                  aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGAAAGTTTAGGTTGAAGAGATCGGGCCGTTTGTAATa  >  2:1273209/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                      
AAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATA  >  NC_000913/3792692‑3792960

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: