Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,024,887 C→T S179S (AGC→AGT xanQ → xanthine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,024,8870CT98.8% 252.1 / ‑4.9 81S179S (AGC→AGTxanQxanthine permease
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (44/36);  total (44/37)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CTTCATTTATCTTGCCCTATTTACGTCGGGTTATTACGCCTACCGTCAGCGGTATTGTGGTACTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGCGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGATCGGCTTTAACTGCTGTCGCAGTCCGTTGCTACGCATGGGAGGGATCGCCATTGGGCTATGTGTCGGCTATATCGCA  >  NC_000913/3024754‑3025022
                                                                                                                                     |                                                                                                                                       
cTTCATTTATCTTGCCCTATTTACGTCGGGTTATTACGCCTACCGTCAGCGGTATTGTGGTACTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTAc                                                                                                                                    >  1:773034/1‑139 (MQ=255)
             gCCCTATTTACGTCGGGTTATTACGCCTACCGTCAGCGGTATTGTGGTACTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTAc                                                                                                                       >  2:273871/1‑139 (MQ=255)
             gCCCTATTTACGTCGGGTTATTACGCCTACCGTCAGCGGTATTGTGGTACTGATGATCGGCTTAAGCCTGATTAAAGTCGACATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTAc                                                                                                                       >  1:1114601/1‑139 (MQ=255)
                   tttACGTCGGGTTATTACGCCTACCGTCAGCGGTATTGTGGTACTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACAt                                                                                                                 <  2:773034/139‑1 (MQ=255)
                        gTCGGGTTATTACGCCTACCGTCAGCGGTATTGTGGTACTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCaa                                                                                                                           <  1:62870/124‑1 (MQ=255)
                        gTCGGGTTATTACGCCTACCGTCAGCGGTATTGTGGTACTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCaa                                                                                                                           >  2:62870/1‑124 (MQ=255)
                        gTCGGGTTATTACGCCTACCGTCAGCGGTATTGTGGTACTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACg                                                                                                                      >  1:1127071/1‑129 (MQ=255)
                        gTCGGGTTATTACGCCTACCGTCAGCGGTATTGTGGTACTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACg                                                                                                                      <  2:1127071/129‑1 (MQ=255)
                        gTCGGGTTATTACGCCTACCGTCAGCGGTATTGTGGTACTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCgg                                                                                                            >  2:998275/1‑139 (MQ=255)
                            ggTTATTACGCCTACCGTCAGCGGTATTGTGGTACTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGtt                                                                                                        >  2:749578/1‑139 (MQ=255)
                            ggTTATTACGCCTACCGTCAGCGGTATTGTGGTACTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGtt                                                                                                        >  2:494310/1‑139 (MQ=255)
                             gTTATTACGCCTACCGTCAGCGGTATTGTGGTACTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAAtt                                                                                                                         >  1:1240675/1‑121 (MQ=255)
                             gTTATTACGCCTACCGTCAGCGGTATTGTGGTACTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAAtt                                                                                                                         <  2:1240675/121‑1 (MQ=255)
                             gTTATTACGCCTACCGTCAGCGGTATTGTGGTACTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTg                                                                                                       >  1:1237873/1‑139 (MQ=255)
                             gTTATTACGCCTACCGTCAGCGGTATTGTGGTACTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTg                                                                                                       >  2:1072854/1‑139 (MQ=255)
                                    cGCCTACCGTCAGCGGTATTGTGGTACTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTAtt                                                                                                <  1:1025317/139‑1 (MQ=255)
                                       cTACCGTCAGCGGTATTGTGGTACTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGaa                                                                                                                    <  2:1074869/116‑1 (MQ=255)
                                       cTACCGTCAGCGGTATTGTGGTACTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGaa                                                                                                                    >  1:1074869/1‑116 (MQ=255)
                                       cTACCGTCAGCGGTATTGTGGTACTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCgg                                                                                                            <  2:945012/124‑1 (MQ=255)
                                       cTACCGTCAGCGGTATTGTGGTACTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCgg                                                                                                            >  1:945012/1‑124 (MQ=255)
                                       cTACCGTCAGCGGTATTGTGGTACTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGTGTTGGTTTATTGGt                                                                                             <  1:998275/139‑1 (MQ=255)
                                       cTACCGTCAGCGGTATTGTGGTACTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGt                                                                                             <  1:1072854/139‑1 (MQ=255)
                                         aCCGTCAGCGGTATTGTGGTACTCATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGttt                                                                                           >  2:1341314/1‑139 (MQ=255)
                                                     aTTGTGGTACTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGa                                                                               >  1:562636/1‑139 (MQ=255)
                                                     aTTGTGGTACTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGa                                                                               <  2:352177/139‑1 (MQ=255)
                                                      ttGTGGTACTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGAt                                                                              <  1:749578/139‑1 (MQ=255)
                                                           gTACTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTAtt                                                                                                >  1:594005/1‑116 (MQ=255)
                                                           gTACTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTAtt                                                                                                <  2:594005/116‑1 (MQ=255)
                                                           gTACTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAAtt                                                                                      <  1:430751/126‑1 (MQ=255)
                                                           gTACTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAAtt                                                                                      >  2:430751/1‑126 (MQ=255)
                                                            tACTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGATCGGCtt                                                                        >  1:834082/1‑139 (MQ=255)
                                                              cTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTAc                                                                                                                                    <  2:383340/77‑1 (MQ=255)
                                                              cTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTAc                                                                                                                                    >  1:383340/1‑77 (MQ=255)
                                                              cTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGATCGGCTTTa                                                                      >  1:1315382/1‑139 (MQ=255)
                                                              cTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGATCGGCTTTa                                                                      >  2:1186080/1‑139 (MQ=255)
                                                                gatgatCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTa                                                                                                  <  2:1264403/109‑1 (MQ=255)
                                                                gatgatCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTa                                                                                                  >  1:1264403/1‑109 (MQ=255)
                                                                gatgatCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGAt                                                                              >  2:148842/1‑129 (MQ=255)
                                                                gatgatCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGAt                                                                              <  2:172435/129‑1 (MQ=255)
                                                                gatgatCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGAt                                                                              >  1:172435/1‑129 (MQ=255)
                                                                gatgatCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGAt                                                                              <  1:148842/129‑1 (MQ=255)
                                                                gatgatCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGATCGGCTTTAAc                                                                    >  2:510646/1‑139 (MQ=255)
                                                                gatgatCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGATCGGCTTTAAc                                                                    >  2:674618/1‑139 (MQ=255)
                                                                 atgatCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGATCGGCTTTAAct                                                                   >  2:753275/1‑139 (MQ=255)
                                                                   gatCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGATCGGCTTTAActgc                                                                 >  1:790237/1‑139 (MQ=255)
                                                                     tCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGa                                                                                                                     <  2:409909/85‑1 (MQ=255)
                                                                     tCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGa                                                                                                                     >  1:409909/1‑85 (MQ=255)
                                                                          ttAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGATCGGCTTTAACTGCTGTCGCa                                                          <  1:1186080/139‑1 (MQ=255)
                                                                                  gATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGATCGGCTTTAACTGCTGTCGCAGTCCGTTg                                                  >  1:1175130/1‑139 (MQ=255)
                                                                                   aTTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTAt                                                                                                 >  1:615744/1‑91 (MQ=255)
                                                                                   aTTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTAt                                                                                                 <  2:615744/91‑1 (MQ=255)
                                                                                   aTTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGAt                                                                              >  1:273294/1‑110 (MQ=255)
                                                                                   aTTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGAt                                                                              <  2:273294/110‑1 (MQ=255)
                                                                                   aTTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGATCGGCTTTAActg                                                                  >  1:685178/1‑122 (MQ=255)
                                                                                   aTTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGATCGGCTTTAActg                                                                  <  2:685178/122‑1 (MQ=255)
                                                                                   aTTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGATCGGCTTTAACTGCTGTCGCAGTCCGTTGc                                                 >  1:1200729/1‑139 (MQ=255)
                                                                                   aTTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGATCGGCTTTAACTGCTGTCGCAGTCCGTTGc                                                 >  1:349327/1‑139 (MQ=255)
                                                                                       aaGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGATCGGCTTTAACTGCTGTCGCAGTCCGTTGCTACg                                             <  2:1315382/139‑1 (MQ=255)
                                                                                       aaGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGATCGGCTTTAACTGCTGTCGCAGTCCGTTGCTACg                                             <  2:1237873/139‑1 (MQ=255)
                                                                                       aaGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGCGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGATCGGCTTTAACTGCTGTCGCAGTCCGTTGCTACg                                             <  1:622299/139‑1 (MQ=255)
                                                                                       aaGTCGACATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGATCGGCTTTAACTGCTGTCGCAGTCCGTTGCTACg                                             <  2:1114601/139‑1 (MQ=255)
                                                                                         gTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGATCGGCTTTAACTGCTGTCGCAGTCCGTTGCTACGCa                                           <  2:834082/139‑1 (MQ=255)
                                                                                         gTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGATCGGCTTTAACTGCTGTCGCAGTCCGTTGCTACGCa                                           <  1:674618/139‑1 (MQ=255)
                                                                                            ggCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGATCGGCTTTAACTGCTGTCGCAGTCCGTTGCTACGCATgg                                        <  1:1341314/139‑1 (MQ=255)
                                                                                             gCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATtg                                                                                     <  1:459114/93‑1 (MQ=255)
                                                                                             gCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATtg                                                                                     >  2:459114/1‑93 (MQ=255)
                                                                                             gCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGATCGGCTTTAACTGCTGTCGCAGTCCGTTg                                                  <  1:916863/128‑1 (MQ=255)
                                                                                             gCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATATCGGCGTTGGTTTATTGGTTTTAATTGTGGTGATCGGCTTTAACTGCTGTCGCAGTCCGTTg                                                  >  2:916863/1‑128 (MQ=255)
                                                                                                       ttttGGTGGAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGATCGGCTTTAACTGCTGTCGCAGTCCGTTGCTACGCAAGGGAGGGATCGcc                             <  2:790237/139‑1 (MQ=255)
                                                                                                           ggtggAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGt                                                                                             <  2:1215281/71‑1 (MQ=255)
                                                                                                           ggtggAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGt                                                                                             >  1:1215281/1‑71 (MQ=255)
                                                                                                             tggAGGATTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGATCGGCTTTAACTGCTGTCGCAGTCCGTTGCTACGCATGGGAGGGATCGCCATTggg                       >  2:896049/1‑139 (MQ=255)
                                                                                                             cggaggaTTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGCTGGTTTATTGGTTTTACTTGTGGTGATCGGCTTTAACTGCTGTCGCAGTCCGTTGCTACGCATGGGAGGGATCGCCATTggg                       >  2:890956/2‑139 (MQ=255)
                                                                                                               gaggaTTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTaa                                                                                        <  2:1240586/72‑1 (MQ=255)
                                                                                                               gaggaTTTGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTaa                                                                                        >  1:1240586/1‑72 (MQ=255)
                                                                                                                      tGCAGCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGATCGGCTTTAACTGCTGTCGCAGTCCGTTGCTACGCATGGGAGGGATCGCCATTGGGCTATGTGTc              <  1:273871/139‑1 (MQ=255)
                                                                                                                          gCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTAtt                                                                                                >  2:999037/1‑53 (MQ=255)
                                                                                                                          gCCAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTAtt                                                                                                <  1:999037/53‑1 (MQ=255)
                                                                                                                           ccAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGATCGGCTTTAACTGCTGTCGCAGTCCGTTGCTACGCATgg                                        <  2:612670/108‑1 (MQ=255)
                                                                                                                           ccAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGATCGGCTTTAACTGCTGTCGCAGTCCGTTGCTACGCATgg                                        >  1:612670/1‑108 (MQ=255)
                                                                                                                            cAAAAGCAGTGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGATTTAATTGTGGTGATCGGCTTTAACTGCTGTCGCAGTCCGTTGCTACGAATGGGAGGGATCGCCATTGGGCTATGTGTCGGCtat        >  2:503965/1‑139 (MQ=255)
                                                                                                                                 gcagtGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGATCGGCTTTAACTGCTGTCGCAGTCCGTTGCTACGCATGGGAGGGATCGCCATTGGGCTATGTGTCGGCTATATCGc   <  2:562636/139‑1 (MQ=255)
                                                                                                                                  cagtGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTACTTGTGGTGATCGGCTTTAACTGCTGTCGCAGTCCGTTGCTACGCATGGGAGGGATCGCCATTGGGCTATGTGTCGGCTATATCGCa  <  1:890956/139‑1 (MQ=255)
                                                                                                                                  cagtGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGATCGGCTTTAACTGCTGTCGCAGTCCGTTGCTACGCATGGGAGGGATCGCCATTGGGCTATGTGTCGGCTATATCGCa  <  1:503965/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                       
CTTCATTTATCTTGCCCTATTTACGTCGGGTTATTACGCCTACCGTCAGCGGTATTGTGGTACTGATGATCGGCTTAAGCCTGATTAAAGTCGGCATTATCGATTTTGGTGGAGGATTTGCAGCCAAAAGCAGCGGTACGTTCGGCAATTACGAACATCTCGGCGTTGGTTTATTGGTTTTAATTGTGGTGATCGGCTTTAACTGCTGTCGCAGTCCGTTGCTACGCATGGGAGGGATCGCCATTGGGCTATGTGTCGGCTATATCGCA  >  NC_000913/3024754‑3025022

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