Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,613,178 C→G 5.3% N70K (AAC→AAG acpT → 4'‑phosphopantetheinyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,613,1780CG 5.3% 422.8 / 8.2 126N70K (AAC→AAGacpT4'‑phosphopantetheinyl transferase
Reads supporting (aligned to +/- strand):  ref base C (53/66);  new base G (5/2);  total (58/68)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 129 total aligned reads displayed.

TTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCGCTGGCTGGCGAACGCCGTATTCAGCCTCGGGGAACACGCTGAGATGGACGCCG  >  NC_000913/3613029‑3613326
                                                                                                                                                     |                                                                                                                                                      
ttACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAAc                                                                                                                                                        <  2:1509054/150‑1 (MQ=255)
ttACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAAc                                                                                                                                                        <  1:28962/150‑1 (MQ=255)
ttACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCGGGCGTGGCGTGCATTGCTTTCGCACACGGTTTCCCCGCTACCGGAGATCATCTATGGCGAACATGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAAc                                                                                                                                                        <  1:1105898/150‑1 (MQ=255)
 tACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACt                                                                                                                                                       <  1:1990380/150‑1 (MQ=255)
  aCGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAAGtt                                                                                                                                                      >  2:2058704/1‑150 (MQ=255)
     cgagacagcaCCGCTAGGTCCACGACGCTAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAg                                                                                                                                                   <  1:423235/144‑1 (MQ=255)
     cgAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAg                                                                                                                                                   <  1:861033/150‑1 (MQ=255)
     cgAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAg                                                                                                                                                   <  1:798955/150‑1 (MQ=255)
      gAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGATATGGTTCAACTTAAGc                                                                                                                                                  <  2:1344396/150‑1 (MQ=255)
          aagcaCCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCTGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATa                                                                                                                                              <  1:948910/150‑1 (MQ=255)
            gcaCCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGc                                                                                                                                            <  2:743346/150‑1 (MQ=255)
            gcaCCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGc                                                                                                                                            >  1:1313591/1‑150 (MQ=255)
             caCCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCg                                                                                                                                           >  1:1019177/1‑150 (MQ=255)
              aCCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCgg                                                                                                                                          <  2:446194/150‑1 (MQ=255)
              aCCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGCTTAAGTTGAACCATAGCgg                                                                                                                                          <  2:649490/150‑1 (MQ=255)
              aCCGCAAGGAGCACGACGCGAACGCGGGCTTGCGGGGCGTGCATTGCTTTCGCACACGCGTTCCCCGCTACCGGAGATCATCTATGGCGAACACGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGCTTAAGTTGAACCATAGCgg                                                                                                                                          <  1:152550/150‑1 (MQ=255)
               ccGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGt                                                                                                                                         <  1:2052876/150‑1 (MQ=255)
                  cAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGAc                                                                                                                                      <  2:1492679/150‑1 (MQ=255)
                    aGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGa                                                                                                                                    <  2:282682/150‑1 (MQ=255)
                       tCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGatat                                                                                                                                 <  2:390088/150‑1 (MQ=255)
                       tCCACGACGCGAACGCTGGCTGGCGGGGCGTGAATTGCTTTCGCACACGCTTTCACAGCTACCGGAGATCATCTATGGAGAACAATCCAAACCTGCAATTGCGCCGGAAATGACGCTATGGTTCAACTTAAGCCATAGCGGTGACGaaat                                                                                                                                 >  2:584616/1‑150 (MQ=255)
                       gccACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGatat                                                                                                                                 <  2:860984/149‑1 (MQ=255)
                         cACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAGCCTGCATTTGCGGCGGAAATGCCGCTATGGTTCAAGTTAAGCCCTAGCGGTGACGATATCg                                                                                                                               >  1:2042288/1‑150 (MQ=255)
                         cACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATATATGGAGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCg                                                                                                                               >  2:1205530/1‑150 (MQ=255)
                              cgcgAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCctg                                                                                                                          >  2:1726175/1‑150 (MQ=255)
                              cgcgAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCctg                                                                                                                          <  1:901594/150‑1 (MQ=255)
                                cgAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCctgct                                                                                                                        >  2:1518402/1‑150 (MQ=255)
                                 gAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCctgctg                                                                                                                       <  2:1664452/150‑1 (MQ=255)
                                  aaCGCTGTCTGGCGGGGCGTGCATTGCTATCGCACACGCTTTCCCCGCTACCGGAGCTCATCTATGGCGAACAAGGCAAACCTGCAATTGCGCCGGAAATGCCGCTATGGTACAACTTAAGCCAAAGAGGTGACGATATCGCCCTGCTGt                                                                                                                      >  1:964040/1‑150 (MQ=255)
                                  aaCGCTGGCTGGCGGGGCGTGCATTGCTTTAGCACACGCTTTCACCGCTACCGGAGATCATCTATGGAGAACAAGGCAAACCTGCATTTGCGCAGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCctgctgg                                                                                                                      >  2:1137007/1‑149 (MQ=255)
                                      ctggctggCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAg                                                                                                                  <  2:2215525/150‑1 (MQ=255)
                                      ctggctggCGGGGCGTGCATTGCTTTCGCACAAGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAg                                                                                                                  <  1:1526677/150‑1 (MQ=255)
                                         gcgggcGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGtga                                                                                                               <  1:113717/147‑1 (MQ=255)
                                          ctggcGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAATCCTGCATTTGCGCCGGAAATGCCGCTATGGTGCAACTTAAGCCATAGCGGTGACGATATCGCACTGCTGTTTAGtgat                                                                                                              <  2:1778069/150‑1 (MQ=255)
                                            ggcGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACGGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGtgatga                                                                                                            <  1:372487/150‑1 (MQ=255)
                                            ggcGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGtgatgg                                                                                                            <  2:816958/150‑2 (MQ=255)
                                              cGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGATGAGTGATGAAg                                                                                                          >  2:1224127/1‑150 (MQ=255)
                                              cGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCCACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAg                                                                                                          >  1:1804679/1‑150 (MQ=255)
                                                  gCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGa                                                                                                      <  2:1862722/150‑1 (MQ=255)
                                                    gTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAg                                                                                                    <  2:656405/150‑1 (MQ=255)
                                                       cATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCg                                                                                                 <  2:499725/150‑1 (MQ=255)
                                                        aTTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCgg                                                                                                >  2:1827502/1‑150 (MQ=255)
                                                        aTTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCgg                                                                                                <  1:1710749/150‑1 (MQ=255)
                                                            cTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGc                                                                                            >  1:789836/1‑150 (MQ=255)
                                                             tttCGCACACGCTTTCCCCGCGACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTCTGGTTCAACGTCAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGGCGGCTGCg                                                                                           <  1:305084/150‑1 (MQ=255)
                                                                cGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGata                                                                                        >  2:2192635/1‑150 (MQ=255)
                                                                  cacaCGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATc                                                                                      >  2:161463/1‑150 (MQ=255)
                                                                  cacaCGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATc                                                                                      >  1:686136/1‑150 (MQ=255)
                                                                      cGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGGTCAACTTACGCCATCGCGGTGACGATATCGCCCGGCTGTTGAGTGATGAAGGCGAAGTCGCCTGCGATCTCGAAg                                                                                  >  1:59562/1‑150 (MQ=255)
                                                                                  tACCGGAGATCATCTATGGCGAACAAGGCAACCCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCcgc                                                                      <  1:2154070/150‑1 (MQ=255)
                                                                                   aCCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCTGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCcgcc                                                                     >  2:1912388/1‑149 (MQ=255)
                                                                                   aCCGGAGATCATCTAAGGCGAACAACGCAACCCAGCATATGCGCCGGAAATGCAGCTATGGTTCAACTTTAGCCATAGGGGTGACGCTATCTCCCTGCTGTTGAGTGATGAACGCGAAGTCGTCTGCGATATCGAAGTGATTCGCCcgc‑‑g                                                                   <  2:641403/150‑1 (MQ=255)
                                                                                    ccGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCcgc‑‑gc                                                                  >  1:720878/1‑150 (MQ=255)
                                                                                     cGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCcgc‑‑gcg                                                                 <  1:22335/150‑1 (MQ=255)
                                                                                      ggAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGACGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCgc‑‑gcgc                                                                >  1:261478/1‑150 (MQ=255)
                                                                                      ggAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCgc‑‑gcgc                                                                >  2:1802137/1‑150 (MQ=255)
                                                                                       gagaTCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGcc                                                               >  1:1281564/1‑150 (MQ=255)
                                                                                            catcTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTg                                                          >  1:1084199/1‑150 (MQ=255)
                                                                                              tcTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGgc                                                        <  1:433642/150‑1 (MQ=255)
                                                                                               cTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGgcg                                                       <  2:896435/150‑1 (MQ=255)
                                                                                               cTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGgcg                                                       <  2:1743880/150‑1 (MQ=255)
                                                                                               cTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGAGGTGACGATATCTCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGgcg                                                       >  1:1496483/1‑150 (MQ=255)
                                                                                                tATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGgcgc                                                      >  2:453469/1‑150 (MQ=255)
                                                                                                 aTGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCgct                                                     >  2:810592/1‑150 (MQ=255)
                                                                                                  tGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCATCTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCgctg                                                    >  1:265556/1‑150 (MQ=255)
                                                                                                  tGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCgctg                                                    <  2:831222/150‑1 (MQ=255)
                                                                                                    gCGAACAAGGCACACCTGCATTTGCGCCGGAAATGCCGCTATGCTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAGTTGgcgcgcgg                                                  <  1:470740/150‑5 (MQ=255)
                                                                                                    gCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCgctggc                                                  <  2:1313591/150‑1 (MQ=255)
                                                                                                       aaCAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCGctggctgg                                               <  1:1900265/150‑1 (MQ=255)
                                                                                                        aCAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCGCtggctggc                                              <  2:1661472/150‑1 (MQ=255)
                                                                                                        aCAAGGCAAACCTGAATTTGCTCAGGAAATGCCTCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCAGCCTGCGCCAGTTGgcgcgcgggc                                                <  2:298791/150‑7 (MQ=255)
                                                                                                           aGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCGCTGGCTGGCGaa                                           >  1:1394928/1‑150 (MQ=255)
                                                                                                           aGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCGCTGGCTGGCGaa                                           >  1:1615597/1‑150 (MQ=255)
                                                                                                             gCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCGCTGGCTGGCGAACg                                         <  1:1827502/150‑1 (MQ=255)
                                                                                                              cAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCTATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCGCTGGCTGGCGAACGc                                        >  1:834279/1‑150 (MQ=255)
                                                                                                               aaaCCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCGCTGGCTGGCGAACGcc                                       <  2:792416/150‑1 (MQ=255)
                                                                                                                aaCCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCGCTGGCTGGCGAACGCCg                                      <  2:1238135/150‑1 (MQ=255)
                                                                                                                aaCCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCGCTGGCTGGCGAACGCCg                                      <  2:2027108/150‑1 (MQ=255)
                                                                                                                 aCCTCCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCGCTGGCTGGCGAACGCCGt                                     <  2:1726868/150‑1 (MQ=255)
                                                                                                                  ccTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCGCTGGCTGGCGAACGCCGTa                                    >  2:2160272/1‑150 (MQ=255)
                                                                                                                   cTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCGCTGGCTGGCGAACGCCGCAt                                   >  1:1712651/1‑150 (MQ=255)
                                                                                                                     gCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCGCTGGCTGGCGAACGCCGTATTc                                 <  2:536379/150‑1 (MQ=255)
                                                                                                                       aTTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCGCTGGCTGGCGAACGCCGTATTCac                               >  2:1598142/1‑149 (MQ=255)
                                                                                                                            cgcCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCGCTGGCTGGCGAACGCCGTATTCAGCCTCg                          >  2:1806457/1‑150 (MQ=255)
                                                                                                                             gcCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCGCTGGCTGGCGAACGCCGTATTCAGCCTCgg                         >  1:124542/1‑150 (MQ=255)
                                                                                                                             gcCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCGCTGGCTGGCGAACGCCGTATTCAGCCTCgg                         >  1:61638/1‑150 (MQ=255)
                                                                                                                                ggAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCGCTGGCTGGCGAACGCCGTATTCAGCCTCGGGGa                      >  2:1000808/1‑150 (MQ=255)
                                                                                                                                    aTGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCGCTGGCTGGCGAACGCCGTATTCAGCCTCGGGGAacac                  <  2:1909433/150‑1 (MQ=255)
                                                                                                                                     tGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCCACTGGCGCTGGCTGGCGCCCGCCGTATGCCGCCGCGGGGGCCACg                 >  2:1382004/1‑150 (MQ=255)
                                                                                                                                       ccGCTATGGCTTAAGTTGAACCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCGCTGGCTGGCGAACGCCGTATTCAGCCTCGGGGAACACGCt               >  1:2103618/1‑150 (MQ=255)
                                                                                                                                       ccGCTATGGCTTAAGTTGAACCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCGCTGGCTGGCGAACGCCGTATTCAGCCTCGGGGAACACGCt               >  1:129350/1‑150 (MQ=255)
                                                                                                                                        cgcgATGGTTCAACTTCAGCCATAGCGGTGACGATATCGAGCTGCTGTTGAGTGATGAAGGCGAAGTCGTCTTCGATATCGAAGTGATTCGCCAGC‑‑GCGCCAACTGGCGCTGGCTGGCGAACGCCGGATTCAGCCTCGGGGAACACGCTg              >  2:137201/1‑150 (MQ=255)
                                                                                                                                        cGCTATGGTTCAACTTAAGCCATCGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCGCTGGCTGGCGAACGCCGTATTCAGCCTCGGGGAACACGCTg              <  2:348999/150‑1 (MQ=255)
                                                                                                                                           tATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCGCTGGCTGGCGAACGCCGTATTCAGCCTCGGGGAACACGCTgaga           <  1:1216094/150‑1 (MQ=255)
                                                                                                                                            tttgtgcAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGTTGGCGAAGTCGGCTGCGATATCGAAGTGATTCGGCCGC‑‑GCGCCAACTGGCGCTGGCTGGCGAACGCCGTATTCAGCCTCGGGGTAGACGCTGAGAt          <  2:2205930/144‑1 (MQ=255)
                                                                                                                                                ttCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCGCTGGCTGGCGAACGCCGTATTCAGCCTCGGGGAACACGCTGAGATGGAc      >  1:1911534/1‑150 (MQ=255)
                                                                                                                                                ttCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCGCTGGCTGGCGAACGCCGTATTCAGCCTCGGGGAACACGCTGAGATGGAc      >  2:262352/1‑150 (MQ=255)
                                                                                                                                                 tCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCGCTGGCTGGCGAACGCCGTATTCAGCCTCGGGGAACACGCTGAGATGGACg     <  1:1246516/150‑1 (MQ=255)
                                                                                                                                                   aaCTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCGCTGGCTGGCGAACGCCGTATTCAGCCTCGGGGAACACGCTGAGATGGACGcc   >  2:2047909/1‑150 (MQ=255)
                                                                                                                                                    aCTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCGCTGGCTGGCGAACGCCGTATTCAGCCTCGGGGAACACGCTGAGATGGACGCCg  >  1:1641320/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACCGGAGATCATCTATGGCGAACAAGGCAAACCTGCATTTGCGCCGGAAATGCCGCTATGGTTCAACTTAAGCCATAGCGGTGACGATATCGCCCTGCTGTTGAGTGATGAAGGCGAAGTCGGCTGCGATATCGAAGTGATTCGCCCGC‑‑GCGCCAACTGGCGCTGGCTGGCGAACGCCGTATTCAGCCTCGGGGAACACGCTGAGATGGACGCCG  >  NC_000913/3613029‑3613326

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